{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "During the Christmass holydays, I spent some time to convert my Matlab [FVTool](https://github.com/simulkade/FVTool) to a [Julia package](http://pkg.julialang.org/). The result, which I call [JFVM](https://github.com/simulkade/JFVM.jl) is a bit Matlabesque, but the package works. In this post I'm going to show how it works.\n",
    "\n",
    "## Installation\n",
    "To install the package, you can use `Pkg.add(\"JFVM\")` in Julia, or clone the most recent version (recommended). You need to have Julia installed. JFVM relies on `PyPlot`, `PyCall`, and `Mayavi` for visualization. Here's the procedure.\n",
    "\n",
    "<!-- TEASER_END -->\n",
    "\n",
    "### Linux\n",
    "You need to first install [Python](https://www.python.org/), [Matplotlib](http://matplotlib.org/), and [Mayavi](http://code.enthought.com/projects/mayavi/). To do it on an Ubuntu-based system, try\n",
    "```\n",
    "sudo apt-get install python2.7 python-matplotlib mayavi2\n",
    "```\n",
    "Then go to your `.julia/v0.4` or `.julia/v0.3` (recommended) folder and type\n",
    "```\n",
    "git clone https://github.com/simulkade/JFVM.jl.git\n",
    "```\n",
    "\n",
    "### Windows\n",
    "There used to be a few issues with 3D visualization in windows. The current version however should work fine. This is the workflow if you want to give it a try:\n",
    "\n",
    "  - Download and install Anaconda\n",
    "  - Run anaconda command prompt (as administrator) and install mayavi and wxpython:\n",
    "    * `conda install mayavi`\n",
    "    * `conda install wxpython` (Not necessary if you clone the last version of JFVM)\n",
    "  - Install github for windows\n",
    "  - open github shell, go to `.julia/v0.4` or `.julia/v0.3` and type:\n",
    "    * `git clone https://github.com/simulkade/JFVM.git`  \n",
    "\n",
    "Please let me know if it does not work on your windows machines.\n",
    "\n",
    "## What does it do?\n",
    "It solves a transient convection-diffusion equation, i.e., $$\\alpha\\frac{\\partial\\phi}{\\partial t}+\\nabla.\\left(\\mathbf{u}\\phi\\right)+\\nabla.\\left(-D\\nabla\\phi\\right)+\\beta\\phi=\\gamma$$ with the boundary condition: $$a\\nabla\\phi.\\mathbf{n}+b\\phi=c.$$ Each term, i.e., transient, convection, diffusion, linear source, and constant source terms can be called via a separate function, which makes it possible to use the code for solving various coupled and even non-linear convection-diffusion equations. The main difference between this code and the previous Matlab version is that in this one, nonuniform grids can be defined. The supported coordinates are  \n",
    "  - 1D axisymmetric (radial)\n",
    "  - 2D radial (r, theta)\n",
    "  - 2D Cartesian\n",
    "  - 3D Cartesian\n",
    "  - 2D axisymmetric (cylindrical, r, z)\n",
    "  - 3D cylindrical (r, theta, z)\n",
    "One other thing that I always had in mind is the possibility to use hetrogeneous transfer coefficients, e.g., a diffusion coefficient that is different on each cell. In addition, I wanted to be able to define different boundary conditions on the faces of the boundary cell. The last feature that I constantly miss in other PDE solvers is the Robin boundary condition, which can be defined quite easily in JFVM.  \n",
    "One last feature that I added in later stages was the periodic boundary condition.\n",
    "\n",
    "## A simple example\n",
    "This is perhaps the simplest example: a diffusion equation on a 1D domain, with a Dirichlet boundary on each side. The equation read $$\\nabla (-D \\nabla \\phi) = 0$$ where $\\phi$ is zero on the left and one on the right side of the domain. The formulation in `JFVM` always follows this sequence:\n",
    "  - Create a domain and mesh\n",
    "  - Create a boundary for the domain (default: Neumann)\n",
    "  - Change the boundary condition (optional)\n",
    "  - Define the initial condition (optional, only if you have transient terms)\n",
    "  - Define the transfer coefficients\n",
    "  - Calculate the matrix of coefficients and RHS vector for each transfer term in the PDE\n",
    "  - Call the linear solver\n",
    "  - Call the visualization tool\n",
    "There can be other steps in between depending on the problems."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "using JFVM, PyPlot"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAi0AAAHeCAYAAABXBztYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xu4T3X+9/HnlmSHHxElNLl11+wJJRqltEtiX6Si0I7kOD9K7E7KT5mau4OaaTL9asykJIOdW05lKOd2xiHRwSGdzFBR/do7iRC+6/5jTe7kEPZh7fX9Ph/X5bpca1t7v/a6XPbLWu/PZ6UFQRAgSZJUypWJOoAkSdLhsLRIkqRYsLRIkqRYsLRIkqRYsLRIkqRYsLRIkqRYsLRIkqRYsLRIkqRYsLRIkqRYiE1p2bp1K4MGDaJVq1ZUr16dMmXKcP/99x/2+V9++SXdu3enevXqVKhQgWbNmjFv3rxiTCxJkopSbErLV199xciRI9m1axft27cHIC0t7bDO3blzJ5dddhnz58/niSee4KWXXuKkk04iKyuLvLy84owtSZKKSNmoAxyu0047ja+//hqA/Px8nnnmmcM+99lnn2X16tUsXryYpk2bAnDJJZdw9tlnM2jQIJYsWVIsmSVJUtGJzZ2WHzvSdzxOmTKFX/7yl3sLC8AxxxxD165deeONN9i0aVNRR5QkSUUslqXlSK1atYqGDRvud7xBgwYArF69uqQjSZKkI5QSpaWgoICqVavud/yHY/n5+SUdSZIkHaGUKC2SJCn+YjOIWxjVqlWjoKBgv+M/HKtWrdoBz9u0aZPzLpIkHYWaNWtSs2bNIv2cKVFaGjRowLvvvrvf8ZUrVwJQv379/T62adMmmjRpwsaNG4s9nyRJyeaUU07hzTffLNLikhKlpX379tx000288cYb/PrXvwZg9+7djB07lvPPP5+TTz55v3M2bdrExo0bGTt2LBkZGSUdOWXl5OQwfPjwqGOkFK95yfOalzyvefEKAhgzBp54Alq3huuue48ePbqyadOm1C0tM2fOZNu2bXz77bdAuOrnxRdfBKBt27akp6fTq1cvxowZw7p166hTpw4APXv25KmnnqJjx44MGzaM6tWr8+c//5kPP/yQOXPmHPJrZmRkcO655xbvN6a9qlSp4vUuYV7zkuc1L3le8+KzaxfcfDOMHAn33AO/+x289VbxfK1YlZabbrqJ9evXA+FuuBMnTmTixImkpaXxz3/+k1NPPZVEIkEikdhnL5dy5coxd+5cBg0axC233MJ3331Ho0aNmDlzJs2bN4/q25EkKda++QY6doT58+G556B79+L9erEqLf/85z9/9s8899xzPPfcc/sdr1GjBqNHjy6GVJIkpZ7166FtW/j0U3j1VWjRovi/ZqxKiyRJit6yZdCuHaSnw+LFUFKjn+7TolIlOzs76ggpx2te8rzmJc9rXnSmToXMTKhbF5YuLbnCApYWlTL+w1LyvOYlz2te8rzmhRcE8Pjj0KFD+Fho3jyoUaNkM1haJEnSIe3eDf37w223waBBMGFC+GiopDnTIkmSDurbb6FzZ5g1C55+Gvr0iS6LpUWSJB3Qp5/CFVfAunUwYwa0ahVtHkuLJEnaz1tvhYWlbFlYtAgO8MabEudMiyRJ2sf06dC8OdSsCUuWlI7CApYWSZL0I08+CVddBS1bwmuvhcWltLC0SJIk9uyBnBy45RYYOBAmTYIKFaJOtS9nWiRJSnHbtsH114ePhZ58MnwBYmlkaZEkKYVt2hQO3H7wAbz8MrRpE3Wig7O0SJKUolauDHe3TSTg9dfhnHOiTnRozrRIkpSCXn0VLrwQqlULVwiV9sIClhZJklLOX/8a3mFp3hzy8qB27agTHR5LiyRJKSKRCN8d1Ldv+GvaNKhUKepUh8+ZFkmSUsD27XDDDTB5MgwfDgMGQFpa1KmOjKVFkqQk98UX4YZxK1fClCnh7+PI0iJJUhJbsyacX9mxI9zhtkmTqBMdPWdaJElKUnPnQrNmULEiLF0a78IClhZJkpLSqFGQlQVNm8LChXDqqVEnKjxLiyRJSSSRgCFDoFcv6Nkz3Jq/cuWoUxUNZ1okSUoSO3ZAjx7wwgvw6KNwxx3xWyF0KJYWSZKSwFdfhauCVqyAiRPh2mujTlT0LC2SJMXc+++HK4S2bIH58+H886NOVDycaZEkKcby8uCCC+DYY8MVQslaWMDSIklSbI0dCy1bQqNGsGgR1K0bdaLiZWmRJClmggDuvz/clr9LF5g5E044IepUxc+ZFkmSYmTnTujdO7zL8uCDMHhwcq0QOhRLiyRJMVFQAO3bw5IlMH48ZGdHnahkWVokSYqBjz+GNm0gPz/cnv+ii6JOVPKcaZEkqZRbtChcFRQE4V2WVCwsYGmRJKlUmzABWrSAjAxYvBhOPz3qRNGxtEiSVAoFATz0EFx3Xbi77ezZUK1a1KmiZWmRJKmU2bUrXCE0ZAj89rfwt7/BccdFnSp6DuJKklSKbN4c3lnJy4Pnn4du3aJOVHpYWiRJKiX+9a9whdDnn4ePgzIzo05Uuvh4SJKkUuCNN6Bp03DzuMWLLSwHYmmRJClikyaFJaVevXBJ85lnRp2odLK0SJIUkSCAP/wBOnaEK6+EefOgevWoU5VelhZJkiKwezf06wd33gl33w25uVC+fNSpSjcHcSVJKmFbtkCnTuF2/M8+Cz17Rp0oHiwtkiSVoE8+gbZtYf16mDkTWraMOlF8WFokSSohy5dDu3ZQrlz4PqGzzoo6Ubw40yJJUgl46SW4+GKoXRuWLrWwHA1LiyRJxeyJJ+Dqq6F1a1iwAE46KepE8WRpkSSpmOzZAwMGwMCBcPvt8OKLcPzxUaeKL2daJEkqBlu3QnZ2OGw7YgT07Rt1oviztEiSVMQ++ywcuP3oI5g+HbKyok6UHCwtkiQVoXfeCZc0p6XBwoXQsGHUiZKHMy2SJBWRmTPhoougRo1whZCFpWhZWiRJKgIjRsAVV8Cll0JeHpxyStSJko+lRZKkQtizJ1wZdNNN0L8/TJkCFStGnSo5OdMiSdJR2rYNunYNN47705/C5c0qPpYWSZKOwuefhyuE1qyBqVPD36t4WVokSTpCq1aFK4R27YLXX4dzz406UWpwpkWSpCMwezZceCFUqRKuELKwlBxLiyRJh+mZZ6BNG2jWLNyDpU6dqBOlFkuLJEk/I5GAwYOhTx/o3RtefhkqVYo6VepxpkWSpEPYvh1uvDF82eFjj8Gtt4a73arkWVokSTqIL7+Eq64Kt+afNAnat486UWqztEiSdABr14bzK999B6+9BuedF3UiOdMiSdJPLFgAF1wA6emwZImFpbSwtEiS9CPPPw+tWkGTJvCPf8Bpp0WdSD+wtEiSBAQBDB0K3btDt24wY0a4F4tKD2daJEkpb+dO6NkTxo+Hhx+Gu+5yhVBpZGmRJKW0/Hy4+mpYtgwmTIBOnaJOpIOJzeOhrVu3kpOTQ61atUhPT6dRo0ZMmDDhsM6dM2cOl112GTVq1KBSpUqcffbZ/Pd//zeJRKKYU0uSSrMPPwwHbteuhfnzLSylXWxKS4cOHRgzZgz33Xcfr7zyCueddx7Z2dnk5uYe8rxXXnmFVq1aAfDss88ybdo0LrnkEgYOHMhtt91WEtElSaXQ66/D+edDmTLhCqELLog6kX5OLB4PzZgxgzlz5pCbm0vnzp0ByMzMZP369dx555107tyZMmUO3L/+9re/Ub58eaZPn056ejoALVq04P3332f06NEMHz68xL4PSVLpMH489OgRFpXJk6Fq1agT6XDE4k7LlClTqFSpEh07dtzneI8ePdi4cSNLly496Lnp6ekce+yxlC9ffp/jlStX3ltiJEmpIQjggQegSxe47jqYNcvCEiexKC2rVq0iIyNjv7spDRo0AGD16tUHPffmm28mkUgwYMAANm3axObNmxkzZgxTp07lrrvuKtbckqTS4/vvw7sr994Lv/sdjB4N5cpFnUpHIhaPh/Lz8zn99NP3O1713/U4Pz//oOc2atSImTNncu211/LUU08BcMwxxzBs2DBycnKKJ7AkqVT5+mu45ppws7ixY8M7LYqfWJSWwli4cCFt27bl0ksv5Te/+Q0VKlRg7ty5DBkyhO3bt3PPPfdEHVGSVIzWrYO2bcOXH86eDRdfHHUiHa1YlJZq1aod8G5KQUHB3o8fzMCBA6lbty5Tpkwh7d87BWVmZlKmTBnuu+8+unTpQt26dYsnuCQpUkuWwJVXQuXKsHgxnHFG1IlUGLEoLQ0bNiQ3N5dEIrHPXMvKlSsBqF+//kHPXb16NV26dNlbWH7QpEkTEokEa9euPWRpycnJocpP9nHOzs4mOzv7aL4VSVIJmTgx3I6/cWOYOhVOPDHqRMkpNzd3v+1HNm/eXCxfKxalpX379owcOZIXX3yRTj/a+Wf06NHUqlWLpk2bHvTcOnXqsGzZsv0Kz+LFiwGoXbv2Ib/28OHDOffccwv5HUiSSkoQwKOPwt13Q3Y2jBoFP1lAqiJ0oP/Ir1ixgsaNGxf514pFacnKyuLyyy+nX79+bNmyhXr16pGbm8usWbMYN27c3rsovXr1YsyYMaxbt446deoAcPvtt3PTTTfRrl07/vM//5P09HTmzp3LH//4Ry6//PK9K5AkSfG3axfcfDOMHAn33BOuEvIdQskjFqUFYPLkyQwZMoShQ4dSUFBARkYGL7zwwj53XhKJBIlEgiAI9h7r27cvp5xyCo899hh9+vThu+++o27dutx3333ceuutUXwrkqRi8M030LFjuB3/c8+Fb2tWckkLfvwTXnv9cGtr+fLlPh6SpFJu/fpwhdCnn4Y73LZoEXWi1FZcP0Njc6dFkqQDWbYM2rWD9PRwhVBGRtSJVFxisSOuJEkHMnUqZGZC3bqwdKmFJdlZWiRJsRME8Pjj0KFD+Fho3jyoUSPqVCpulhZJUqzs3g39+8Ntt8GgQTBhQvhoSMnPmRZJUmx8+y107hy+nfnpp6FPn6gTqSRZWiRJsfDpp3DFFeG7hGbMgFatok6kkmZpkSSVem+9FRaWsmVh0SI4xNtblMScaZEklWrTp0Pz5lCzZvgCRAtL6rK0SJJKrSefhKuugpYt4bXXwuKi1GVpkSSVOnv2QE4O3HILDBwIkyZBhQpRp1LUnGmRJJUq27bB9deHj4WefDJ8AaIElhZJUimyaVM4cPvBB/Dyy9CmTdSJVJpYWiRJpcLKleHutokEvP46nHNO1IlU2jjTIkmK3KuvwoUXQrVq4QohC4sOxNIiSYrUX/8a3mFp3hzy8qB27agTqbSytEiSIpFIhO8O6ts3/DVtGlSqFHUqlWbOtEiSStz27XDDDTB5MgwfDgMGQFpa1KlU2llaJEkl6osvwg3jVq6EKVPC30uHw9IiSSoxa9aE8ys7doQ73DZpEnUixYkzLZKkEjF3LjRrBhUrwtKlFhYdOUuLJKnYjRoFWVnQtCksXAinnhp1IsWRpUWSVGwSCRgyBHr1gp49w635K1eOOpXiypkWSVKx2LEDevSAF16ARx+FO+5whZAKx9IiSSpyX30VrgpasQImToRrr406kZKBpUWSVKTefz9cIbRlC8yfD+efH3UiJQtnWiRJRSYvDy64AI49NlwhZGFRUbK0SJKKxNix0LIlNGoEixZB3bpRJ1KysbRIkgolCOD++8Nt+bt0gZkz4YQTok6lZORMiyTpqO3cCb17h3dZHnwQBg92hZCKj6VFknRUCgqgfXtYsgTGj4fs7KgTKdlZWiRJR+zjj6FNG8jPD7fnv+iiqBMpFTjTIkk6IosWhauCgiC8y2JhUUmxtEiSDtuECdCiBWRkwOLFcPrpUSdSKrG0SJJ+VhDAQw/BddeFu9vOng3VqkWdSqnG0iJJOqRdu8IVQkOGwG9/C3/7Gxx3XNSplIocxJUkHdTmzeGdlbw8eP556NYt6kRKZZYWSdIB/etf4Qqhzz8PHwdlZkadSKnOx0OSpP288QY0bRpuHrd4sYVFpYOlRZK0j0mTwpJSr164pPnMM6NOJIUsLZIkIFwh9Ic/QMeOcOWVMG8eVK8edSrp/7O0SJLYvRv69YM774S774bcXChfPupU0r4cxJWkFLdlC3TqFG7H/+yz0LNn1ImkA7O0SFIK++QTaNsW1q+HmTOhZcuoE0kHZ2mRpBS1fDm0awflyoXvEzrrrKgTSYfmTIskpaCXXoKLL4batWHpUguL4sHSIkkp5okn4OqroXVrWLAATjop6kTS4bG0SFKK2LMHBgyAgQPh9tvhxRfh+OOjTiUdPmdaJCkFbN0K2dnhsO2IEdC3b9SJpCNnaZGkJPfZZ+HA7UcfwfTpkJUVdSLp6FhaJCmJvfNOuKQ5LQ0WLoSGDaNOJB09Z1okKUnNnAkXXQQ1aoQrhCwsijtLiyQloREj4Ior4NJLIS8PTjkl6kRS4VlaJCmJ7NkTrgy66Sbo3x+mTIGKFaNOJRUNZ1okKUls2wZdu4Ybx/3pT+HyZimZWFokKQl8/nm4QmjNGpg6Nfy9lGwsLZIUc6tWhSuEdu2C11+Hc8+NOpFUPJxpkaQYmz0bLrwQqlQJVwhZWJTMLC2SFFPPPANt2kCzZuEeLHXqRJ1IKl6WFkmKmUQCBg+GPn2gd294+WWoVCnqVFLxc6ZFkmJk+3a48cbwZYePPQa33hrudiulAkuLJMXEl1/CVVeFW/NPmgTt20edSCpZlhZJioG1a8P5le++g9deg/POizqRVPKcaZGkUm7BArjgAkhPhyVLLCxKXZYWSSrFnn8eWrWCJk3gH/+A006LOpEUHUuLJJVCQQBDh0L37tCtG8yYEe7FIqUyZ1okqZTZuRN69oTx4+Hhh+Guu1whJIGlRZJKlfx8uPpqWLYMJkyATp2iTiSVHrF5PLR161ZycnKoVasW6enpNGrUiAkTJhz2+dOmTSMzM5PKlStTsWJF6tevz8iRI4sxsSQdmQ8/DAdu166F+fMtLNJPxeZOS4cOHXjzzTd55JFHOOOMMxg3bhzZ2dkkEgmys7MPee6wYcO455576NevH0OGDOHYY4/lvffeY9euXSWUXpIO7fXXwzss1auHK4Tq1Ys6kVT6xKK0zJgxgzlz5pCbm0vnzp0ByMzMZP369dx555107tyZMmUOfNNo+fLl3HPPPQwbNow77rhj7/FLL720RLJL0s8ZPx569AjvskyeDFWrRp1IKp1i8XhoypQpVKpUiY4dO+5zvEePHmzcuJGlS5ce9Nwnn3yS8uXLc8sttxR3TEk6IkEADzwAXbrAddfBrFkWFulQYlFaVq1aRUZGxn53Uxo0aADA6tWrD3puXl4eGRkZTJw4kTPPPJOyZctSp04dBg8e7OMhSZH5/vvw7sq998LvfgejR0O5clGnkkq3WDweys/P5/TTT9/veNV//5ckPz//oOd+9tlnfPXVVwwcOJAHHniAX/3qV8yZM4dhw4bxySefMHbs2GLLLUkH8vXXcM014WZxY8eGd1ok/bxYlJbCSCQSfPvtt7zwwgt0+vcofmZmJtu2bWP48OHcf//91HPiTVIJWbcO2rYNX344ezZcfHHUiaT4iEVpqVat2gHvphQUFOz9+KHO/fLLL2nduvU+x7Oyshg+fDhvv/32IUtLTk4OVX6yDWV2dvbPrliSpJ9asgSuvBIqV4bFi+GMM6JOJBVebm4uubm5+xzbvHlzsXytWJSWhg0bkpubSyKR2GeuZeXKlQDUr1//oOeeffbZzJo166AfT/uZbSaHDx/Oueeee4SJJWlfEyeG2/E3bgxTp8KJJ0adSCoaB/qP/IoVK2jcuHGRf61YDOK2b9+erVu38uKLL+5zfPTo0dSqVYumTZse9Nxrr70WCJdN/9jf//53jjnmGM7zdamSilEQwCOPhBvFtW8Pc+ZYWKSjFYs7LVlZWVx++eX069ePLVu2UK9ePXJzc5k1axbjxo3be7ekV69ejBkzhnXr1lGnTh0Aunfvzl/+8hduuukmvvrqKzIyMpgzZw5//vOf6dev394/J0lFbdcuuPlmGDkS7rknXCXkO4SkoxeL0gIwefJkhgwZwtChQykoKCAjI2Of4VoIh24TiQRBEOw9VrZsWWbPns1//dd/8dBDD1FQUMD/+l//i0ceeYTbbrstim9FUgr45hvo2DHcjv+558K3NUsqnLTgxz/htdcPz+OWL1/uTIukI7J+fbhC6NNPwx1uW7SIOpFUsorrZ2hs7rRIUhwsWwbt2kF6erhCKCMj6kRS8ojFIK4kxcHUqZCZCXXrwtKlFhapqFlaJKmQggAefxw6dAgfC82bBzVqRJ1KSj6WFkkqhN27oX9/uO02GDQIJkwIHw1JKnrOtEjSUfr2W+jcOXw789NPQ58+USeSkpulRZKOwqefwhVXhO8SmjEDWrWKOpGU/CwtknSE3norLCxly8KiRXCIN4lIKkLOtEjSEZg+HZo3h5o1wxcgWlikkmNpkaTD9OSTcNVV0LIlvPZaWFwklRxLiyT9jD17ICcHbrkFBg6ESZOgQoWoU0mpx5kWSTqEbdvg+uvDx0JPPhm+AFFSNCwtknQQmzaFA7cffAAvvwxt2kSdSEptlhZJOoCVK8PdbRMJeP11OOecqBNJcqZFkn7i1VfhwguhWrVwhZCFRSodLC2S9CN//Wt4h6V5c8jLg9q1o04k6QeWFkkifAw0aBD07Rv+mjYNKlWKOpWkH3OmRVLK274dbrgBJk+G4cNhwABIS4s6laSfsrRISmlffBFuGLdyJUyZEv5eUulkaZGUstasCedXduwId7ht0iTqRJIOxZkWSSlp7lxo1gwqVoSlSy0sUhxYWiSlnFGjICsLmjaFhQvh1FOjTiTpcFhaJKWMRAKGDIFevaBnz3Br/sqVo04l6XA50yIpJezYAT16wAsvwKOPwh13uEJIihtLi6Sk99VX4aqgFStg4kS49tqoE0k6GpYWSUnt/ffDFUJbtsD8+XD++VEnknS0nGmRlLTy8uCCC+DYY8MVQhYWKd4sLZKS0tix0LIlNGoEixZB3bpRJ5JUWJYWSUklCOD++8Nt+bt0gZkz4YQTok4lqSg40yIpaezcCb17h3dZHnwQBg92hZCUTIqstOzatYuRI0fy6quvUq5cOS6++GJ69OhBxYoVi+pLSNJBFRRA+/awZAmMHw/Z2VEnklTUiuTx0Pfff0/Lli257777+J//+R/eeustBg4cSP369XnjjTeK4ktI0kF9/HE4cLt6dbg9v4VFSk5FUlqGDRtGkyZN2LhxI4sWLeKjjz7in//8J1dddRVZWVm8/fbbRfFlJGk/ixaFq4KCILzLctFFUSeSVFyO6PFQt27dOPnkk2nevDkXXXQRJ/x7um358uVMmzZtnz/7i1/8gj/96U+0bt2abt26sWLFCsqWdYRGUtGZMAFuvBF+/WuYMgWqVYs6kaTidEQtYsyYMaxfv55FixYxdOhQCgoKqFy5Mtu3b+eTTz6hTp06+53Tpk0b3n//fUaPHk3v3r2LLLik1BUE8PDD4XuEunSBZ5+F446LOpWk4nbEtz5+8Ytf8Itf/ILsfz80/vbbb7nyyit54okn+Ne//kX16tVp0aIFl156KdX+/d+e/v37c8MNN1haJBXarl3Qt2/4pubf/jb85QohKTUU+nlNpUqVOPHEE/n9738PwFdffcX8+fO59957+Z//+R9+8Ytf0KJFC4455phCh5WU2jZvDt8blJcHzz8P3bpFnUhSSSqSIZOzzjqLl156iSuvvJITTzyRjh070rFjRwA2bNjA3Llzeeedd7j++uv53//7f3P55ZdzwQUXWGQkHbZ//QvatIHPP4fZsyEzM+pEkkpakawe6t+/P4MGDeKTTz7Z72OnnnoqPXr04Fe/+hXjx4+nS5cuvPvuu/Tr168ovrSkFPDGG9C0abh53OLFFhYpVRXJnZYTTzyRe++9l8aNG/OXv/yFDh067PPxadOmUaVKFQDOOOMMzjjjjKL4spJSwKRJ0LVr+A6hadOgevWoE0mKSpGtQe7SpQv5+fl07tyZU089lczMTGrWrMnq1at58803efPNN4vqS0lKAUEAjz0GgwZBx47hDEv58lGnkhSlIn1h4oABA1i4cCFnnXUWU6ZMYcSIEaSnp7N48WJOPvnkovxSkpLY7t3Qrx/ceSfcfTfk5lpYJBXDCxObNm3KSy+9VNSfVlKK2LIFOnUKt+N/9lno2TPqRJJKC7eolVRqfPIJtG0L69fDzJnQsmXUiSSVJpYWSaXC8uXQrh2UKxe+T+iss6JOJKm0KdKZFkk6Gi+9BBdfDLVrw9KlFhZJB2ZpkRSpJ56Aq6+G1q1hwQI46aSoE0kqrSwtkiKxZw8MGAADB8Ltt8OLL8Lxx0edSlJp5kyLpBK3dStkZ4fDtiNGhC9AlKSfY2mRVKI++ywcuP3oI5g+HbKyok4kKS4sLZJKzDvvhEua09Jg4UJo2DDqRJLixJkWSSVi5ky46CKoUSNcIWRhkXSkLC2Sit2IEXDFFXDppZCXB6ecEnUiSXFkaZFUbPbsCVcG3XQT9O8PU6ZAxYpRp5IUV860SCoW27ZB167hxnF/+lO4vFmSCsPSIqnIff55uEJozRqYOjX8vSQVlqVFUpFatSpcIbRrF7z+Opx7btSJJCULZ1okFZnZs+HCC6FKlXCFkIVFUlGytEgqEs88A23aQLNm4R4sdepEnUhSsrG0SCqURAIGD4Y+faB3b3j5ZahUKepUkpKRMy2Sjtr27XDjjeHLDh97DG69NdztVpKKg6VF0lH58ku46qpwa/5Jk6B9+6gTSUp2lhZJR2zt2nB+5bvv4LXX4Lzzok4kKRU40yLpiCxYABdcAOnpsGSJhUVSybG0SDpszz8PrVpBkybwj3/AaadFnUhSKrG0SPpZQQBDh0L37tCtG8yYEe7FIkklyZkWSYe0cyf07Anjx8PDD8Ndd7lCSFI0LC2SDio/H66+GpYtgwkToFOnqBNJSmWxeTy0detWcnJyqFWrFunp6TRq1IgJEyYc8ee55557KFOmDA0aNCiGlFLy+PDDcOB27VqYP9/CIil6sbnT0qFDB958802Y+ajGAAAgAElEQVQeeeQRzjjjDMaNG0d2djaJRILs7OzD+hxvv/02jz32GCeddBJp3t+WDur118M7LNWrhyuE6tWLOpEkxaS0zJgxgzlz5pCbm0vnzp0ByMzMZP369dx555107tyZMmUOfdNo9+7d9OjRg759+/L222+Tn59fEtGl2Bk/Hnr0CO+yTJ4MVatGnUiSQrF4PDRlyhQqVapEx44d9zneo0cPNm7cyNKlS3/2cwwbNozNmzfzwAMPEARBcUWVYisI4IEHoEsXuO46mDXLwiKpdIlFaVm1ahUZGRn73U35YS5l9erVhzx/zZo1PPjgg4wYMYIKFSoUW04prr7/Pry7cu+98LvfwejRUK5c1KkkaV+xeDyUn5/P6aefvt/xqv/+b+ChHvXs2bOHnj17cs0115CVlVVsGaW4+vpruOaacLO4sWPDOy2SVBrForQUxuOPP87HH3/M9OnTo44ilTrr1kHbtuHLD2fPhosvjjqRJB1cLEpLtWrVDng3paCgYO/HD2TDhg0MHTqURx99lLJly7J582YgHMrds2cP33zzDccddxzly5cvvvBSKbVkCVx5JVSuDIsXwxlnRJ1Ikg4tFqWlYcOG5Obmkkgk9plrWblyJQD169c/4Hnr1q1jx44dDBgwgAEDBuz38RNOOIGcnBz++Mc/HvRr5+TkUOUn+5VnZ2cf9jJrqTSaODHcjr9xY5g6FU48MepEkuIqNzeX3NzcfY79cJOgqKUFMVhK88orr9CmTRteeOEFOv1oh6usrCxWr17Nhg0bDrjvyjfffMM777yzz7EgCMjJyWHLli0899xz1KpVi3oH2IRixYoVNG7cmOXLl3PuuecW/TclRSAI4NFH4e67ITsbRo0CbzRKKmrF9TM0FndasrKyuPzyy+nXrx9btmyhXr165ObmMmvWLMaNG7e3sPTq1YsxY8awbt066tSpQ+XKlbn4AA/pK1euzO7duw/4MSlZ7doFN98MI0fCPfeEq4TcY1FSnMSitABMnjyZIUOGMHToUAoKCsjIyNjvzksikSCRSPzsPixpaWnuiKuU8s030LFjuB3/c8+Fb2uWpLiJxeOhKPh4SMli/fpwhdCnn4Y73LZoEXUiSckupR8PSTo6y5ZBu3aQnh6uEMrIiDqRJB29WOyIK+nITZ0KmZlQty4sXWphkRR/lhYpyQQBPP44dOgQPhaaNw9q1Ig6lSQVnqVFSiK7d0P//nDbbTBoEEyYED4akqRk4EyLlCS+/RY6dw7fzvz009CnT9SJJKloWVqkJPDpp3DFFeG7hGbMgFatok4kSUXP0iLF3FtvhYWlbFlYtAgO8lYLSYo9Z1qkGJs+HZo3h5o1wxcgWlgkJTNLixRTTz4JV10FLVvCa6+FxUWSkpmlRYqZPXsgJwduuQUGDoRJk6BChahTSVLxc6ZFipFt2+D668PHQk8+Gb4AUZJShaVFiolNm8KB2w8+gJdfhjZtok4kSSXL0iLFwMqV4e62iQS8/jqcc07UiSSp5DnTIpVyr74KF14I1aqFK4QsLJJSlaVFKsX++tfwDkvz5pCXB7VrR51IkqJjaZFKoUQifHdQ377hr2nToFKlqFNJUrScaZFKme3b4YYbYPJkGD4cBgyAtLSoU0lS9CwtUinyxRfhhnErV8KUKeHvJUkhS4tUSqxZE86v7NgR7nDbpEnUiSSpdHGmRSoF5s6FZs2gYkVYutTCIkkHYmmRIjZqFGRlQdOmsHAhnHpq1IkkqXSytEgRSSRgyBDo1Qt69gy35q9cOepUklR6OdMiRWDHDujRA154AR59FO64wxVCkvRzLC1SCfvqq3BV0IoVMHEiXHtt1IkkKR4sLVIJev/9cIXQli0wfz6cf37UiSQpPpxpkUpIXh5ccAEce2y4QsjCIklHxtIilYCxY6FlS2jUCBYtgrp1o04kSfFjaZGKURDA/feH2/J36QIzZ8IJJ0SdSpLiyZkWqZjs3Am9e4d3WR58EAYPdoWQJBWGpUUqBgUF0L49LFkC48dDdnbUiSQp/iwtUhH7+GNo0wby88Pt+S+6KOpEkpQcnGmRitCiReGqoCAI77JYWCSp6FhapCIyYQK0aAEZGbB4MZx+etSJJCm5WFqkQgoCeOghuO66cHfb2bOhWrWoU0lS8rG0SIWwa1e4QmjIEPjtb+Fvf4Pjjos6lSQlJwdxpaO0eXN4ZyUvD55/Hrp1izqRJCU3S4t0FP71r3CF0Oefh4+DMjOjTiRJyc/HQ9IReuMNaNo03Dxu8WILiySVFEuLdAQmTQpLSr164ZLmM8+MOpEkpQ5Li3QYggD+8Afo2BGuvBLmzYPq1aNOJUmpxdIi/Yzdu6FfP7jzTrj7bsjNhfLlo04lSanHQVzpELZsgU6dwu34n30WevaMOpEkpS5Li3QQn3wCbdvC+vUwcya0bBl1IklKbZYW6QCWL4d27aBcufB9QmedFXUiSZIzLdJPvPQSXHwx1K4NS5daWCSptLC0SD/yxBNw9dXQujUsWAAnnRR1IknSDywtErBnDwwYAAMHwu23w4svwvHHR51KkvRjzrQo5W3dCtnZ4bDtiBHQt2/UiSRJB2JpUUr77LNw4Pajj2D6dMjKijqRJOlgLC1KWe+8Ey5pTkuDhQuhYcOoE0mSDsWZFqWkmTPhoougRo1whZCFRZJKP0uLUs6IEXDFFXDppZCXB6ecEnUiSdLhsLQoZezZE64Muukm6N8fpkyBihWjTiVJOlzOtCglbNsGXbuGG8f96U/h8mZJUrxYWpT0Pv88XCG0Zg1MnRr+XpIUP5YWJbVVq8IVQrt2weuvw7nnRp1IknS0nGlR0po9Gy68EKpUCVcIWVgkKd4sLUpKzzwDbdpAs2bhHix16kSdSJJUWJYWJZVEAgYPhj59oHdvePllqFQp6lSSpKLgTIuSxvbtcOON4csOH3sMbr013O1WkpQcLC1KCl9+CVddFW7NP2kStG8fdSJJUlGztCj21q4N51e++w5eew3OOy/qRJKk4uBMi2JtwQK44AJIT4clSywskpTMLC2Kreefh1atoEkT+Mc/4LTTok4kSSpOlhbFThDA0KHQvTt06wYzZoR7sUiSkpszLYqVnTuhZ08YPx4efhjuussVQpKUKiwtio38fLj6ali2DCZMgE6dok4kSSpJsXo8tHXrVnJycqhVqxbp6ek0atSICRMm/Ox5kyZNolOnTtStW5fjjz+eunXr0rVrVz766KMSSK2i8OGH4cDt2rUwf76FRZJSUazutHTo0IE333yTRx55hDPOOINx48aRnZ1NIpEgOzv7oOf9/ve/p0aNGgwdOpTTTz+dDRs28NBDD3HuueeyZMkSfvWrX5Xgd6Ej9frr4R2W6tXDFUL16kWdSJIUhdiUlhkzZjBnzhxyc3Pp3LkzAJmZmaxfv54777yTzp07U6bMgW8cvfzyy1SvXn2fYy1atOC0007j8ccfZ+TIkcWeX0dn/Hjo0SO8yzJ5MlStGnUiSVJUYvN4aMqUKVSqVImOHTvuc7xHjx5s3LiRpUuXHvTcnxYWgJo1a1KrVi0+/fTTIs+qwgsCeOAB6NIFrrsOZs2ysEhSqotNaVm1ahUZGRn73U1p0KABAKtXrz6iz7du3To2bNjAWWedVWQZVTS+/z68u3LvvfC738Ho0VCuXNSpJElRi83jofz8fE4//fT9jlf993+/8/PzD/tz7d69m549e1KpUiVuvfXWIsuowvv6a7jmmnCzuLFjwzstkiRBjEpLUUkkEvTq1YtFixYxadIkatWqFXUk/du6ddC2bfjyw9mz4eKLo04kSSpNYvN4qFq1age8m1JQULD34z8nCAL69OnDuHHjGD16NO3atSvynDo6S5bA+efD7t2weLGFRZK0v9jcaWnYsCG5ubkkEol95lpWrlwJQP369Q95fhAE9O7dm9GjRzNq1Ciuv/76w/q6OTk5VPnJHvHZ2dmHXGKtIzNxYrgdf+PGMHUqnHhi1IkkSYcrNzeX3NzcfY5t3ry5eL5YEBMzZ84M0tLSggkTJuxzvHXr1kHt2rWDRCJx0HMTiUTQq1evoEyZMsEzzzxzWF9v+fLlARAsX768ULl1cIlEEAwbFgQQBNnZQbB9e9SJJElFobh+hsbmTktWVhaXX345/fr1Y8uWLdSrV4/c3FxmzZrFuHHjSPv3C2h69erFmDFjWLduHXXq1AFgwIABjBo1ip49e1K/fn2WLFmy9/Med9xxNGrUKJLvKZXt2gU33wwjR8I994SrhHyHkCTpUGJTWgAmT57MkCFDGDp0KAUFBWRkZPDCCy/Q6Ud7uicSCRKJBEEQ7D02ffp00tLSGDVqFKNGjdrnc5522mmsW7euxL4HwTffQMeO4Xb8zz0Xvq1ZkqSfkxb8+Ke79lqxYgWNGzdm+fLlnHvuuVHHSRrr14crhD79NNzhtkWLqBNJkopacf0MjdWdFsXbsmXQrh2kp4crhDIyok4kSYqT2Cx5VrxNnQqZmVC3LixdamGRJB05S4uKVRDA449Dhw7hY6F586BGjahTSZLiyNKiYrN7N/TvD7fdBoMGwYQJ4aMhSZKOhjMtKhbffgudO4dvZ376aejTJ+pEkqS4s7SoyH36KVxxRfguoRkzoFWrqBNJkpKBpUVF6q23wsJStiwsWgQ/83YFSZIOmzMtKjLTp0Pz5lCzZvgCRAuLJKkoWVpUJJ58Eq66Clq2hNdeC4uLJElFydKiQtmzB3Jy4JZbYOBAmDQJKlSIOpUkKRk506Kjtm0bXH99+FjoySfDFyBKklRcLC06Kps2hQO3H3wAL78MbdpEnUiSlOwsLTpiK1eGu9smEvD663DOOVEnkiSlAmdadERefRUuvBCqVQtXCFlYJEklxdKiw/bXv4Z3WJo3h7w8qF076kSSpFRiadHPSiTCdwf17Rv+mjYNKlWKOpUkKdU406JD2r4dbrgBJk+G4cNhwABIS4s6lSQpFVladFBffBFuGLdyJUyZEv5ekqSoWFp0QGvWhPMrO3aEO9w2aRJ1IklSqnOmRfuZOxeaNYOKFWHpUguLJKl0sLRoH6NGQVYWNG0KCxfCqadGnUiSpJClRUC4QmjIEOjVC3r2DLfmr1w56lSSJP1/zrSIHTugRw944QV49FG44w5XCEmSSh9LS4r76qtwVdCKFTBxIlx7bdSJJEk6MEtLCnv//XCF0JYtMH8+nH9+1IkkSTo4Z1pSVF4eXHABHHtsuELIwiJJKu0sLSlo7Fho2RIaNYJFi6Bu3agTSZL08ywtKSQI4P77w235u3SBmTPhhBOiTiVJ0uFxpiVF7NwJvXuHd1kefBAGD3aFkCQpXiwtKaCgANq3hyVLYPx4yM6OOpEkSUfO0pLkPv4Y2rSB/Pxwe/6LLoo6kSRJR8eZliS2aFG4KigIwrssFhZJUpxZWpLUhAnQogVkZMDixXD66VEnkiSpcCwtSSYI4KGH4Lrrwt1tZ8+GatWiTiVJUuFZWpLIrl3hCqEhQ+C3v4W//Q2OOy7qVJIkFQ0HcZPE5s3hnZW8PHj+eejWLepEkiQVLUtLEvjXv8IVQp9/Hj4OysyMOpEkSUXPx0Mx98Yb0LRpuHnc4sUWFklS8rK0xNikSWFJqVcvXNJ85plRJ5IkqfhYWmIoCOAPf4COHeHKK2HePKhePepUkiQVL0tLzOzeDf36wZ13wt13Q24ulC8fdSpJkoqfg7gxsmULdOoUbsf/7LPQs2fUiSRJKjmWlpj45BNo2xbWr4eZM6Fly6gTSZJUsiwtMbB8ObRrB+XKhe8TOuusqBNJklTynGkp5V56CS6+GGrXhqVLLSySpNRlaSnFnngCrr4aWreGBQvgpJOiTiRJUnQsLaXQnj0wYAAMHAi33w4vvgjHHx91KkmSouVMSymzdStkZ4fDtiNGQN++USeSJKl0sLSUIp99Fg7cfvQRTJ8OWVlRJ5IkqfSwtJQS77wTLmlOS4OFC6Fhw6gTSZJUujjTUgrMnAkXXQQ1aoQrhCwskiTtz9ISsREj4Ior4NJLIS8PTjkl6kSSJJVOlpaI7NkTrgy66Sbo3x+mTIGKFaNOJUlS6eVMSwS2bYOuXcON4/70p3B5syRJOjRLSwn7/PNwhdCaNTB1avh7SZL08ywtJWjVqnCF0K5d8PrrcO65USeSJCk+nGkpIbNnw4UXQpUq4QohC4skSUfG0lICnnkG2rSBZs3CPVjq1Ik6kSRJ8WNpKUaJBAweDH36QO/e8PLLUKlS1KkkSYonZ1qKyfbtcOON4csOH3sMbr013O1WkiQdHUtLMfjyS7jqqnBr/kmToH37qBNJkhR/lpYitnZtOL/y3Xfw2mtw3nlRJ5IkKTk401KEFiyACy6A9HRYssTCIklSUbK0FJHnn4dWraBJE/jHP+C006JOJElScrG0FFIQwNCh0L07dOsGM2aEe7FIkqSi5UxLIezcCT17wvjx8PDDcNddrhCSJKm4WFqOUn4+XH01LFsGEyZAp05RJ5IkKbnF5vHQ1q1bycnJoVatWqSnp9OoUSMmTJhwWOd++eWXdO/enerVq1OhQgWaNWvGvHnzjjrLhx+GA7dr18L8+RYWSZJKQmxKS4cOHRgzZgz33Xcfr7zyCueddx7Z2dnk5uYe8rydO3dy2WWXMX/+fJ544gleeuklTjrpJLKyssjLyzviHAsXwvnnQ5ky4QqhCy442u9IkiQdiViUlhkzZjBnzhxGjBhBnz59yMzM5Omnn+byyy/nzjvvJJFIHPTcZ599ltWrV/N//+//JTs7m8suu4wXX3yRM844g0GDBh1RjtxcuOwyaNAAFi2CevUK+53pp36uhKroec1Lnte85HnNk0MsSsuUKVOoVKkSHTt23Od4jx492LhxI0uXLj3kub/85S9p2rTp3mPHHHMMXbt25Y033mDTpk0/+/WDAB54AK6/Hq67DmbNgqpVj/770cH5D0vJ85qXPK95yfOaJ4dYlJZVq1aRkZFBmTL7xm3QoAEAq1evPuS5DRs23O/44ZwLsGtXuELo3nvhd7+D0aOhXLkj/AYkSVKhxWL1UH5+Pqeffvp+x6v++3ZHfn7+Qc8tKCjY++eO9FyA/v3h3Xdh7Fjo0uVIUkuSpKIUi9ISpQ8+gNmz4eKLo04iSVJqi0VpqVat2gHviBQUFOz9+KHO/eHPHem5APfe+x4VK8KKFUeSWEdr8+bNrPBilyivecnzmpc8r3nJeu+994rl88aitDRs2JDc3FwSicQ+cy0rV64EoH79+gc9t0GDBrz77rv7Hf+5c2vWrMkpp5zCnXd2LUx0HYXGjRtHHSHleM1Lnte85HnNS9Ypp5xCzZo1i/RzpgVBEBTpZywGr7zyCm3atOGFF16g0492csvKymL16tVs2LCBtIPsn/+Xv/yFm266iSVLlvDrX/8agN27d3POOefwH//xHyxatOigX3fTpk2HtbpIkiTtq2bNmqlZWgBat27Nm2++ySOPPEK9evXIzc3lmWeeYdy4cWRnZwPQq1cvxowZw7p166hTpw4A33//PY0bN2bLli0MGzaM6tWr8+c//5m///3vzJkzh+bNm0f5bUmSpMMUi8dDAJMnT2bIkCEMHTqUgoICMjIy9rvzkkgkSCQS/LiHlStXjrlz5zJo0CBuueUWvvvuOxo1asTMmTMtLJIkxUhs7rRIkqTUFovN5YpSaXrxYqo42ms+adIkOnXqRN26dTn++OOpW7cuXbt25aOPPiqB1PFWmL/nP3bPPfdQpkyZvZsx6uAKe82nTZtGZmYmlStXpmLFitSvX5+RI0cWY+L4K8w1nzNnDpdddhk1atSgUqVKnH322fz3f//3IV8Lk+q2bt3KoEGDaNWqFdWrV6dMmTLcf//9h31+kfwMDVLM5ZdfHpxwwgnB008/HSxYsCDo06dPkJaWFowfP/6Q5+3YsSOoX79+cOqppwbjx48P5syZE1x99dXBscceG7z22msllD6ejvaaN23aNGjXrl0watSoIC8vLxg7dmzwq1/9KqhUqVKwevXqEkofT0d7zX/srbfeCsqXLx+cfPLJQYMGDYoxbXIozDV/+OGHg2OOOSbo379/8Oqrrwbz5s0LnnrqqeCpp54qgeTxdbTXfObMmUFaWlrQokWL4KWXXgrmzp0bDBgwIEhLSwsGDhxYQunj55///GdQpUqV4JJLLtl7re+///7DOreofoamVGn5+9//HqSlpQUvvPDCPsdbtWoV1KpVK9izZ89Bz33qqaeCtLS0YMmSJXuP7d69OzjrrLOCpk2bFlvmuCvMNf/yyy/3O7Zx48agXLlyQe/evYs8a7IozDX/wa5du4JzzjknyMnJCS655BJLy88ozDV/8803g2OOOSb4/e9/X9wxk0phrvn1118fpKenB999990+x1u3bh1Urly5WPImm6+++uqISktR/QxNqcdDUb94MRUV5ppXr159v2M1a9akVq1afPrpp0WeNVkU5pr/YNiwYWzevJkHHnhgn8F2HVhhrvmTTz5J+fLlueWWW4o7ZlIpzDVPT0/n2GOPpXz58vscr1y5Munp6cWSN9kc6b8LRfUzNKVKS5QvXkxVhbnmB7Ju3To2bNjAWWedVWQZk01hr/maNWt48MEHGTFiBBUqVCi2nMmkMNc8Ly+PjIwMJk6cyJlnnknZsmWpU6cOgwcPZteuXcWaO84Kc81vvvlmEokEAwYMYNOmTWzevJkxY8YwdepU7rrrrmLNnaqK6mdoSpWW/Pz8o355YmFfvJiqCnPNf2r37t307NmTSpUqceuttxZZxmRTmGu+Z88eevbsyTXXXENWVlaxZUw2hbnmn332GR988AEDBw4kJyeHuXPn0r17d/7whz/Qo0ePYsscd4W55j9sezFx4kRq1apF1apV6dWrFw899BA5OTnFljmVFdXP0Njs06LUlkgk6NWrF4sWLWLSpEnUqlUr6khJ6fHHH+fjjz9m+vTpUUdJGYlEgm+//XaffacyMzPZtm0bw4cP5/7776devXoRp0wuCxcupG3btlx66aX85je/oUKFCsydO5chQ4awfft27rnnnqgj6iBSqrRE+eLFVFWYa/6DIAjo06cP48aNY8yYMbRr167IcyaTo73mGzZsYOjQoTz66KOULVuWzZs3A+Edrj179vDNN99w3HHH7TcHoML/2/Lll1/SunXrfY5nZWUxfPhw3n77bUvLARTmmg8cOJC6desyZcqUva+AyczMpEyZMtx333106dKFunXrFk/wFFVUP0NT6vFQw4YNee+99/Zbh1+cL15MdYW55hAWlt69ezN69GieffZZrr/++mLLmiyO9pqvW7eOHTt2MGDAAKpWrbr316JFi3jvvfc44YQT+K//+q9izx9Hhfl7fvbZZx9yqPFg71VLdYW55qtXr6Zx48b7XdsmTZqQSCRYu3Zt0QdOcUX1MzSlSkv79u3ZunUrL7744j7HR48eTa1atfaZaj7QuWvXruWNN97Ye2z37t2MHTuW888/n5NPPrnYcsdZYa75D3dYRo8ezdNPP82NN95Y3HGTwtFe80aNGrFgwYJ9fs2fP5+zzz6bunXrsmDBAm6++eaS+BZipzB/z6+99loAZsyYsc/xv//97xxzzDGcd955RR84CRTmmtepU4dly5btV3gWL14MQO3atYs+cIorsp+hR7AsOym0atUqqFq1ajBy5Mhg3rx5B9yMqGfPnkHZsmWDDRs27D22c+fOfTbGmT17dtC+ffugXLlyQV5eXhTfSmwc7TXv379/kJaWFvTq1StYsmRJsHjx4r2/VqxYEcW3EhtHe80PJDMzM6hfv35xR469o73mu3btCho3bhxUqVIleOKJJ4LZs2cHd911V1C2bNnglltuieJbiY2jveYjRowI0tLSgjZt2gTTpk0LZs2aFdx1113BscceG7Rq1SqKbyU2ZsyYEUycODEYNWpUkJaWFnTq1CmYOHFiMHHixL373hTnz9CUKy1bt24NBg4cGNSsWTM47rjjgnPOOSeYMGHCPn+me/fuQZkyZYL169fvc/yLL74IbrzxxqBatWpBenp60KxZs2Du3LklGT+Wjvaan3baaUGZMmWCtLS0/X7VrVu3pL+NWCnM3/OfcnO5w1OYa15QUBD07ds3OPnkk4Ny5coFv/zlL4PHHnusJOPHUmGu+bRp04KLL744qFGjRlCxYsWgQYMGwYMPPrjfhnPa12mnnbb33+Ef//v842tcnD9DfWGiJEmKhZSaaZEkSfFlaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEkSbFgaZEUC3369KFMmTK0atXqgB//7W9/S5kyZTj77LP5/vvvSzidpJKQFgRBEHUISfo527Zt4+yzz2bdunU8/vjjDBw4cO/Hli5dyoUXXsixxx7LsmXLqF+/foRJJRUX77RIioUKFSowZswYjjnmGAYPHsyaNWuAsMzccMMNJBIJ/s//+T8WFimJWVokxUazZs0YNGgQO3bsoGvXruzatYvbb7+djz76iMzMTO64446oI0oqRj4ekhQru3fvpmnTprz11lu0bNmSOXPmULlyZd59913q1KkTdTxJxcjSIil23nvvPRo3bsyOHTsAeP7557nhhhsiTiWpuFlaJMXO999/T4MGDfjwww+pUqUKn3zyCRUqVIg6lqRi5kyLpNgZMmQIH374IWXKlGHz5s3ceuutUUeSVAIsLZJiJS8vjz/+8Y9UqFCBOXPmUKVKFZ555hmmT58edTRJxczSIik2tmzZQrdu3QiCgN///vdccsklPPXUU0C4+dxXX30VcUJJxcnSIik2BgwYwIYNG8jKyqJv374AZGdn06lTJ7744gt+85vfRJxQUnFyEFdSLEyZMoVrrrmGqlWrsmrVKk4++eS9H/v661IeuuUAAACwSURBVK9p0KABGzduZNSoUXTv3j26oJKKjaVFUqn3xRdf0KBBA/Lz88nNzaVTp077/ZnZs2eTlZVFpUqVePfddzn11FMjSCqpOFlaJElSLDjTIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSIkmSYsHSov/Xbh2QAAAAAAj6/7odga4QABakBQBYkBYAYCF18hklt1P3awAAAABJRU5ErkJggg=="
      ],
      "text/plain": [
       "Figure(PyObject <matplotlib.figure.Figure object at 0x7fa9bcbbd5d0>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "PyObject <matplotlib.text.Text object at 0x7fa9bb423310>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "L= 1.0 # length of the domain\n",
    "Nx= 20 # number of cells\n",
    "m= createMesh1D(Nx, L) # create the domain and mesh\n",
    "BC= createBC(m) # create the boundary condition\n",
    "BC.left.a[:]= 0.0 # left side Neumann term equal to zero\n",
    "BC.left.b[:]= 1.0 # left side Dirichlet term coefficient equal to one\n",
    "BC.left.c[:]= 0.0 # left side Dirichlet term value equal to zero\n",
    "BC.right.a[:]= 0.0\n",
    "BC.right.b[:]= 1.0\n",
    "BC.right.c[:]= 1.0\n",
    "D_val= 1.0 # value of the transfer coefficient\n",
    "D= createCellVariable(m, D_val) # assign D_val to all the cells\n",
    "# harmonic average of the transfer coefficient on the cell faces:\n",
    "D_face= harmonicMean(D)\n",
    "(Mbc, RHSbc)= boundaryConditionTerm(BC)\n",
    "Mdiff= diffusionTerm(D_face)\n",
    "M= Mdiff+Mbc # matrix of coefficients\n",
    "RHS= RHSbc # off course!\n",
    "phi= solveLinearPDE(m, M, RHS) # solve the linear PDE\n",
    "figure(figsize=(6.0,5.0)) # unit is [cm] for me\n",
    "visualizeCells(phi) # visualize the results\n",
    "# decorate the plot\n",
    "xlabel(\"x\", fontsize=16)\n",
    "ylabel(L\"\\phi\", fontsize=16)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## A colorful a bit playful case\n",
    "The other night, my daughter wanted to see what I do with my computer. So I decided to solve an equation with fancy colorful results. So I created a diffusion process on a radial coordinate, Dirichlet boundary condition near the center of the circle (with alternationg zero-one values on the left boundary). I reduced the value of the diffusion coefficient in the direction of increasing $\\theta$ to increase the contrast in the result. Let me recreate it here."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAj4AAAIzCAYAAAAArBNyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzsnXmcT9X/x1+zMrbsQ0YGWeorqWw1yQySb5hfIn1jvkWWQr60IWXpmxBJKUsSEvENLTOVZBkSlUS2iDD2GWtmRraZ+fz+mLnjzp27nHPvucvn83k/H4/Pg/l87j3n3O2c132/3+d9Qnw+nw8EQRAEQRBBQKjbDSAIgiAIgnAKEj4EQRAEQQQNJHwIgiAIgggaSPgQBEEQBBE0kPAhCIIgCCJoIOFDEARBEETQQMKHIAiCIIiggYQPQRAEQRBBAwkfgiAIgiCCBhI+BEHYxubNm/HAAw+gXLlyKFu2LNq0aYNNmzapbrt161a0a9cOZcuWRYUKFdC1a1ccOnSo2HZ79uxBXFwcbrjhBiQmJiIjI8PuwyAIIoAg4UMQhC388ssvuO+++3DlyhUsXLgQH3/8MS5fvoy2bdvip59+KrLt3r17ER8fj5ycHCxduhRz587Fvn370KpVK5w5c6bItt26dUOLFi3wxRdfoGzZsujfv7+Th0UQhJ8TQmt1EQRhBx06dMCOHTtw8OBBlCxZEgCQnZ2NOnXqoH79+vjhhx8Kt+3evTvWr1+PAwcOoEyZMgCAI0eOoF69enj22WcxceJEAMDZs2dRv359nD17FgCQmZmJmJgYZGZmOnx0BEH4K2TxIQjCFjZu3Ij4+PhC0QMAZcqUQatWrbBp0yakp6cDAHJycvDVV1+ha9euhaIHAG666SYkJCTg888/L/yuQoUKCAkJweLFi5GZmYlZs2ahYcOGzh0UQRB+T7jbDSAIIjC5du0aSpQoUex76btdu3ahWrVqOHDgAC5fvozGjRsX2/a2227DqlWrcPXqVURGRiI0NBRz5szB448/juzsbNSoUQNffvml7cdCEETgQMKHIAhbuPXWW/Hjjz/C5/MhJCQEQL51R4rvkdxV0r8VK1YsVkbFihXh8/lw/vx5REdHAwAeeughZGRk4OjRo6hTpw7Cw6kbIwiCHXJ1EQRhC4MHD8a+ffvwzDPP4MSJEzh69CiefvppHD16FAAQGmq++4mKikL9+vVJ9BAEwQ0JH4IgbKF3796YOHEiPv74Y8TExKBWrVr4448/8MILLwAAatSoAQCoVKkSAODcuXPFyjh37hxCQkJQoUIF5xpOEERAQ8KHIAjbGDZsGM6ePYtdu3bh8OHD2LBhA86dO4cyZcrgrrvuAgDUrVsXUVFR2LFjR7H9d+7ciXr16iEyMtLpphMEEaCQ8CEIwlYiIiJw6623ombNmjhy5Aj+97//oV+/foVBzuHh4ejcuTM+++wzZGdnF+535MgRpKam4uGHH3ar6QRBBCCUx4cgCFvYvXs3li5dimbNmqFEiRLYvn07Jk6ciDp16iA1NRWlSpUq3PaPP/5As2bNcOedd2LEiBG4dOkSRo8ejb/++gu//fZboTuMIAjCKiR8CIKwhf3796N///7YuXMnsrOzUatWLfzrX//CiBEjEBUVVWz7rVu3Yvjw4fjxxx8RHh6Otm3b4s0330Tt2rVdaD1BEIEKCR+CIAiCIIIGivEhCIIgCCJoIOFDEARBEETQQMKHIAiCIIiggYQPQRAEQRBBAwkfgiAIgiCCBhI+BEEQBEEEDSR8CIIgCIIIGkj4EARBEAQRNJDwIQiCIAgiaCDhQxAEQRBE0EDChyAIgiCIoIGED0EQBEEQQQMJH4IgCIIgggYSPgRBEARBBA0kfAiCIAiCCBpI+BAEQRAEETSQ8CEIgiAIImgId7sBBEEEB3l5ecjOzkZ2djYuXryI3Nxc5OXlIS8vr/D/Pp8PoaGhRT5hYWEIDw9HmTJlUKZMGZQqVQohISFuHw5BEH4KCR+CIJjx+Xy4cOECTp48iRMnTuDkyZM4efIk0tPTceHCBWRnZyMrKwtZWVnF/n/x4kUhbQgNDS0UQWXLli38yP8uX748qlevjurVq+PGG28s/H+pUqWEtIEgCP8lxOfz+dxuBEEQ3uDChQvYt28f9u3bh2PHjhWKG7nIuXTpUpF9ypUrh2rVqqF8+fK6YkT+d+nSpREeHl7MugPkiyvJEiR9rl27pimq1P4+f/48Tp48icuXLxdp6w033FBEDEn/1qpVCw0aNEDdunVRsmRJx843QRDOQ8KHIIKMK1eu4MCBA4UCZ9++ffjjjz+wb98+nDp1qnC78uXLFxMIyv9Xr14dpUuXdvFotPH5fPjrr78KBZtcvCkF3d9//w0ACAkJQa1atVC/fn3Ur18fDRo0KPx/zZo1ERYW5vJREQRhFRI+BBGg+Hw+HDhwAL/++iu2bNmCXbt2Yd++fUhLS0NeXh4AoGzZsqqDfP369VG2bFmXj8AZfD4fTp8+XUQASp8///wTV69eBQCUKFECN998Mxo0aIDbb78dTZs2xV133YXo6GiXj4AgCB5I+BBEAKAUOb/++iu2bt2KCxcuAABuuukm3H777UXETYMGDRAdHU2Bwjrk5ubi8OHDRUTR3r17sW3bNpw/fx4AUKNGDdx111246667SAwRhB9Awocg/JCTJ09iw4YN+OWXX1RFjjQQS58qVaq43OLAwufzIS0trVBkSh81MdS8eXPExcUFjQWNILwOCR+C8APS09Oxfv16rFu3Dqmpqfjjjz8AXBc5kqXhzjvvJJHjEpIYkkSQJIrOnz+PsLAwNGvWDPHx8YiPj0dcXBzKlCnjdpMJIigh4UMQHuTUqVNYv349UlNTsW7dOuzZswcA0LBhw8LBs3Xr1qhWrZrLLSX08Pl82LdvH9atW1coWjMyMhAeHo5mzZohISEB8fHxuOeeezwbJE4QgQYJH4LwAH///TfWrFmDlStXIjU1Fb///jsAoH79+oWDY+vWrVG9enWXW0pYwefz4Y8//igUQevWrcOpU6cQERGB5s2bIz4+Hh07dkTz5s1pBhlB2AQJH4JwifT0dHz99ddITk7GqlWrcOnSJdStWxdt2rRBQkICWrdujRtvvNHtZhI24vP5sHfv3kIRlJqaijNnzqBq1aro1KkTOnfujPvvv5+sQQQhEBI+BOEQPp8Pu3fvRnJyMpKTk/Hzzz8jNDQUcXFxSExMROfOndGgQQO3m0m4SG5uLn7++efCe2TPnj0oUaIE2rVrh8TERHTq1InEMEFYhIQPQdjItWvX8P333yM5ORkpKSk4dOgQSpcujQ4dOiAxMREPPvggKleu7HYzCY+yf/9+pKSkIDk5GT/88ANyc3PRtGlTJCYmIjExEY0bN6Z0BATBCQkfghCMz+fDjz/+iIULF+J///sfzp07h5iYmEKrTnx8PC2LQHBz7tw5rFixAsnJyVixYgWysrJQr149JCUloWfPnqhbt67bTSQIv4CED0EIYu/evVi0aBE++eQTHDx4EDExMejRowceffRR3HHHHfRmTgjj6tWrSE1NxeLFi7F8+XJkZ2ejZcuWSEpKQvfu3SmlAUHoQMKHICyQnp6OJUuWYNGiRdiyZQvKlSuHRx55BElJSbjvvvsKF94kCLv4+++/kZKSgoULF+Lbb78FADzwwANISkpCYmIirUhPEApI+BAEJ9nZ2fjiiy+wcOFCrFq1CmFhYejYsSOSkpLQsWNHcmMRrnH69Gl8+umnWLRoEX788UeUKVMGDz/8MJKSktCmTRuaIk8QIOFDEMzs3r0bM2fOxIIFC5CVlYVWrVohKSkJ3bp1Q8WKFd1uHkEU4cCBA1i0aBEWLVqEffv2oVatWnjqqafQp08fVK1a1e3mEYRrkPAhCB2uXr2KL774AjNmzMD69esRHR2Nfv36oU+fPoiNjXW7eQRhiM/nw+bNm/H+++9j8eLFyM3NxSOPPIKBAwfinnvuodgzIugg4UMQKhw7dgyzZ8/GBx98gPT0dNx3330YOHAgunTpgsjISLebRxCmOHfuHObPn4+ZM2fizz//ROPGjTFw4ED07NmT1g4jggYSPgRRgM/nw5o1azBjxgwkJycjKioKjz/+OAYMGIBGjRq53TyCEEZeXh5Wr16NGTNmICUlBaVLl8YTTzyBAQMG4NZbb3W7eQRhKyR8iKDn0qVLmDdvHt555x3s27cP//jHPzBo0CAkJSWhbNmybjfPbwkJedWxuny+MY7VFWgcOXKk0Lp56tQpxMfH4/nnn0fHjh3JDUYEJCR8iKDlwoULmDFjBt5++22cOXMG3bp1wzPPPIN7772XOnwVnBQydkIiSZ2rV6/is88+w7Rp0/Djjz/itttuw4gRI9C9e3eEh4e73TyCEAYJHyLoyMjIwNtvv40ZM2bg8uXL6N27N1588cWgz3wbKMLGKsEujHw+HzZs2IAJEybg22+/RZ06dfDiiy+iV69elKqBCAhI+BBBQ1paGiZPnoy5c+ciPDwcAwYMwLPPPovq1au73TTHIZHDTzAKom3btmHixIlYunQpoqOj8eyzz+Lpp59GuXLl3G4aQZiGhA8R8OzevRtvvPEGPvnkE5QvXx5Dhw7FoEGDUKFCBbebZjveFDhRAK6Z3DdHZEMsEyxiaP/+/Zg0aRI++ugjlC5dGoMGDcKQIUNoaQzCLyHhQwQsO3bswOjRo/Hll18iJiYGL7zwAvr27YvSpUu73TRbcF7kSHEfESb2DQzho0Ygi6Hjx4/jrbfewvvvv4+8vDz0798fI0eOpISIhF9BwocIOA4dOoTRo0dj0aJFqFu3LkaOHImePXsGXP4d+4UOa0ArCR89AlEInT17Fu+99x6mTp2K3NxcPP/883juuefIBUb4BSR8iIDh1KlTGDduHGbNmoVKlSphzJgx6NOnDyIizAzM3sM7QkdJoAof6XyIrTOQhNDZs2fxxhtvYNq0aShbtixeeeUVPP300yhRooTbTSMITUj4EH5PZmYmpkyZgilTpiA8PBzDhw/Hf/7zn4BwaTnrviLhcx2jcyGuHYEghI4dO4axY8di3rx5qFmzJv773/+iZ8+etCgq4UlI+BB+y5UrVzBr1iyMGzcO2dnZGDx4MEaMGOHXC4a6G4xMwuc6Zs6F9bb5uwjau3cvXn75ZXz22Wdo1KgRJkyYQIkQCc9BwofwO3Jzc7Fo0SKMHj0aR48exZNPPokxY8YgJibG7aaZwlszr8wM+E4JHzfcXLzIz8Uly63wVyG0efNmjBgxAqmpqYiLi8PEiRNx7733ut0sggBAwofwM3766ScMGjQIW7duxcMPP4zXX38dDRs2dLtZ3HhL7MgRLXy0yjMjCvxN+CixJoT8TQT5fD589913eOmll7Bt2zY8+uijePPNN/32BYUIHEj4EH7B6dOnMWLECMydOxd33nkn3nvvPdx9991uN4sL74odOWYG/CgT+wSi8OG1fJkXQv4kgvLy8rBo0SK8+OKLyM7OxujRozF06NCAm2VJ+A8kfAhPk5ubi1mzZuGVV15BSEgIXn/9dfTv399vgiatiR17ZhWx1ckDCZ98zM4evAYrx+YvIujChQsYO3Ys3n33XdSrVw/vvvsu2rVr53aziCCEhA/hWTZt2oRBgwZh+/bt6NOnDyZMmIDKlSu73SwmzAkercHW6wM+CZ98rAgfOeaO018E0M6dOzFo0CBs2LAB3bp1w1tvvYWaNWu63SwiiCDhQ3iOU6dOYfjw4Zg/fz6aNm2K6dOno3nz5m43yxCxYkeO1wd8gF/8eFn4ODm7DdAP9A5MEeTz+fDJJ5/ghRdeQGZmJl555RU899xzlP+HcAQSPoRnyMnJwcyZMzFq1CiEhYVh/Pjx6Nu3r+fdWvyCx52p0uw4YfUJROETBf428sxu4z9+rwugzMxMjB07FtOmTUPdunXx7rvvon379m43iwhwSPgQnmDXrl3o3bs3fv31V/Tr1w/jx49HpUqV3G6WJvxiRxIGXs9fA5DwEXX8LO117n7wsgjatWsXnnnmGaxfvx49e/bEtGnT/DofF+FtSPgQrpKTk4NJkybh1VdfRd26dTF37ly0bNnS7WZpYl7wSASq8DEzSGVxbm89Lw4bdgk/tWvo/P3gVQHk8/mwYMECDB06FCVKlMD777+P//u//3O7WUQAQsKHcI1du3ahV69e2LZtG1588UWMHTsWJUuWdLtZqvAJHqNB0OviR23gN4pfKWuinmATPkrMHo+Y+8CrAujEiRPo378/vv76a7L+ELZAwodwHKWVZ968eWjRooXbzVJFrOCR8LrwMTOIB5LwiQL/NfI/4SPhRQHk8/nw8ccfY8iQIWT9IYRDwodwFH+w8tgjduSQ8MnHi8JH79i1rpuZ8wV4RfhIeFEAkfWHsAMSPoQj+IOVh0/whMP89GXA+2tV8Q7mwSB81LhmYh/Aa6JHjtcEEFl/CNGQ8CFsZ+/evUhKSvKslYdf8MgxI37C4eWBLx8SPuzIj521jV6e1XYdL4kgpfVn+vTpuOGGG9xuFuGHkPAhbOWjjz7CwIEDUbNmTXz00UeesvJYEzwSPMJHXkYwCx/pnGVy1nEJ9h+/aDef1nX2+vUvilcEkGT9GTx4MCpXrowlS5agWbNmbjeL8DNC3W4AEZhkZ2fjiSeeQK9evdC9e3f8+uuvnhE9ISGvcoiecOjP8jFyWYVDvQyzcSFeQzq2CI6PiPp4Pm4SpfhYwR3RA3hngd2QkBA8/vjj2LZtGypVqoS4uDhMnToV9P5O8EDCx0ays7MxbNgwtG/fHlWqVEFoaChefbV4BzJt2jS0bNkSVapUQcmSJVGrVi089thj+P3334ttGxoaqvqZNGlSsW1nz56NmjVronr16hg/frwtx6jG9u3b0bRpUyxfvhwLFizAvHnzULp0acfq14OtAxcxaIoedCXR4LRgYhUVvNYbJ7FTGPG6+KTszlYFoPPwvTDYS506dfDDDz9g8ODBeO6555CYmIizZ8+63axC5syZg9DQUJQtW/T+6NWrl2r/feuttxYrY8+ePYiLi8MNN9yAxMREZGRkONX8gMft16GA5syZM/jggw/QpEkTdOnSBXPmzEFISEix7c6dO4dOnTqhcePGqFChAg4cOICJEyeiRYsW+PXXX1G/fv0i2z/yyCN4/vnni3ynXORv8+bNePnllzFt2jRERUXhueeewy233IIuXbqIP9ACfD4f3n//fQwdOhQNGjTAr7/+igYNGthWHw/sgscKvPtHQd/l4eTgqGx7MHUN0rEqRaXdFhb59TWyHLpn7VEiPUtuu78iIyMxZcoUJCQk4IknnkCTJk2wePFi3Hvvva626/jx43jhhRdw4403IjOz+AtBVFQUUlNTi32npFu3bnjggQcwbtw4zJkzB/3798eXX35pW7uDiWDq3RwnNjYW58+fBwCcPXsWc+bMUd1u7NixRf5u1aoVWrZsiVtvvRWLFi0qZiWKjo42XLRzw4YNePzxx/HYY48BAI4ePYrU1FTbhM9ff/2Ffv36YdmyZRg4cCCmTJniiQBmZ95Qzc7sUeKU0GF57C8hcNxxZlE7T3YJEB4R5A28IoA6deqE7du3o0ePHoiPj8err76KESNGuLbG39NPP42EhASUL18ey5YtK/Z7WFiYYf999uxZpKen46233gIA3HXXXYiJibGlvcEIubocgtcHXblyZQBAeHjxzpelrLp16+Kbb77B/v37cezYMXz66ado2LAhVxtY2bx5M+644w6sWrUKy5Ytw/Tp010XPebM8m68VYuKfdHCazEv/o7yXPJeN5aZbHbfEyyw3y9ecIHFxMRg7dq1eOmllzBq1Ch06NAB6enpjrdj4cKF2LBhA6ZPn67ZT7P03xUqVEBISAgWL16MzMxMzJo1y7b+Oxgh4eMhcnNzceXKFezduxd9+/ZFlSpV0Lt372LbLVq0CKVKlULJkiXRtGlTzJ8/v9g2Dz30EJo1a4YGDRrgpptuQqVKldC/f3+h7fX5fHj77bcRFxeHqlWrYtu2bejatavQOnhxrxPmnaUjBbuamQauZ4khkeMc5Qr+FRm8rSQC7ru52AWQm4SHh+O1117DqlWrsHPnTtx+++1Yu3atY/VnZGRg6NChmDhxIm688UbN7S5duoTq1asjPDwcNWvWxODBgws9AxKhoaGYM2cOnnrqKZQvXx7vvvsuZs2aZfchBA3UM3qI0qVL4+rVqwDy3WSrV68uZt7s0aMHOnXqhJo1ayIjIwMffvghnnzySRw8eBD//e9/i2y7YMECTJo0CTk5OcLNpJcvX0b//v3x8ccf47nnnsOECRMQGRkptA5exHS8ObDvsbDDdeT1R1jvmL0cEG0F/3NbFUXrnpK+1xZiXnB/tW3bFtu3b0fPnj3Rvn17TJ06Fc8884xqfKVIBg0ahFtvvRVPP/205jZNmjTBHXfcgUaNGgEA1q1bh6lTp2LNmjX45ZdfikwCeeihh5CRkYGjR4+iTp06qtZ/whyUx8chzpw5g6pVq2Ls2LEYPXq06ja//fYbrl69ij///BNTp07FkSNHsGbNmsKHRIvExESsXLkSx48fL3SR2cmJEyfQpUsX7NixAx9++CF69Ohhe516iH/TtNLBqA30RoKHN5GflQSIPPAKtWjO7Xlnqdh9zLzHW854k2KcM7GPU2uUSbDe/8aWKDcFUE5ODoYNG4apU6eib9++mD59um0vZ8uWLUNSUhK2b99eOKGjV69eWL58ObKy9J/vzz77DN26dcPUqVMxZMgQW9pHFIVcXR6iSZMmaN68OXr06IHU1FSEhITg5ZdfNtyvZ8+euHbtGrZu3Wp7G3/55Rc0a9YMx44dw/fff++q6PFCbIE2PLlbWNxdSvdVsAce+yOZ4HdDelX0sG3r5vMZHh6Ot956C3PnzsWCBQvQtm1bnDp1Sng92dnZeOaZZ/Cf//wH0dHR+Ouvv/DXX38VWu8vXLiAixcvau7fpUsXlC5dGj///LPwthHqkPDxKGXKlEGDBg2wf/9+t5tSyMKFC9GqVSvUrFkTW7ZscTVjqr0dqpWYihyIEyUUq+McbgjJQLi+xu13+wWld+/eSE1Nxf79+9G0aVP89ttvQss/c+YMTp06hTfffBMVK1Ys/CxZsgQXL15EhQoV8O9//1tzf5/Ph7y8PKFtIvTx5ycuoDlz5gx27tyJVq1aGW67YMECREZG4q677rKlLbm5uRg5ciQmTZqEJ554ArNmzXJt1pZzHShvrI+ooFYvPZI0pV0bM24uPeTX3a1gZqsJOyOgZ6EKCXnVNdfXPffcg19++QVdunRBXFwc5s+fj0ceeURI2dWrVy+00Ev4fD5MnDgR69evx7fffqsbgrBs2TJcunQJd999t5D2EMZ4qZcNSFasWIGLFy8W+nl3795dmNuhY8eOuHr1Ku6//3707NkTN998M6KiorBv3z688847uHbtGsaMud5RvPnmm/jjjz/Qtm1bVKtWDadOncKHH36IVatW4dVXX0WlSpWEt//ChQvo0aMHvv32W0yZMgXPPvus7UGCWnjTrSVC8EQAqAj+YF+jBIhOkwH+OB9/RS4IWa4Bz7WVumUvXVsjpOdA/550M/i5Zs2a+P7779GnTx90794do0aNwtixYxEaas3xUaJECbRu3brY9/PmzUNYWBjuu+8+AMDhw4eRlJSEHj16oE6dOvD5fFi/fj3eeecdNGrUCH379rXUDoIdEj42M3DgQBw+fBhA/jozS5cuxdKlSxESEoJDhw6hWrVqaNKkCWbPno2jR4/i8uXLqFatGhISEvD5558Xyd3QsGFDJCcnIzk5GefPn0dUVBTuuOMOLFmyBN27dxfe9v379yMxMREnT57EN998gwceeEB4HSy4J3i0rD5GYicLxnE7/rdkAaGF0iomUrA4ZQmyMhQo72XpfHhPAJUqVQqffPIJGjdujJdffhk7d+7Exx9/jDJlygivKyQkpMhL4g033ICqVati8uTJyMjIQG5uLmJjYzFkyBCMHDlSNXszYQ80q4tQZdOmTejUqROqVq2K5OTkYstmOIX7Vh75gMAjVrSEj14ZZqZ3e2mmUwSAGhzbH+dsi53T3+2cwWZ2jSWta2uHABIpfJTo36Nuub9SUlLQo0cP1K1bFytWrED16tVdaQfhPCR8iGKkpKSge/fuaN68Ob788kuUL1/e8Ta4L3jkmH0Tk8QPj2DiHdydFD4sxxEMwseMO8+uqfsiRJCdokfCm+Jn586d6NChA0qUKIGVK1eiXr16rrSDcBaa1UUUYe7cuejSpQsefPBBrFy5kkSPJZxYdsBu87gXlk+wm3Kw9xjtXFVbxKwws8fNs59+age3Zn7ddttt2LRpEyIjIxEXF4ctW7Y43gbCeUj4EADyZyGMHz8effr0Qd++ffHpp5+6MnPLm6KH16oiDUaBkJk4EI6BB+XyE24IPjNWPLMCSBIjZo7TTJ36Qt2N579WrVr44YcfUKdOHcTHx2PVqlWOt4FwFhI+BPLy8jBkyBC8/PLLGDt2LGbOnOn4ysbGb3xez3ciKidLsMyK4sVNAaYUQl4OQrVyD/IIPfk5MCN+tPdxw/pTuXJlrFmzBvfddx86duyIxYsXO1o/4SwkfIKcK1euoEePHpg+fTpmzpyJMWPGOD5d3biTC9f4v5NovYWLEDs8WZ619iecJQps181ON5ceLPelVrvNWLp4ngNp/TLjxIdOUrp0aXz55Zd47LHH0KNHD7zzzjuO1k84h5dfoQmbyczMxMMPP4wffvgBS5cuxcMPP+xo/XyCR/m9F1asZiET2snutAaeaLg3YDpFBPgTRHphwU+t4G3evD5aiA5WtzIdXhI/yvOuJ/SMFzItvr3+oqdOBj5HRERg/vz5qFatGoYOHYr09HSMHz/etdxlhD2Q8AlSMjIy8M9//hMHDhzAypUrVRNw2Yl50eMml8C2rpYeoq0z0uDkZrI7s1P+eTFTthtiSTkTzguCDSgqSnjTFAB8x6ElaNTK0BdLTuf8CQkJwRtvvIFq1arhueeeQ3p6OmbPno2IiEAO8A8uvDi6EDZz7NgxtGnTBllZWdiwYQMaN27saP1iRI/TVh+1oVj2AAAgAElEQVQR7ixRON0B+3M3YXSunHATytvgBREkuaXMWIDMLGIq1WN07N6y/jz77LOoWrUqevXqhb/++gv/+9//bFvdnXAWyuMTZBw5cgQJCQnIyclBamoq6tSp42j9xgHMvNgtftTaZGawtJr7hUXsiAgA1qqHx9IVw7jdMY4yszi25UHtuPTuKZF5ipRCwCmrnfz+5Xl+5OeKt62822u3y+mcP19//TUefvhhPPDAA1i6dClKlCjhaP2EeCi4OYhIS0tD69atC9eIcVL0sM3a8hoi22QmZsfuKdVemLrtRcJVPgCf6GFBft7dclVanQZvF/qzvpykY8eO+PLLL/Hdd9/h4YcfxuXLlx2tnxAPWXyChIMHDyIhIQERERFYu3YtbrrpJsfq1u+oomDd/C/a6sMyENhl9ZGXKzqzsRVhY4fFB2C3+jhp8dFCflws96yXMlNLGN23Ws+S3nkyEm1WRJ3+s+2k9Wf16tXo3LkzWrdujc8//5zW1vJjyOITBBw8eBDx8fEoUaIE1q1b5zHRIwJRlhk3cwWpTYvmtTAoZ4+RNcc+jM6tF0UPC2aTEtolArwz5b1du3b4+uuvsWHDBiQmJuLSJTcnFBBWIOET4KSlpSEhIQElS5bEunXrEBPD8zZuDXbR4/agbKazN9PpKd1dVvP3KAmHN4WOMueN/KPmbvNHF5w/tJfnPpM/E6xWMbXyrYoDY2uuk+KnTZs2+Oabb7Bx40Z06dKF3F5+Crm6ApgjR46gdevWCAsLw/r161GjhugYBXWMOyKtDthpl5cbM7WiTeynZz1QHoOTLiGtAT6Wo9w0xu3OcZTJA+ugzvvCIHfhGd2XTll7zNyv5WDuubwEJ0SPHCfdXmvXrkXHjh3Rpk0bfPbZZxTw7GeQxSdAOXbsGBISEhASEoLU1FQ/ED2AtTdl3qUERLi0eDp2yRIgYgAXtTwGD/5ohTGiIuw5DmXckhvXS4kVq6KZ6+18/IvTlp/k5GSsWbMGjzzyCK5evepY3YR1SPgEIOnp6UhISEBubi5SU1NRs2ZNR+o1dm3Z0RnydspuCQarsA6eVhMsamGXxcULOCnqvCCCWFGLGWMlC9bOpbkJC06u83X//ffjiy++wMqVK/Gvf/0LOTluZ5MnWCHhE2BcuHABHTp0wN9//421a9eiVq1ajtQrLojZTEZZ1v3tGGy0rD56nT5r4KtURixnmwjrOLEwaTjy7x+7rWeik2catVfpbuU9PusCwinx06FDByxfvhwpKSl4+umnQZEj/gEJnwDi8uXLeOihh3D48GGsXLnSsTw9YkSP/C3YaB+jzle5v91v2HLxI8JiECiupECDdWFSM3jtmmutLyfB217WbcVZTZwSP506dcKHH36IDz/8EKNGjXKkTsIa/mBvJRjIzc3Fv//9b/z4449YtWoVGjVq5Ei94kSP2r5q1hQzq0Y7gQixQ3iTWI3vlff3JfBlo1ZD9PIWZgOaWVGuRaYXXG9mzS9rOLXMxeOPP46MjAwMGzYM0dHRGDx4sO11EuYh4RMA+Hw+DBkyBJ999hmWL1+OVq1aOVKvtuiRplWzoHcLKsUPjzgoC2ez4eqtwq7GcfC5sGJgPKhK56ci3IvJCUd+O+txbB8oSNPzAbbB3WgGngihYGaRVK2XDr06APZ7TqtN9sTIOCV+XnjhBaSnp2PIkCGIjo5G9+7dba+TMEcg9TpBy+uvv47p06dj9uzZeOihhxypU1/0APkdm5FQYb39zFp5eDtwpxD52NGCpeaoaHP5Ii03ZsuShDiPgJKyi0uWIp7nh2fhU6X4sTcw2AnxExISgsmTJyMjIwNJSUmoVKkS2rZta2udhDkoxsfP+eCDDzBq1Ci89tpr6NevnyN1GoseCa2OlifmhmeAcnu2jFE+FmX7zLpF3JpSHkgzu9TW4rJKmsb3Iq+XlXLM7MfqKpPufZ5zKR2LM7OhnIj5CQ0Nxdy5c9GmTRt06dIFW7dutb1Ogh9KYOjHfPHFF+jatSsGDBiAd999FyEhIbbXyS565Mg7XJ6OUd7pGrkF9Mp10uqj5u7SaxtrYjzpXKRxtMUOoVKVcTtWV9d+xu1OMW7HCstx5MCeZIyAuGuj9XJh5HZV289oLTmt50hL8LMKGmetsk64vbKzs9GmTRscPnwYmzZtQt26dW2vk2CHhI+fsmHDBrRv3x6dOnXCkiVLEBYWZmt91ldWj2DcTkLtTVNN/LCW6Yb4YW2blvhROwdpHO0g4aONmeMwuofSOOoXfW2UQoY13ky+H8siukDx86Bn6TQSP+64op0QP6dPn0ZcXBzy8vKwceNGREeznl/CbsjV5Yfs2bMHiYmJuPvuu7Fw4UI/ED0Ae2wCz1Rhr8abWHWf2LnoI2EeO6ezW0XuAuOdlRUBdtEDFD12XveuHPfi75xwe1WpUgXfffcd/v77bzz44IO4ePGi7XUSbJDw8TPOnz+PxMRE3Hjjjfj8889tXyNGjOiRMOrojAYTeVZiMytI2400iJhZL4t1QI1lLC8c7FYNN+BdpNSLIld+zXjiteyMlTJzrmqA//mIAt8aY8o2uT/pwAnxExsbi2+++QZ79+5F7969KcGhRyDh40fk5OTg0UcfxdmzZ5GcnIwbbrjB1vrEih4JrQ6PteOtaKFuu8SP1YBT5artZnFzOQS1TMd6HzMoA5KNPlrYJQi9sByF9HJgph0818bMPSu1x33RI+GE+GnSpAkWLFiApUuXYvz48bbXRxhDwsePGDZsGNauXYulS5faHixnj+jRgrWzNWPKtxM9wcOSo8XJNbxEEI58a4W/LFiqJYhY28sap2RUtxynZ8YZ3Rdqixfb6cpzWxgWxwnx07VrV4wdOxavvPIKvvzyS9vrI/Qh4eMnzJs3D1OnTsXbb7/tcm4Ikaua83SwysHKrPgR0aFbtfCo7ZvGUcYtYB9AzFo3eCwo/ohI8WYUoO3UOdRaoNZs3VrPihULpfTse+t+ckL8jBo1Cl27dkVSUhJ27txpe32ENjSryw/YtGkTEhIS8Pjjj2P27Nm2T1s37gREdFoRYF9JXGtQ4okxUGLW3M47QErHyLJfrM5vykGIdTYUYDwjivd6sgbC3sKwzR7GskS5AwG+9rPcJyKvhRVYnydpppWatUcL6TyIED1qbXES7USLds/2ys7ORlxcHLKysrB582ZUrlzZ1voIdUj4eJyjR4+iWbNmqFevHtasWYPIyEhb62N/87EifuQiwKizNhIMZsVPBOe+Zi0CvBmCYxV/61mozE4Ftypcg0n4KFEbvN2aki/BKnjkxIA/o3SaiXok9ASkk+JHfu+7I37S0tLQrFkzNGrUCN999x0iIrzqKg5cyNXlYf7++2906dIFkZGRWL58uYdED2C+s1I+5FqxMHbFjphZVdpKO8zEdIieMh3ILiu70BJtXpzSbvYe5d3HrmfSqXtSWY96vXa7vWJjY7F8+XL88MMPePbZZ22ti1CHhI9H8fl86NOnD/bs2YPk5GRUrWrv1GRzDzuv+GENBObpXKPB/gbPGyfkxjpYVlf3lpDHrridOE1vdpfVGVpukwbvBHiztkOeMJN1n+OKfXhgcRfafZ21ynZH/Nx3332YPn06pk+fjvfff9/WuojikPDxKJMmTcKSJUswf/58NGnSxNa6rD3kPIsS6iGJH9ZOlTcwWqtcNfEjeiAzsvqY7fSVs47cmGWlJ1a8ZBWRsFME6p1/u9xcaq5UM/eAGUsoC7yxdKLFD8uz5Y746d+/PwYOHIhnnnkGGzZssLUuoigU4+NBNmzYgPj4eAwfPtz2vA/iHm7WtbqMymCZrcUz24Slbnmsj12iQW2A0jpnPNOo0xi3ExkjA7AFxjZm2GYHwzbHjTdhRt5uPdHOEp8EsMUoXYOzwkdZtxyWteGU+xidf61YISv5ekTF/fAIKedjfq5du4a2bdvi4MGD+O233yjY2SFI+HiMs2fPokmTJqhduzbWrl2L8HD7zL/i32jU2sojeiT0xI+RBUE+wPOIGCeSqkmDFMs11RM/8nPAGhwMiBU/gSB81JAGP5GB2UD+uRcdxMsTOH8N7AviyvdhPfdK8SPiebJ6vnj7Tu367BQ/x48fx+23346WLVsiJSXFkcWmgx1ydXkIn8+HXr164dKlS/jkk0/8TPQAxTsOM6JHDxGuLa36zcyM4cWONbxYLRMEG+EA7oRYN50kOEXGLPHOFpSyafPuY8e2rFg5T2b21d7HTrdXjRo18NFHH+Hrr7/G1KlTbauHuA4JHw/x9ttv46uvvsL8+fMRE8P7dsaOvb5rSfxYET1q08xZO+1Yxu2kup0MnOVxd0hTpL02gyjYsGsWl9V7j3f/2IJ/eY7hGGc90guHSOupWAFj777m6dixI55//nkMHz4cmzdvdqUNwQS5ujzCL7/8gri4OAwePBhTpkyxrR4nMpTmY9TBsnQw5RjK0arTyESvVb+ZBUZ5YJ2dFwE+a46e20U6ViddRoB/uroAtnYD+W03csfwuBd5XDvK+8ho31iV7/QEitrsQpb2WUkqqgfruRElXJyP97l69SpatWqF06dPY+vWrShfvrxtdQU7JHw8wIULF3DHHXegcuXK+OGHH2zN1+Oc8AG0RQtP58Q6C0etLrUBk6VuN8WP0lJmNshW6zhFiQgjAREBoDlDOSxvt2kM27DAmqnYrGBTGyzNxFUZDfJ694/avrEG5akJIK20Cnpts0v0sNQNXO8DeBMzqqFtubJT/Bw6dAh33HEH7r//fnz66acU72MT5OpyGZ/Ph/79++Ps2bNYsmRJAIkeQL3j4H0jMxo49FwQyoHOy/lgRExDdzLvzXGor3clPw5/NNmzWnvUEHX+rZSh3C+WYR/l86OXS8rNvEp69cqPwepzJO2v3q/Y2Y/Wrl0bc+bMwbJlyyi/j42Q8HGZ2bNn49NPP8WcOXNQp04d2+pxXvRIyMWP2Q5TS/ywuMEk8cNTN2/gKC/yWB+rgkcSfndaalFw4OSq8uGwPotOTWSwuEql/WI56pLuI9YEmsp22W3t0aoXUO8HrCwibFS2vf1pt27dMGDAAAwdOhTbt2+3rZ5ghlxdLrJjxw60aNECvXr1wsyZM22rxz3RI0fErCm524tnhhdPNmR5x2dmuQkWwsGXSE/N3aV2/CzxMoA4d1cswzYi3F1pDGWwEMuwTXMYuzqdPs8SvIJcSonAE2ycBj5XkV3PCAs5MO4HeI5FTyw56/a6fPkyWrRogStXrmDLli0oU6aMLfUEK2TxcYkrV66gR48eqFevHt566y3b6vGG6AmHmJke0hs0j+gB2PKXOJXxWHpjNWMNoBlezlBW8fEKZu9R1vslTVYPK24u18FybVjb5oUlR65TsmRJfPrppzh27Bief/55t5sTcJDwcYnXXnsN+/btw8KFCxEVZc9A5h3RIyFC/PCKHgk98aPV6Yl0eVmJjdgPtuNmiU/xVgfvH8hFkFvWHjksYkOZANNIMKeZqENuGXP6vuLpM43aJq0Zx1+fnX1sgwYNMHnyZMyePRurV6+2rZ5ghISPC2zZsgUTJ07EqFGj0LixlWBKr6PWmZgVP1JHzDKgaHV0auLHqFO0Kn70BI+R1UdUIKkba3ixUKrgoxckLardsQzbsLjlnIoRUqLmGjXTBl6LodtJDNWQt58nx5BRWTz1XsdO8fPUU0+hTZs26NOnDzIznYqjCnwoxsdhrly5grvuuguRkZH4+eefERFhT2fhvrXH3BuUOmrnSGtqstH5lOJ9eM67mTgG1g5ZbUDT2pdnmrXR8aUxlGVErOJvtTrjGMrZaPB7GktjZKjFdcQy7Gc1Hkler2iLD0tMmFQ/61pv0ktIGuP28uPTi4MSMZ1cC61+gzXPj9Q2K+U4G+9z6NAh3HbbbejZsyfN9BIEWXwcRnJxzZ8/P4hFDw9a50g5sLC+/cYwbicRDvakg9L2PMcvcjmDssgXGna/ecutb160JtlpOWKp1w3RI9XPKnoAvllcUvmAcfC3XedZ72WJZ9kbq+U4a/mpXbs2ubwEQ8LHQZxwcfmP6LkEY7eXUQcqDTA8HW0U2Kf5yo+FdRqxGawKRTuDcJ0WD3bjb8fDM/sPMBf8znP/sSb3FH1uWY7JvGixpxxxkMtLLCR8HOLKlSvo1asXGjdujBEjRjhcu5tJx4zQEj+snSbr7CjlgBCrsy3v+bJ6flktBPLAWquzjmJVvvM3USCKCADxuB5zpIUbCRl57i152gOWgXm/7P929Q8i7iEekaF3HNKzIkr8FMeuF8/Q0FDMmTMHZ8+exYsvvmhLHcEECR+H8IaLy24BZLZspfhhPT9SfUbmeq2OLlanTDXUrD5OC0oWscMSqyJht8gxit/xIqUUHx7ShLcmHzPPrp71Z7/Kd0Z1mM3ZY+XeMmNZUTsG5TNjVfxc0/ydXF7eh4SPA7jr4lJ7OO0YrK2WKYkfXtEjoSV+RL0tSkjiR7SINLL6SOLEqtUhGC06VpHPPnMStQB+rftOb003pQBSEz3KOpRYTVRoZ/4hNeTHoPWiYLZ8efC2s+KHXF5iIOFjM1euXEHv3r1ddHGZ+U1kPTywTnXXqk8pflg6tliYc205gRWBIrf6eEHoKJMCyj/hBh+rxAooQ7Jc6QnHNAH1sCI/N6wL2UbBWPTIy5cQlZ2ZNxbPKuEwto4a1aO8/9RmrDln9SWXlxhI+NjMlClTsGfPHsybN89hFxfLwyjigRX90Bu9xfDM3mDdjnUmjDTY8QacsiIP1ta7V1isPvL8OFrEsjXLECPh4qXsx6KxS1SyriTP+/zxtDUc4pekYKm7HON2TpUTjnzBwzdN306X16RJkzB79mxs2rTJljoCHcrjYyNHjhxBw4YNMWDAAEyZMsWWOvhcXHqw5sGwUgcP5UzWFwH2AV0pjvTehtU6UKuLUKoRy7idVhyPPCaFJb4mjbE+PW42+L2dwe9G8Qp/crRFDWX7lPe6iFxDwPVzaTWPDavoAYrmdjKymO5R/M3SzlPGm5hCq27lc2/mXKr1HUbl6J07Vku0eh9qR36f3NxctGjRArm5udiyZQvCwsKE1xHIkMXHRp577jmUL18eY8bYk9hKnOgxs4/d5l2l5YdV9ABsgznP0hdab412WH7SGLeTW320AnFZBvRgxA6LVJrs/1atQGb35XUPGdVhl+jRqltNsPCeB7UyWNA6dzyZ5p1zeYWFhWHGjBnYvn07Zs2a5Vi9gQIJH5tYtWoVli9fjjfffBPlypl9GM1g5eFj3depBzwT7HEerPEWejNdlC4vrwf/mpl1pCRWQDsCAb2FSa3MTLNTAKlNlNC6v5XWHqM67BQ98rolRPSRemWwXAPleROxtqB9Lq/mzZujT58+eOWVV3D69Glb6ghUSPjYwNWrVzF48GC0bt0ajz32mC112Jeo0Eurf0eALWGaVqeWpvib5bgk8cM6WDlt9REZU+Ll/E5uI3p1dtZrFsu4r9HsUPlzrCV6lHVIOCF65PUaiR6W88byHLKKH5bkqlo4O8trwoQJCAkJcWHijH9DwscGpk6dij///BPvvfceQkJChJcv1sUlR94xaIkENwZKPfFj1JmlFfzLE+zMk3IgHHwxGWZRG/xYLBFKd5fI2VL+gNX4I8D4nKVxtMdpFxjvmnhOih7g+qw+I/SOWxI9LMdqdP5YM1Pr4dyzVblyZbz++uuYO3cufvrpJ8fq9XcouFkwR48eRcOGDfHUU0/hrbfesqUO87O4jFDrFORvPk4OlmptUb5988xOYZm5pew49d6UleeCNfsyD7EwPkaWOJ6fGbYRHUAMXG+71eBmlinYWsH5dgde58B6gLgUeBvLuH1zsA/QUrZv1skLUsC+mckOvCifZ5Y6lUHKapYeFmuNWrBzlsHvvBQ/HrsCnZs3z5/ssHnzZgp0ZoAsPoJ5/vnnUa5cOYwdO9aW8tVFTxSsuz609pfEgNuiByjaMfEmOjQaPK0mOhRt9WF9u9ez+ti5hpeEZAlxe7kLran0drdD1AsHbztZrq18iRMWK598lqLdz7ta23ktP1ruLV7LTxaKC0m3suvzExYWhunTp2Pr1q2YPXu28PIDERI+Alm9ejWWLl2KyZMn2xLQbPzQmB1wjPYpC+fiflhM0WazO2uhdWxqieH0yhQhfuTXMM3E/mqxKSzuHCOsJBd0O72+naLMqqVMIhZs7VKmMeAVtlrXTC01g13iR6/NLHXGwjimh6W/UhM8cqzcJ05YzK7TsmVLPPnkk3j55Zdx5swZR+v2R0j4CEIKaG7VqhV69uzpYM1qDzjPA+t0NlU9WNvCklRNrQNVs/oYHZMkfpyIi7HS0Yq27ojOnuw13LRO6aHXHq3cTWrXfofahgUor6lePirR174iwzYsLxcsfRFvVmZR6IseuwKdJ06cCJ/Ph5deesmW8gMJEj6CeP/997Fv3z4XApq1iIC4tO9y7BI/vIOPlvgxGqgl8cMze+1O1kbB3PRwvcEuTWc/6VhZYnhEBPEGGvJzoiaE7E6sKBGr8p0V95ee6JETDrYknKLuB0n0sBwXizvZSl8k5QkzOjaea5ADVkuPHeKnSpUqGDduHD788EPs2MF6DwQnJHwEkJWVhddeew1PPPGELYuQaj8kVt96jB5qrU5B9JR3s2/cSvHD2kGncdTBOysGYBc/Zq0NosSJfKBnXe/Jn2AJaNfCK9YgqQ1a1h6tfURj9Z5TWnp4xY+WG9no+VT+ngn+5KgsbXXWtaVF//79cfPNN2PkyJFuN8XTkPARwNSpU5GZmWlbQLM6PGZcNaFiVvTwtIEFq520JH5YO2WpPqPcJspzpidoeQdI3u3TCv7VG3xY42jsGtBZF7+UUXg4Rmt9uY2dIiiWcTvWRJVSRm+WtpqZiWjmemi5t1jFj9XYuSioCx4etNrKbuVRYofVJyIiAuPGjcPXX3+NDRs2CC8/UKDp7BY5ffo06tSpg/79+9uyHpd5a49WByVN9RQhfJRl8iJyIKlqsj41S4feuVWakPWOIY2xDSwYTcsG9N1ZrO4bEwKmEMnConLvhMdr75azTqdMluR7WhhZfPTOF8t5EvGWH8uwjTJlwd8a22ktXqs2Ndtq+gWWY2eJ5wH0p47HMGwDGPdBLC49o2NStkGMlUf0FPe8vDw0a9YMJUqUwMaNG20JvfB3yOJjkfHjxyM0NNRh06KVoL0oGHdIvG91Ziw/ot+ejRKvadWnHFiNjkVu+TE6hljFtlaOmSWeRDlYO+2u8YqFRkLPaiR6ppsZbja5L+8yJXa5vkRsA2i3L0bjezW0ntsMsC8mzOryMm/lcYLQ0FBMnDgRP/74I1JSUtxujich4WOBw4cPY8aMGRg2bBgqVaokvHztnD16GD280v5aM4DMduJeWOpCS/wYdfyS+GFtfxxDmTz1i8YLsSleQ+E+E67PrAggvX21ElQq12nTsvZIyO8HUck29Y63KsM2eihFD8v9rHx+WQUPD+IFjx0ur/vvvx9t27bFyJEjkZubK7x8f4eEjwXGjBmDChUqYOjQocLLtmfKo7JjsCPBHUv+IjsHZaX4Ya2LVfTwnK8IWAuwVaJn9ZHEznqDMkRYOgIE06FEeu5A1sLUXJfKfVmycpcCsI1hOyD//hCdYVztWJVuZ57EhDHQtvTwiB+zokcvREDMoqVOMWHCBOzevRsLFy50uymeg4SPSXbt2oUFCxZg9OjRKF26tEO12pWXQhrMrb4GSx2T1VWSrSKJH55EhyzTP+Wix2iWjbzuWMZ2sKAUPyKtO+EInNldJo6jUHesE1C/FTeYmf1Y7oE0/qYwIW+rVqwdy/GwvCSwHKdVS4+yrfYLHjtedJs1a4Zu3bph9OjRuHz5svDy/RkSPiYZOXIkateujb59+wovW3stLr0AP1YXlxYVwWatYaWcSnlOul9YkxzKz5ue+FGz9KiJH6fcTCLq8drsKS8h8rwoy2ENVGexLm6U/d9N92Y4jCcY6J3L2IJ/rSw0erzgI+I8hMNpK48d4mfcuHE4fvw4Zs6cKbxsf4aEjwk2btyIlJQUjBs3DpGRkQ7XLmLxPCXyjsKs+NHqbMoZ/G4HUgdrJiOtUvwYuQPl4kfvGGN1fuMhHMZv7gYzksLbwR6hozMDK2dj0dU0iqysobakhI0CUm+GmfoOso+VWW9mhJTe/bdR5Tut85bGWS8v0Rr1KlE7/ljF32bEj9KFZ/Xe8S+3lhYNGjTAk08+iddffx2ZmVam8wcWJHxM8PLLL6NJkyZ49NFHhZfNtvL6NRQVQFatPUrUrDV6GHUyrNNaRaA8F2bWIJLED288jxGxHOUpsWh94DJgWHF3OShg7EB3ar0I6sH4nKjFYfHG48nLT+PYzwzydbN4xU+sxjY8fZZW3JKZ+05u5XHeCmqH1WfMmDG4ePEipk6dKrxsf4WEDyebNm3C+vXrMWbMGISGun36lAJIDaMORK9zEBmobPcq4YB2R2Vm1WnW/DGlAMQzbmsGLbXCYHnwFy9WBZYgXl78QXyptcko+Lws8i09atYetfLTONvEi9pioaziJ9ZgG6O+KwNig7XVrDxef3iMqVGjBvr374933nkH2dnZbjfHE7g9cvsdEyZMwK233orExEThZbNZe9SwyywbhfyOzWglZD3k7bdjFplaPWpkMGwjIXXcRlOE5dOJjQZvSYGwxHfI9zFDOAwDdLndPYFAgQAyfRtaSaYIqAfvmhFlrPdFmsnyWdHrF4zqrAdr6S+OFfzLmntHD2/N2LLD6vP8888jKysLs2fPVv09Ozsbw4YNQ/v27VGlShWEhobi1VfZ2rF8+XJ0794dtWvXRqlSpVC7dm0kJSXhzz9FrWUnHhI+HOzcuRNfffUVhg8fLtzaY/1mV3t4rVh7lJh5s9PqlESLH9aBgOXtUHlMWuJHLYmclvhRts9I/LCaaZRWH8V+trtt/BC5dUkZa+Tq+erYSfYAACAASURBVIoA8E+G7aT4LV5Tnmjxw7KMhFqd9VBUAPKutQVcFz0SVsQPi+Dxf6vPTTfdhKSkJEyZMgVXrlwp9vuZM2fwwQcf4Nq1a+jSpQsAMGd8njx5Mi5fvozRo0dj5cqVGDduHLZt24Y777wTv//+u9DjEAUJHw4mTpyIWrVq4bHHHnOoRjMPnCSArIoetf15fPlOdRa89eiJH61jUoofvcy5SvHD0z4npz9DxeojdwuJX2zXOcy2Xc8tZoe1R6sNPGhd9zQBZWvBs3aWvE6zOa2kvugYioseM+2Q8I6VR4kdVp/hw4fj5MmT+Pjjj4v9Fhsbi/PnzyM1NRUTJkzgKjclJQXJycno3bs3WrVqhZ49e2L16tW4cuWKZ+OKSPgwcvDgQSxZsgQvvPACIiLEvj3Zk6wwy4YygXzxEyugnIqyj1nMiis18cN6TVmXCzASI0qrj9ljsTDLqNiMKhGw5EPyB+TnxaroYUEe26N1PbRm6ynvtTSdeqxeZ7no4XEd64keFpeXUW4enufHjGvL/60+DRs2RJcuXTBp0iTdbM68y3dWqVKl2HfVq1dHjRo1cOyYWaFqLyR8GJk8eTIqVaqEJ5980qEarTxoUueWBXUBZMbao0TvrY/X9GxG/FjtiOTih2Uw2AZ20dOBox02Rx+ruW/sDLUqRtT1T3lof+TbFfnoTXV3KmjZqvXrFuhfZ62AZvnxGS2aCrClOpDKNYPaM89y77LMEtTrc9IK/rVqZY6ANSuPs+LHjhfil156Cfv378fy5cuFly3n4MGDOHLkCP7xj3/YWo9ZSPgwkJ6ejnnz5mHo0KEoVYp3gUB9zAc08yAXQFZFj/z3GijeGZptu5NT3uWwDALSNiwzaSRFYTQ7B8gfEMya/+UwWCSK5M2RoTujimfAlwsRuXBxAqUQMlH3eb1rK9KKZeb5sEPc8ZZp9kVHEj0s10NtmzSG/VjILPgEN02bNkW7du0wYcIEbssOKzk5OXjyySdRtmxZPPvss7bUYRUSPgxMnToVkZGRGDhwoNtNYUCvQ8uCPQ+/WZ+/koowzv4KiBOGx2HcsSrbqzVAqqkKljd5u4nITx4oukzPTxUHilmPatsxbd4INWuH3PrDIpBXg+0c886iYb1usQzbqD2TymPnET9pMBenpNYOkX1eYFh9fvvtN6xcuVJ42Xl5eejTpw82bdqEBQsWoEYNnrHBOUj4GPDXX39h5syZGDhwIMqXLy+0bGesPWpla3UEZmZYSNSAcQfJ2mnpiR87zk+axvda7VUKCauJDu2w+ogUJIJEzqGtgtpjgWIuNjcJZ7id5S4uvfNvduqw0TWNlW1nhPxgtNxbLOInzeB31n7ELiuPf8f7JCQkoHnz5txBzEb4fD7069cPixYtwvz589G5c2eh5YuEhI8B06dPx9WrV21ZgV2dHAv78gxKyk7BiuiR121lZWU5auLHzg4nTfG3UXsl8WMkepRBq1qIED9SHSr16LlyVN1d8hgbf+JOvs3LAzi/Fdr3txU3l0FsizSrjvu2Vl4TEflS1K5zLMM2SsJhHNOj15fsh5jn3LsztngRbfUJCQnBSy+9hO+//x4bN4qxBvt8PvTt2xfz58/Hhx9+iB49eggp1y5I+Ohw9epVvPfee+jduzeqVasmtGz9mzkH/ALIbOCfHW9EMdAWQFqota+qwe+iSSv4l3Ww/5lxu39ylGkFq7OP1GJzOIUEAMADlh1TiIxNMoiPUqYSUI191gtolgSuyCRx8ns0lmEbJY0LPmbiefaj6AxFs7l5tCZ0iMau/sgZa1JiYiJuueUWvP3225bLkiw98+fPx+zZs/HEE08IaKG9kPDR4YsvvkB6ejqeeeYZl1pgRgCZIRP600XN5gSKMfhdQu9hrwq+GCKr6OX5kSO12Wi2jXTsrMHODlNkVpULVCr4FFspXvBst9q8Ak6ExcuEa7DwsFlmcUk7iCQC5lzWSqHHcz9ppWTgFT9KweNP1kr5/V78uEVbfUJDQzFo0CB8/vnnOHHiBABgxYoVWLZsGVJSUgAAu3fvxrJly7Bs2TJcupRvQevTpw8iIiJw9OjRwrL+85//YO7cuejduzcaNWqEn376qfCzbds2oe0WRYjPrtDuACA+Ph4+nw/r168XWq75m1irI7A6zVOJPFGhqESIaRq/s7RNqsMol4dVWN54geJt1gtilqM3mMnLtGq50bE2VIjTjm3RjMExY8Fprv51U422bdFyKRktGyJHR9zoCR/m4zZaF09C7Rhl+7IsF5KzDsYvPXLBIOoFKbbgX5ZjlbbRs24ZuZyMXIksx3VO5zfWa2YWq+ddq/8rXq7PN8ZiXUXJzMzEjTfeiBdffBFjxoxB7dq1cfjwYQD57jBJGoSEhODQoUO46aab0Lt3byxYsKDwbwCoXbs2jhw5ojpLLDY2FgcPHhTabhGQxUeD3bt3Y/369R6byWXGAmTmjTAD4kVGrMn95OLBypphPPUAfEJNTdCoCUI1geTkKqJRJgN6zbi73GAruC1GXAHXVgK8C/ZjFj2AfvvVliuxSqzs/6xpHqzkONrDUI/RcemJHjCU7xZG96f9Vp9y5cohKSkJs2fPxrVr13Do0CHk5eUhLy8Pubm5Rf4viZx58+YV+RsADh06VGR7+ceLogcg4aPJzJkzER0dXbhuiSjE3LxyAWTng221U1Fai2JRtHM148e3umgqaz1AcfGj115J/BgNjHLxo1WeVZeX8i1aFreiN9Bzu4P8gYIBptClxoqRINK6zgZCwCh5ZLGEk2oDpFnXkB6xKt8ZPd/NYRzcr2Yx3oOiVk0z4uccjPsn1vKtwHvOnXzRMWbAgAE4ceIEkpOT3W6Ko5DwUSErKwsLFixAv379EBkZ6XZzdMiBvilZxAN2HOaWeNBzkcVCzKrKIjCqJw3snRWPS9SoPBHxPqICdb3TUeuj4V6TU0nxsYxJCxB35mzWa2BmYI3l3B5gOteFyO9Bs25c+TGxCh4vYea62G/1uf322xEXF4cZM2YILdfrkPBRYdGiRbh48SL69+8vtFx71uQCzK09Y4SyM9cSQGapB2tTuNWyRvPCOmCxzJyRyjIISg0HEN6OzeVhmgjoDjBM7h2lu4hnoAP4YnNsRjeuSFQgdQSAFvqbKNMGKK0/hqvDh4N9bTbW44k1+F3tGVHeCywqLgr6ooflWeSx8pgp3yx65zpK5fcrAL4FMA/AeACjAKxlru3UqVPo1asXqlSpgtKlS+Oee+7B2rXs+ysZOHAg1q5diz17nFiTzhuQ8FHg8/kwY8YMJCYmombNmm43hxO5ABKx1pcaxw1+B/hyAqmJHx7Tt1nxw9sR6okfxiBmrktiJn7CYoLBSgCa3gn/sfCIRhJAVqbja1jZ9JYGKQsG0QPwW0uMruPNDNsA1++p5tAWwHriZ0fBx4qV95TBviw4KX7k94Gy3r8BbAGQB+BWznJzEB19G1JTUzFt2jQkJycjOjoaHTp0wPfff2+q5V27dkWVKlUwa9YsU/v7IyR8FGzcuBE7d+70WFAzL3ZYgOQcK/ioYSYRotz6Y6ZzEmH9YUFN/Gi1VyF+1Pp8TauPvEwW8aMleBgGcKFuHyMidBYpVWaHFpl52gmU8VG8bkaj4xXhIpJzM8M2cuIZtlETP1bXOTsFMaLHKViuewUArwDoA6A9Z/lbAZzCp59+isceewxt27bFsmXLUL9+fQwbNsxEe4ESJUoUJh+8ePGiqTL8DRI+CmbMmIF69eqhbdu2Qsu1z82lhtSRmRFAPNYWPQFkF3qdNKv4sTKgysWPUTkF4kd38ka84guetlkQB4d2OCB2VETMajvN6V4USlEGC8FCkVVb7RisnjPlDXgzwzZypPazLNAsFz9qoofH6mOH4LHrHomCvtVLq16WbDLyc7YXQGW0bPlt4TdhYWFISkrC5s2bcfLkSYbyivPUU08hKysLn3zyian9/Q0SPjIyMjKwbNkyDBgwAKGhgXRqRFmAtDotSfxYXfaiMfSDelk6zViGbazyJ0c56xi30xMxalYf1vqVVh+GeBatmBjuOB83UIitdnbmDZIwmA2nty6Y6lIidli7pGuuJnqU28hRijYW8SO5tnjqkRMB/7HyiJpAwMopqM1qve222wDkp2ExQ61atdCpUydMnz7dtlXbvUQgje6W+eSTTxASEiI85bY71h41jASQlc42A/qLC/KIoltQXADxxJ3Ewtz0XFZYg0wZ21whnmEjaQA3OyiqCB7NwT/AcHNhUnmKAO66DYLUuWGZKSi/R7QsVXriZzOs9yMZ4HveeRHVD/AKHiv1yq34xeutWLEiAODs2bOma3jqqaewfft2bN++3XQZ/gIJHxkLFy5Ex44dC2+iwMWMBYi1I0qD8erKrPBM6VbrVGINfjeD/Dzo5VORbacXuCpZx43cIQDMH4PVtbP8IeaG5fwVUCSmyCo6VjC1vEiS+Dq/UX/hWADW42Pk8MTQhcP4fKqJH7m1zExuHrszs4tCz8pjl2Arfj5Fv1C3b98elStXxsKFC4WW60VI+BSwZ88ebN26FUlJSULL9Y61Rw25ALIyCKjtm4brAsiKC+wWGGcO1mt7rMHvPKidX8YMukrxo5bITlP8SB2tDQkGtaw+7RpDXRxwCAwYDewusVqZPE/6mJ2Cb2Y6vNEzYYfokTB6HtqBP9GQ2rljFT9ameK9ZvUR4dbirVf5DEZB7aX13Ln8af6VKpkP3IuIiMC//vUvLF68GLm5uabL8QdI+BSwaNEilC9fHg8++KDbTXEJvVXarXRARsHPrKJIK+6EpSOJZdjGCJblAwzOkyR+9MaUYuJHeX7Mih+GQd0td5AmFoVTO7sXfVVae2QCiCkLthPxIVq5stSem3YomlncSPyUQv59ZSVnk5tWHlYRwnudRAg2tZeOKOTH91w/Z9KL9c6dOwEAjRo1slRrUlISTpw4gXXr1lkqx+uQ8EF+7p5FixahW7duKFmypNvNMYmVh016g8iEvgBSg6Xz2A/2xGt6NIa53DaAdtwPCyznlvH4cngGc4eCJrfssVnsKNbPqo7in2JrbLF+vIhBu4olkFReZxHWHpYEofJnt53mVtpshPHzr/U7a0JUN69xOdjzDGqdE0lg6Z3ThgDOQP5CmZOTg4ULF6Jly5aoVq2apZY1b94cdevWDXh3FwkfAJs2bUJaWpqfu7lEIhdAVjoe5b5KAWTWBSYP9NXDKO6HBZ7jNwpELWiPUWxH7TjonxtRVh+DOBdui4kb4kRUfbzWJZ3YnqaNtXMjaWbNjoLxbCgeWM9DBPRFj5bVx2gKvrIOOSyCR47TLq9yBR8rGLV5H4BdyJ+eDgBnAWwv+Egryi8F8BKAv2T73QOgasFvOwEcQPfu3bF//3688cYbFtucvxJ7UlISli9fjkuX7MwF5y4kfJAf1FyzZk20atXK7aaYxM6OIUvnN7OxM/thfep7HPSnVxvF/bBg5rxqiR9Fe7TEj2R5sXXBUKuBvZI7zqvWF4VlSY6teYRQPBUAV8iFqIBrSbCyXBMWS49S/JhxQUZA/LI3olETPHYF9KcA+B+ALwr+3gHg44LP5YLvfLKPRDiAx5Hfh30DYDE+/3wLVqxYIWz86tmzJ7KyspCSkiKkPC8S4guGSfs6XL16FdWrV0e/fv0wceJEoWU7Z/FxaiaBsgPU6xRYE5VpWRWMhI+8LVYCK/WWorByXuXHpdMWeUyPmrtJ00IQDvOxFRoByloWniJiQTonrIOfSl3/VnHBfKzmKrRQh249OSbrkdd3Tf0nzRxIALboza5T/qZRviFq11B5vBJy0cNS37cGv+uVkcZQvhFax2GVsgZlm70WUClX3m/pWVT02lN8P59vDEebjGnRogWio6MDdtX2oLf4fPvttzh37hy5uZjIwnULkKhZYHtQ3ErCI3qA4msI8Uyl1UroZlVMSsfEeJ60YmyKWX7k1hWnEgpaseqoiApVkSMYNdFTiAgrlYp1Rlf0sKxVpSyfF62XCLV6lZYeo/pWa5RjVEYaxKW3EI3a1Eo1RFh9WOuS0DvXxftH0eNNUlISVqxYgTNnzggt1ysEvfBZuHAhbr/9dsvR8O7hXN6IouitkmzFRcQreuQ0h7n8IUrxI+qcMrhVzm/lCCxWa5cZ8aNh4ZBbdooEHqvBM63dI6gKLl4RpDzuAgFkKHqU9cnRsgTxDLhG8VjyOrXcW1r1aSy6a0iayf20EPVcqokQu/rRiip1STiZ8ZmPRx99FD6fD0uXLnW7KbYQ1MLnwoULSE5ORs+ePYWWG5jWHjmS1ecc9AWQGnqd+R5YD+6Mg7lBWRI/ojtAveMp6Pg03VkF1HZwxXQ1saNrPbFITMGncEaL/OMEchFowQrEPStOqscouSSL+JFyLrHUaRTTIy9nNYqLHharTxrss/KYuT5XACQDmID8xUHHI391dBb+APAegDEARgJ4E8AG5K+sroeIAGk97H0+qlativbt2wfs7K6gFj7ffPMNrly5gkcffdTtppjELWuPErkAEjELTGt2C0teEQk18WPUtltgPAXYDGrHoui4tMSPNDtIz6IgxOpTMOjb4oZSETWFgqeAl2LUdnQJLQFkEBullgNJd1kQ1ozaes+j/L7Q2S48vuDDUp+a4ClSmM5vafrtcIUlALYB6AxgEIBaAOYB+EWxnfK4fgcwo+D/jwHojfwXpC+RL6TklFN8WNETMHx9qegX7sceewybNm3CiRMnhJbrBYJa+CQnJ+POO+/ETTfd5HZTAgQj6w9vhygfNHhEj4Tc+sPTidgtfjQ6O6X4Uc4I0hU/ZqxcXpyRZQaFqHqJNXjaCJPnRhI/QtdCU5vxxbh4bXi84m+WunhJg3OxPKzXpCLyZ5HtAZAE4D4A/wDQE/n5cD6HvuXml4K6+gO4Dfn9wkMA6qOoaIqGvoDxmhBko2PHjggLC8NXX33ldlOEE7TC59q1a1ixYgUSExOFlquuuv1pcLH6kJ6CuZWVtc6PUwtpKo/bLvFjYKKWxI/WNGhV8SO1nUf8hINrNpOqu8tM0jvRaLQhRvFRhfX4wwG01qheI7amcD0wLczOyJPK1BPBsnqVoscQydLDEyuXpt8Gx6mKfNED5Ft6SgK4S7HN3QAuQF+sRQAIw/Vjkf4tWfD/aKitlC4WviBnkVSsWBH33ntvQM7sClrhs2HDBly4cEG48NFH5Bu214WUUgBZ6Qj3QD81vt5q0QDQAcXT8cvRahtL9lseIsDk3jjEI/Z4zyvD/adrHTETg2M2OFYkImKHOJ7d1Xqz+qws8QAALRi2idAXPcUOQ829ZXRvpcFYONiF2nWoWvCRcxxANRQd6iIA3Fjw/5M65bZCfg6d5cgXSX8jP/HgbgD/x9levXNh5p68BrVp9qLdXYmJiVizZg0uXrwotFy3CVrhk5ycjJiYGDRp0sTmmrQ6Si+6GezoqE7B2AXGcw6UAshI9ChdZLzTeAG+VeK1kNejJ34KzoWem6Sp1gKigLbVR+teY7B6FAlADgR4BZDynBacS6a8R0DRa2VV9DAmttRc8FZGOGAcz6OFXv4rp1ETPNK9ng2gtMo+0nd6A3pNAE8j32o0GvlZlOcA6I78lykldj8f16AleOyic+fOuHz5Mlav9sLLiziCUvj4fD4kJyejc+fOCAkJEVauebXNK4K8JJZYMbsCs5YwsTKA6Fl/lEjtsyJ+1I5BTfwozoWW+CkPg+Uk5IOeBXHNFIDsBXeXAROOqXwpoCPXnOqvhUjRozPISqKHKW2MmeePR/TYafWpgeKChxW945Z+OwDgA+THA70E4FUAXZC/XMQX6ruaRk806Ykde8eCevXqoWHDhgHn7gpK4bN7924cOnTIYTcXK25ZgezqoJTHoiWAzLAN+hYLo57/nwa/K9tuRvzonVe5+GG85vJZQ4ZrabHeR/JzWGAN+VhNLIgiongMTkzB90U+JhA+O8wgK7Ry+r9tS2KoWXpUBkulpUfrEchZl/8BwC5+/oQ5S4/ovoU1tiYcQBmoW3Wk79SsQTHIF1UpyBdWw5EfI/QP5M/u6gTgM6jHMooKcr4k+/Bhh7srJSUFubm5Qst1k6AUPsnJyShTpgwSEhJsrsmsgImA8/lMnIZFALF2FBtRXAAZiR6pbC2Lhda1uwXsAoi1/Tr3iWT1UZsqrVueRjCuJir3mpqlxLSwsChoDMvlwaK1Ry3QuzoMRI+Z9a0k9NxbBdesQpy2e0v5KBQKHh684NpSEzxG1z4G+XE8ytlb0pphNyq2le7vCABHAdQBoPQKNER+7I8yPsgq5sWOnSQmJuL06dPYvNmqxdI7BKXwSUlJwQMPPIASJUoIK9PepIVyEeSPbi49tBYuNOrQ1H43O7jwuL4k9GbV8LAVhi6QLToDajGrj/z+4Jnl5bAPf5Aoq59EgQBSuudU3Vys8N4TFuKoNGHI0VSb4TpL4kdT9Gj1K/sLPlb7HSuCtwbMz55qivwEhr8qvt+E/DeJe6E99a8S8gWfUjT9UfBvRfCjPA9Zso/Zc2TvmNCyZUtUrlw5oBYtDTrhk56ejp9//tmjbi4W7LAEOeXm0kNLAKmh196NAH42ub800LGa/hmnFDOhJX4KytGzJhSKH6dFsWRtMXIZ2oVOvcXcZ2ZRPGtaWayLzISTXwcRokfnukprubFYA3OM2qKsRzm7z+n7q0bBxwi963sb8i20i5CfcXkvgAXIT074NK5bc6YhP37ntGzfLsi3+kxAft6e7chfPT0Z+c9+TY06jfpmudixB5Ev4mFhYejUqVNAxfkEnfD5+uuvERISggcffNDmmqy4uVjLDER3GI8A0kNrtorBIBjeziD3iVoSOaUAMjvQKsWPohwt8VMdwL+13G8WrT5Ki0kMgOkxxdtmyArO7U0wXc8qIAkgVsuW0tqj85xprgEm0tKj0p8oF7DVEj/nN+Z/mNgv+4iG9Z5REzxWRNd/ALREfsblaQBOAHgZQBtcv6a+go+cfyI/vucygJkAJiHfQtu94HteMgs+/kdiYiJ2796NAwcOuN0UIQSd8FmxYkWh6U4U3liby6wA8oK1R400mEuMpqzXpAuHO/GbJH6snk9J/GiUoxQ/8sBazTW1BCwoqpsIMEh4KYbzHJidqKDl3pKVpRQ9EkrxU0zwWL0/7bT6xILNwqOG3nGVADACwDIA3wCYheIxcEOQP1OriuL7ZgBeAzAX+VajqQC6FpSph9QXZ6K44NE7h950d91///2IiIjAt99+a2s9ThFUwicvLw/r169H27Zt3W6KjQSaFSgN7KnwtR5+yfpjZO1R/h2vEEBGnRJLYjkWDGJ+JPHDPZ2aBblQLLCSTGCJx3HL3aWD6TgildgeeVC3XARqJnuUCw6eQckopidcW/Qo0bTyqN3Hewo+LFgZZNXqji34iK63HtiSkJqdiaW1nzRxw13rjsgX8jJlyqBFixZYt26dsDLdJKiEz+7du3HmzBnEx8fbXJMTbi4WomD/KsFq2PH2kQZrLjDJzaFhAdJN6xEP9jcxxkHJEAPXxGqNAbeY1UcSwTyBugxxMbpuJaNyCz4NUfTDNZ3djNAS7G7TnPKvdu1Ynok4GB5308bay5lInN/K4doC2AWPaGJRXPBYsUhJ+6qJHbtfBDMgNlWH94iPj8e6deuQl2e0Mr33CSrhs27dOkRGRuLuu+8WVqY33FwsqK0cbGdyMbtIg7oFiEdsmXB/6U0XBlC0Y70TYgSQ1uBlcKyF4kfZ2bOKHxWBIHwWllVWQEjen2KoCCrdKfxK66qe4NC7bvJ7S+NYdBepLaBwoVujgT4CfFYeJVZebuqBzcLDW+8tsGeNPT3OwZzY8T93V0JCAs6cOYPff//d1nqcIKiET2pqKlq2bImoKLvVv105S0RhpxXIrodPeS7TwO4CU7sOMuuPUZPleVBUxY/W/WSH+JE1Vm9NLbUVym0lAoCTMyUVda2pWNR6NMhomRQ9GJ7dYtPko8AWzKx2s6ndU7L6mzYuLnrUrD6HlNnA9fo5uxb/vYL8ta1GAXgewBu4nqgzFsaCZzOA/wLoDeAJAC8AWGOwD2tuLb3zwePuOg5xkzDsReSLecuWLREZGYnU1FRhZbpF0AgfKb7HfjeXElYRJNrNxVJmWdnHX0kDuwBSY51+Qje1U1PE+uNELNVGaAbKKsWPYRCyQKvPmmjwCXynp8PKnz0WN5fS2hNRMINNgaXcQABzriU1wSNHEj+HtqqIHi12QJzoUeuX5iB/6vc/kT9dvBaAjwAcVmynds98BmBywT5DAQwD0B7A/7P35mFSVPf6+Ns4MExgQAeVVdMubEFwAdFAEvGCiSTKEzXxxgXvNTEmoN98Y676EzERTKJGkzwxCclVkYgbcc19MBf9Kl7QBAMjEEW5CKh03BDjjMigA8ww/fuj+nSfOnX2OrV0T70+ZQ/dVadObee89X42NmNwHSpkhyU8UT2PaSc7naUlOnzqU5+qGT+fbkN8XnnlFbS2tsaQrVn1VpFW5+M0kyCdiVWUTl9zUval8HeJidRiC0VoNCE/7BwdunyDgCyUfXJqHMSPSes0mppP66CMttPJzbNNRWLocSYqlYdgI7zq5ecCmATP7HQNgOMA3I1gIkAarwNYAuBCAJfASzw4Bh7x+SK1Hi99hCuwY/Lb1BK3aUrVZiezxINTTz0Vzz77bNX7+XQb4rNixQrU19fj5JNPdtZmeBkxrSSoiVpMkIas0ib1hDj9pcmPigMeMdmLsBFG2bDthyE/gJz8CFQILvmxcQymHJJdYWoYc5RrKM4JUdKcFTydDOnzQpJSyshPuYit6rnbhOhID73vDfDCvI8DcHRpAYCp8HxhtjDb0pP7EwB6ARDlVzMhPLYmLcBPdtKG+IkOiylTpqClpQWvvPJKYn1wgW5DfFauXImTTz4ZvXv3TrorAjSCr7okTSZsSZAr2L5NvQZrYASwFAAAIABJREFUE1jnSnWWW3ZCCpAf0XUzJT9sO7x+lQZzkQnGR37I+dQlP8SJuASWrDzD3hcGfj6lbXuM+Rg9xnyMijnPMP8N24cAobIwsfGi1i7bATd+e7TSwzlOtgwJj/ys1SUyuiawMCDHsB2emWoE81u+9Pebkjb+F14On+cB/B8AXwdwBTxSyfPfcf2yWEA4k3kUcKPqROHnU+3mrm5BfIh/T7JmLtlAztuOECBbJ2RbnyHZdoeWFtM2k4QoC61Gf40y3kKh/tCYCL2kgqI+0n3SKCwKCLItJ5F7hyI2U3eVCI9tO65g6pRN+zWZqj28604dS6D2WgmE/KzdICA97PlYDz/piTrQ4mgA+8CXSfuWPnklGki/WuERp0UALgPwR3jk5254WZZdgjwzBQQJTxTjpgxsm+ksVErQ0NCAz372sxnxqQZs2LABH374YQKOza7QD/xw9KRwKOQkyBXCDNa8ukO6afiZ/bLkR+V7MUFFfkj7MvKjmthXQUi0eWUmQoFRTJSqT9Q4u/TpESiiGPFJlK7aUyIzCwZK1B4WLkgPQZ2Y9BDIitWSNgCIVZ4oyA/JmRM2OKMIb7K/BcBFAD4L4CoAFwD4b/DVIhuTFhkHChp9igPtSDvZYTFlypSq9/PpFsTn2WefRX19PU46yVVm3aTz9xDyY1ux2CUI+UlDX3RgWYeIqD8q0kMibYQROewEoVAApJBMvIT80KSHm3TQlepDm6jOlqz3mKP9CTB1F8xMZYzaw3Pc5pIe02SIkmi6mcO9RZaJW1aktoxmxGPaAvhJAvuBn62YKD08NWgygCnwTOk5BEtJnFb6fNWqlx6irD9mg+oiOiymTJmC1tbWqvbz6RbE54UXXsDxxx+fUv8e23TpA6nPgfATjyjkWtl2pL4Ory9JQGfSEyVvU7wVG4UOQy/pnMrnQwoR+enJV3q45EfH1BNUTuRKSxoQwhxGyI+TxI2E9HCedTbTNkt+lm9iSI/o/lSUOVFurwtZCYgj4FUzZ5WAQunzcOq7yfDf96MRLBJKIyf5TYQCxGQninFQtF3yqo7LF/UTTzwRPXr0wAsvvOCszbjRLYjPunXrMH78+Ij3Yuvf4xKEdPCqG8eNsAQoruSPJtlrqWvMIz+iUgJl8iM7pskIvu3qgiU/pf2IJu0FuiRZhLqSslJBj/cOsGhHBj1lSL5fVRucJIg0ymU0WJioPbwq7yUIi8qWIFR52D7pkp4wIDlzZGPZ5+FVMv8r8/1yeA/HNxEkPARfKX3+D/P9/wA4AMAxgn2y4+4mhMtI7QomRCet/pF89OnTB6NGjcK6deuS7oo1ap74tLW1YfPmzTEQnzRiKJInQXH3wXYQsRgsafVHVT9p2ji1DwcA+0rqhPwwk6JUsaDX1VV90jtIdw1iE93JoHG8ZV8mNuxaFyLzVoOc9AyGhmmrmVpMYUJ2dbMiA14l8+MBLADwJICXAPwanvntl6iY6L8H76WITgb4b/DIzXUA/gDgLwBuAnAPgJmQ2wHTQnZ2oZpNWCYYP3481q5dm3Q3rFHzxOfvf/87isWiU+KTfH0uG6e+ofDCSkWpfeMwj+mSoCRLfYgy20rO+bb1+uHFQvKjmc1XCoHZiyU/o1DKuszCLLpJTjRkfj4VTBi4rrzEC8YPSBkSr5v9GfAIj8SnR5VY8t6tiJ5gqp4xGeGR9W0egH8BcF/p7+2lv79BrdNVWmjTVh2ABwCcCeC3AP4dwFMA5sArYcFiA9xmoTYBOf5dzNJ9MH78eGzYsAEdHR1Jd8UKNU981q1bh969e+Mzn/lM0l3hwNa/xwWGQaO+QYTIQ692jwlcThYmg2ppvyLywzpETxvNECBR/SYbAiSYnAn5oR13udFYKvIjNh+5N3eFga4TdZ1mVJpJLiDJs0tID++xu3crU4KEd1+sgl5NMFuQRIG2408DgIcBvAfPqfkFeGHpNH4L4J8InoQh8FSedfAyOa8AcGnpt0aIyU6c/jptpaX6iI7LF/bx48dj79692Lhxo7M240S3ID7HHnss6uqifIOqA9BBLexvaccwRENEdJHkvlXYAKNkaWyeFVkU2DSVz4QtROSHkyk5ZCi6mXnJBSpK0sTii5aKkV+NCjhoh8ooTaLkOBMnq/TQ/5QVnC3DNeGh++iiFARRumRkw6Ykjq1JzwXamCUqVMM8UcFxxx2HXC5XtX4+3YL4TJgwIea90iQoLfZe2WDETux5aonzbQrwIkbykt+TAJsMjobgOHVNXySUWQiF2UQIlvyUzjlvUveRH1U4OsBTU/jkQ9ROZfvm3HGKfenBth1apeox5mMB6dFVezhFTglE5q1hUJAeUqctKpVHRnh0nu9psL9HZQjjwxQWraUlSqJT3ejbty9Gjx5dtcSnumimIYhj89VXX51wT3jkxzblepg6NDZgJ+W4cmHkqb8LGutHdSvTgz8hPzrZmeFPOMfz7aH9NWcO50yA9PWcBvOEeU/Am4yZCWxqq0Dpoc/h2dA1F3UN2l9215hYfNEZmbGHrN8qUkfOASkRYEt6SlgwTFz2qZxosgHiF6RVpT7xShbsAfBf8ExDHwMYVOrHiRr9pe/h2wDcCy8c/R6NbQHUTQlRRaERfFJByJ3Ib6Sn5DfROdLZ7n3BbzLI2uweqGYH55pWfNLt2NwJsYSapHOvCrI8HlEhX1qOlq8WG4gCZEC22CgdXpCKT/3hkVibN2s9h9we79UbthtxIsJIwJi4GJ8kv9muDuFJj6S6u6i0SBk6Ks/vAfwNnkPw/4X3jNwJYI1kmxPgJz2b4ZWHIAkEWTBjUd0UbwEUt7/uGLYK0fst8fA+tWSwQTU7ONc08Um3YzMLmgSlJUpApiCNQyXyg1UzbM1cOtL60fBXflbtLwxUg7dMiudsG0hGZwtT8sOZwEsmHdq/JeicrBeZxWJi8UWr7WzhV5jsCJm9r9IMcK81myiSJj9c0kM/azwSwN7jL8ML374QwBcAjIRX6mE0gEcQTCLIEh7Ae/maD+Br8CcXZHc9xU94nCBusvM+0k92ojfAZA7OHmqe+MTj2BwVXIVLmvj32IImQWGdJHXBEqCkYOiLsJyN3mEwZ7jnEyIMe+4J81ITPPITvKfMyI+Orw67fbDS9DQsxzQsx7XFn3PW41eoNidXcrUnCB0CRUfAUc8YNzs2PPIjVXpUZIAea/4OoDcAVs2eDOAjANtK/55YWnjj1N3wXrRmg5s1+aDJ3iKDkeqznFqifDki2EEtNe3VkQiq2cG5pu+GDRs2YNKkSUl3wyFYAkQPsFH494QBTX6izrVBm95c+SDZmBsJ+dEMQ793a9CxmeU6c4YxkyXdL0J+dPPLLEVlsvYe/a5B+2MIQ7d2BmFwbvmvaSV/pxtzVxq3QpOmtTvGc9QeU9JD0BNYIIiSU5a+MPXfegeeTw/77kpyZNXDIzwivAGvEvrP4Y0d1D3Akh2RS44WTI8rLFyUGNFF9/bz6du3L0aMGIGXXnop6a4Yo2YVn66uLmzduhWjRrEVB9MOEy5Kv9G8A38m1KTAI2AkcmQyxIOxqwix4UjGD4mG7M2dOZZA7hYOpOoPYKb+BLMvsxN/WNUnqN7w8JDidz349yMiLOL+uyM98BzGecONj/SIVBBdkGu3G0Afzu+fK33ulGy7H56J618AlF4Me/cFDviUWuHRAq3sRI13mCVDnBg5ciS2bk1L8Vd91Czxeeutt7Bnzx6MGDHCWZvu7KNROi+7GgiijB6bSC1RQuSDpIKr68MSIAmBu3crcK/C4XWBK/ITnNxZAhA0Jen5+0zjTnbncr5zA7/iEzSJsZCbyHRIj+A8iPIhcZUecn+5JAYk4aXOvXs/vFCz/wCOOMFbVJCl36kDxGQnipQXGdFJC0aMGIHNmzcn3Q1j1Czx2bJlCwCPkXZv7IA3yNELQRrMY4QA2WYq1oUtCXIBA0dOkQ8I4TwLBop9SIxKTojJD8nHo+9HU1FvbExPplCrSUCFAIkJm3nYPWmLKXXBK2562Q6JeesJhCM9dQD6wgthZ58bovSIMmf+E8B/AgN+DBw+Hti/01uKnUCP/UDxI6C4R78rnSu9JVIUmMUW1eDZUV0OziNHjkShUMDevXudtRkHapr49OzZE5/+9KeT7koMMFUoCAEqINxA4gqfov6eDP9gHkUCRRdZak2hSkRHkdCb3laHOwfIDzkXM6BPgNTKjx9qR2UaegTFFdSmM5rIBUmPSu1hjp1HeAimtkJ8b5oUOZVhPLzoJDZ6i5gdOE7/E8YBI+oB7ANavge82VRZ9j4PdGwCdh4EtF8r3zUhO5ESngLSMz7JkObUI9FjxIgRKBaLeP3115PuihGqgQJbYcuWLTjqqKOqOKIrThSYf+ct25EpSDap6icz27n2GeA5YMcxkBHyQ8id4LwR8iMycS0YWFIVeH2eAb08NI9BporwExLyyM5DkJu0zoWInFyCheW/b8Q3S3+RhH7ie0qtLlX6c23x51guTANgSHrgpQHo2sjxsfFlfqadX10RHnLP7IcXtv7/4Ddz/he8oIexla8mUPf5p44HBq1k2iwCrd8HunYBB/8B+IRTRPhDVXLBMAkEO8EnON3bebgaQFxJtmzZUiVpYzzUyswdwJYtW5z698SDtFwOntogS6wXF+iJKyoS1ABxeQpbiK4rmUxkeXl6euRGZN4iVda5pRZ0yY8fzbnjfOqIn/ycDZW6cmPuSonaQ2co9vp8ZO55vFEURV+S9c8sr28LfiSYBekpOX8HyA/3GrgmPARfgOeYfAM8k9dhAJbBu6duAdALmDAaKHwLWHcPMPYNoNdhQF1/YNgXgBamuR79PXNXwxcAYrX4UNM8a43XIm4/Q9QYOHAgGhsby64l1YK0zLTOsXnzZnzta19z1h7fLkrecKrxNJoqG0RCp7dLwl+GgH7LZUmQrXmMkDra2dM1CeKB9F9CgIi/iJAANQkmXjJxyyZ48ltlkmfJjx8i9Sao+rxRHIUjc8+X/sWWZXgcHqExA1GIPHVIRID9ag+BmQ+SmPSU/z3mY3QNitK/Qeb3tgDALwD8Bl7uniOBsX8EBpaOfSeAYldp4eTp8SGHcubmyAhPGoiOTJnKYIpcLoeRI0dWnYNzNc7YSuzduxeFQiFGxYd9kMjDxTu9cdqEo768dBZiWxL0KclvuuYxljCstOsKF4QEJUGAOPcKTYDY0Gnic1ImQKa1t/xmL0J+pmE5phVZtURsuvLW+6ZExbFHsM1p8AgUYK5GmtX04uU78nyieJOpudLmB31Pi+p4fQrAHG/h1YIDgCP+4C0sBsCv+uz5hfe5bT3ktcNM62VVX6hz8nBJztjac+4xYsSIqlN8atK5+fXXX0exWEzQ1NVJfbJLO9xWbE+Lc91WeH4y7ALY+feEQCQp9m1D8G3Ip4YZ7zKJ2eeZJsF+dULS/WTAPPJpMcQmKXN1R4SKikSDPFs6as9DEDtnm5AeHlySHgVmjhaTHlFgF8G29ZXFGbYyS1ikZXxLO+oECx8uI7uqkfjUpOJDLkK6fXx45KcO6QgxF8FmENrAbBd17h4GNPlphETGl5139jGhj8GgVIURaPIjyNFDlB1OdBGZqIOTs7nyQ/vsXFv8uUD1Cd7PxHfnjeIkAVFxgcc533nn5driovI3QdJDg1Zng6SHmPzW7vDKQ/AJTx3CFW6VkR2OAsNm/NbF2iiyqLuoP1fNiNMJuyHGfelhxIgReP/997Fz504ceKCKbacDNan4bNmyBY2NjRg4UJTvJM1oFyxxIWou3Ax/bSsbfwIJAVN1n9QfcpKhFgiG35t0RhdPwO8kyxz/1Faffw8pOgqIalKdDbX64ycHNHHwOy5PA993JylUVCU6WqwCkXkumPNnYvHFMunZv/sTHHLTv6Nr0HB4UVOfA/AotbaI9KwF8DsAVwH4TunzdvhLK0yH8qVi5nD/oosD4ZEdsvgguz9lLwI94ZEdslQb0v6+3yBZ0geSK6+aVJ+03wFWePPNN5HP55HL5ZLuigFUl0JEftoB9DPcV1rk47+XPlnyw4axRwSa/BwIB5I/TX6icBB9AtIMzVNbAdQBrAOusfpDzD9+Z+VgtJZ+lBVf9WkHmTiPOH976btN1G8E55bbAGgzl5hgkXUvwUJNtcnvlM06dm89+4fY9fRr8Eo9HA3v3HwLwBoAJ0vafQLe8zkDHmFqAfBnAPPghaJrJFgl5UpEqZ0GA9jO+X65S1LiQinKHIv9SCeRMQXJlffWW29h4sSYFX1L1CTxeffddzFkyJCkuxEjeJXbTclQVLAhWasQvDUN/B5sQafuL5Mg2wgxwCNCrgkQrfzwExWaFR5lyQ87MfHJjz9a60wESYgXseUPVd8EjzS8rNk3ApYciRD0Iar0sRU6EWBAkPTsXLYau55eB+AuAOeUvv0cPDXnIXimT5F4/n/hfxanA7gMwIkA7oAXlUXAmEyk9dkESB3ZcYFqz+eTBoITHelsampCr1698O6770bSfhSoSeKzffv2KixO6hoy85gg46w14riN6Iin5eASIVk3ZAISzyxNSNAAWPpFnFL6jFIFIk60/orrQEXdoQmQOHMxMfGIzEB+8kNMSH4Vh0d+PMgJy8vwJdvTAq0KBSeVoDIEiCO//KSHTXTonae74N1AZ5W+JUTx8/DMVm+AmykZgJ/0kOs0EF5ldc5EYUN2ykVuwyQQBOKJXKxVyF40bUhb9ahjuVwOgwcPxvbtqheT9KBmic+pp57qrD2XHvDpgMhEISNEaTGPPVv65L29R6QK0ZlvgRAOonSmZpcJGJdCpv4AQRWDn5FZFqL+EMj5pVUcuePyQwDGIrdkGYrnfRnF876M3JJlqoPRBFGNCMH3m8P8oAkZTZaCpAfwEh36naE3ARgBT9Wh1TFCUt6BmPgAwWvzj9I2X658RednEpm0hlG/3esqTFzmnF/tSksUkL1BVQdRiQIZ8UkYxWIR27dv72amLleQ+Wxw0thHBsvbkq0d5DScncKEcR75mTY6hGnBdRZqOoQ66LjMS0hI/u0nQCLy0w464SCf/BDVR+WMPxasuYuQIxXExMl/HeR+Pe0A/s33jbyuWCu8nDmsPxTJ2rxbsB2PjHYC+AG8IqPfkRScFSA04YkqClEX1aBkyF4A00AE00dIhwwZkpm6kkRrayv27duHwYMHJ92VGoOo0rRt5FwMChIhQnVTKmHsbDSXLPpygOS3aSUViM2hsnwTzKNlCAlypQIFMzEDFYLDEiB/iDogV3745McDz4fHM2XpEhsR5NuOraxzPpBbso3pLwtS/sKb4NSZnXdDnmiTBTnv7ATfBeB7QF0zcMmjwFjNHEm+grWy5IIiUhF16YlqRB3EBMc1MUsfUQE8S0axeL2TtgYPHoznnnvOSVtxoOaID2GdmeITF0SECLAvduoYndTAz+bxOdBVWHsJ00Z7UTaAxds5HbHlosYTvwApTYBI/Soy+ftz9AAeARKXmiDkp3jeEZQaE1R0CNyYu1TO0UT94ZVIqDhAv1EchYW4pPzvIOl5CN4Q2RdePSwW5Lu+pU9ZioAicPr/Bzz1CHD1PcDUM4FXBasOA3CZyN6lC1dkJ52TthlEL2fVflzpQab4JAxiZ8wUnzSgIPg+H8G+LG9lXgg7Hd0VBmy+FSMi5IoE0eaZoAJESlEQBNUfuiQEDZr83KXw4RGpPoRQeeaX+Usqf/vhfcdvm1J7fOsQckQrbDTpocPdeQyEKF6d8JjIGniKDR29RcjJhfCitFjUAc/08+pk/fwS4Km7gSsXAVPP56yLSjkSKxCik3Yzki10CBjPHG9zPqrBHJcuDB48GC0tLdi3bx969eqVdHeUqFniM2jQoIR7YoKauwwKFMA3deUV28XkYL1tvUd+iBMz69wsM4/J+PYcigjdxL7Ry46N9hUJUw4haAIjJIc291TID/H5EkVsPQ7i1Bt0YOarPmJyI4d/O354cJD0ABVn5iDpAYgvUMXsxTfxHQfgOXjh63Sekq0AhgCYENiinDvplSLwi297pOeKO4AvUb5FowBMFZEd2URPzF3d0YSVF3yfqTdJgggN7733Hg4//PCEe6NGzc247777LpqamtC7d++ku5Ig0hKBZYqC5DdZ1EwE2EZFbrFRXNMYImQDNmz5Jt23fRckKGgCYwnQG0UvHYQ/XB0IRkhVwtHl5KdCduaV/tNFcF2iFHkkRO07VLl3gqSHQFxwFTgGwGcA3FPa9zkAHgHwPwAWolzVHJcBdX9Ebs0rKEd8PTALeHIRcPo3gSOOAf53tff9//kIQD3Mivvq+IDVilqR53yXBnJTzaa/6O4N4lry7rvvZsQnCWzfvj0zczlHGm6T18Dvh2XNIgDWicV4kVyiQpG6oKN7tE0eM2BfI4pM9P6QbjY7czAnDi9yS5aLx+8fxJIev3Pl9WBdLeeXsq/ziVJz+Xu+2gOqX5vwRvFb5W/1cvzQmAVgI4A/AVgCL7x9EQiB7PHeAej6vwAe7vJMWwRPPwnkch75eWIR0+bh8ExoIrhMeZAm0BMw+0LjemKuZiJo2ndy/8q2ieZckDm3WkLac8Ui/ZRWP77+9a9j586dePrpp521GX0eH9fEwrXik5b+mfSDEBHZvkwKk9KQtVnajvXvkeWlk/0mq8LO7aMuCWIHQJoAeftkc+JUyMIm8B2MWX+bedw9h4kkYZ/FCpGayOlThYyRPm17YLCA9BDw7olzOd/5o+NIAVMWXh4l2WTDqgcunNqjmNzCqBy8lxNZe6L+y7Zx3Z7JNQu7LxlEY5QsZYSs7/7tXEV1dXV1ob6+Hrfddhtmz57tpM0okYZXeaf46KOPqqZCbIYoIcqv48BMJQT1OAUcmRvssvKyldenyogQ4DdhiUgQb2AMKkCEIPjz9dxV+pWQCl7o+jzQ/xG4GGRJG7ncfEYFaoZ/kgiSntySZcASWplTFVNVEx6CCQPXSaq3y97ceyKc31aawKqeaVdakjRb0ZFmMhJTHWa1Hj16oH///vjoo4+S7ooWao74tLW1YdgwiwkmQzcB7a/DKjc60Vwh1DTWoTkMEVISIECPBLFYXPqsOAPTBCgYPcU6MXu+PDTxuT4CUblYvB7BGsTtoJ2Pad+fSn83gR/mTrfxb9xf2HIWLIKER4YoyU7U5h325SGuyTnt/jXDIO6fKqln9aOxsRFtbW1Jd0MLNUd8du/ejcZGd5W9a69cRQYxRLWKIqo4XCZCFJmifX0k5eZ6vFcf+E4+8arqcbEgaghNgO4qkwl/6DohP0F/HldSOg+0+lNBM4CJAtID8MPcCYJFTgF5Vudg2Q8RqlHVoV8Ean/i9oMmj3nmN5tzIUs8aQPX7YVHY2Mjdu8WZTFPF2qO+LS1taFv377qFTNUGZK8VUXh1y6SHzIKks+xuWfQ1CUBCaGeMHCdZEKmzTcsCeINpH4CREgEid6iv2MRJelh90OTn3ng9UuW9JAlPK24tsg6InsI1vGSwcbpPG7QpL47ORYTs5yMPKSLWKQdffv2zRSfpNDW1uZU8clQq3DhAM7mUaEfJ5oUWUaPAXyTlgYZYv1Q+ERIRoJYkN8roesAqPD1eb614yI99P5ovx9/f3ikpx08Hx7iz0RndAZEpSxYpJno0PejTqX2akNPiFMDZAQmDmSmrgTh2tSVIYMe2EdJlFzOQQX5qa3oovZHlB7AU3t4oMtT8Cdxul8yp19CIioEaB6+XCYcUZu3ZCgWr8f8HO1QzTo8EwRD73lV3auT7BCCk3YCY+KvQ5O2TyTrVcekW6vITF0JYd++fdi3b19GfDKkGKLcLNMF3/Pgf2xp355m8AuR0mB9VoITPG36EVVbrxCgNJAeAnliRD/hKZ53BLY9EMz5Ja/sTpSvJIdOmqTGpWa49imZBjHxyQhMXHBZqLSxsbFq8vjUFPEhMlvm45Oh+iDL3zJD8hsfrGlLRYQuwUIAvEn/TIjz9gBsRuZkaU/Q34eN9AL4mZ71yE7cIGTYJpIpIbNV3RT/v9MuPKUCrkllMtc+8/FJCOSkV5fiU1OXIEMkEEUEyaqB+9GcO66sBhEQ5YeQHoBX0oHkQ6LVEkHl9YTVHoIg+Sl9zxCeYF4fAp1Mzi5BiG1KQ7VJ0d6dgt9TP9fVfkRVGpD5+CQEYl+sLuLjGmnP2pzBHVj/krPBq8MlAjFx3YhvAgj6ubxRnIQjzq9I1/7wdcCm2GjcIKav68/z/h2MQNsEzymW59dECpy6ALkmKZBAJozz6s+xxXcBoMWivUY4Jj/OG8wQA6rJx6eH7oq7d+/G1VdfjS9+8Ys45JBD0KNHD8yfz3mrKhbxn//5nzjhhBPQv39/HHzwwZgyZQqWLeOHvP7mN7/BqFGj0Lt3bxx55JG44YYb0NkZHBzuuOMOHHbYYRg8eDBuvPFGblvxKD51mksGPqq1gGo14DHqk10egspkc2Tu+cBCg0Rw8UxFaVF7CNis0bklyzik5+XSosrgzIOIwJwtWGLCzOGVZdpo/nIg+KQnQ4YQcKX47N69G9///vcxdOhQNDQ04Pjjj8eDDz7oW6ejowOXX345BgwYgFGjRgn5hQjaM/QHH3yAO++8E8cddxzOOussLFy4ELlg6lTMnTsXN998M2bNmoVbb70V7e3t+M1vfoMzzjgDjz76KM4666zyuj/96U/xox/9CHPmzMEXv/hFNDc347rrrsM777yD22+/vbxec3Mz5s6di1//+tdoaGjAD37wA4wePdrXFgB8/PHHAIA+ffoYnYRowJ5a0YRvS5Jk7WVvSxlE4JGfdoiS97FkgZAekiungnQRn+uLRcbcJcvjYwN+OYtIsGAY8LbgN9H3GbohaDIev9m0T58+2Lt3Lzo7O1FXZ//yf/bZZ2Pt2rX42c9+hhEjRuD+++/Heeedh66uLpx3nifd/upXv8Lzzz+PBx54ANu2bcOFF16Il19+GUOHDtXah3bv8vkMGCgyAAAgAElEQVQ8PvzwQwBAS0sLFi5cyF1v8eLF+NznPocFCxaUvzvttNMwaNAgLF68uExWWlpa8JOf/ASXXnopfvKTnwAAvvCFL6CjowPXXXcdvv/972P0aM/+/pe//AUXXXRR+aDfeustrFixIkB8iFIU5qTXBkSKl0ZxTSfbvCPZJkM6ITL1+COhKkTI79ScftiYrc4EoJ9A0hhsPqZXo9tVhjiwC/aFT0U+Q7btxQ8y7+7fv996Dl62bBmWL1+OJUuW4F//9V8BAKeccgr+8Y9/4KqrrsI3vvEN5HI5/OUvf8ENN9yAL33pSwCAp59+GmvWrMHZZ+upq1a9kxV0b2hoQP/+/X3f1dfXo76+Hg0NlYHnySefxN69e3HxxRf71r344osxd+5c/Nd//VeZ+Bx11FFYuHAhvvvd76KhoQEPPfRQmQTR6OrqAuAVTMuQJIZCTJhsiJRsO942ogKlGYJQOWmqlZIoanHFi6MhTn5nBhI9J6zWvrGfk/1k4GGH5DfRfW5bhT1zbmZB5l0yD9vgT3/6ExobG/H1r3/d9/3FF1+M888/H2vWrMHJJ5+MI488Evfccw8++9nPYtu2bVi1ahVuuOEG7f04l0auvvpqXH755Vi0aBHOOuss7NmzB7feeiva2trwve99r7zeK6+8AgAYO9b/Rjlo0CAcfPDB2LhxY/m7r371q3jssccwcuRIAMCMGTNw6aWXBvadEZ8MHsgkJiNZord/G/LFbidKXpghTgSju4KJC3XBq9clK1jabdHJZhrfJVnZplxESiPfMjghPq+88gpGjx4dmMMJT3jllVdw8sknY+7cuTj99NPL/sY//vGPMWbMGO39OCc+l156KTo7OzFr1ixccomX9r2pqQmPP/44PvvZz5bXa2lpCahABAcddBBaWvzhBffccw9uueUWdHZ2CquvZ8QnQzqgW65CRqTC+ISlLZtwOkE7afMSGQLB0hWpwzD4/Xzo/Jii8HPAUrCwISr9ICc/pnBdod11zhvX/Uv78Vbggvi0tLTg6KOPDnzf1NRU/h0ADjnkEKxduxbbtm1DU1MTDjzwQKP9OCc+t912G6666ip8//vfx/Tp07Fv3z4sXrwYM2bMwGOPPYYvfvGL1m0PGjRI+ntGfDJkAMyiiDi1wMqQmQzTH8oO+P2Prg9ax5W4BAuNyE/zz77g/fEPyUoih+QPJNvISIzstwwZYoIL4mOCXC6HI4880mpbp8Tngw8+wNVXX43vfOc7uOWWW8rff+lLX8Kpp56K7373u3jjjTcAAAMGDMDevXuxZ88e9O7d29dOa2srTjzxROP9kxPOizbLkCFD9wNNfOYtmShekZvIkMDEwVlH3RC9ccve7G38TWz9V0TbyRQfm/ZsHXdF27luz/X5s1VuXLcXDVwQnwEDBgSsPYDHCcjvLuCU+Lz22mvo6OjAxInBAWb8+PF49tln0d7ejoaGBowb5+WR2LBhg2/99957Dy0tLTjmmGOM909OvMz5OkOGWgevPMU0QY0wOmszCzp5IYv5S+h/pSuUPUOGDPHDhcVl3LhxWLJkCbq6unztvPyyF2Rhwwt4cEp8iO/N6tWrMXPmzPL3xWIRq1evRlNTU9mn5/TTT0fv3r1x9913+4jP3XffjVwuh69+9avG+49basuQgYuZwyt/893R1L+N4n/dY8zHwk1EldldIJj8r2Lqqh7aI4tQe03wvcxHizd8qsyMj8Gm9lqGDGmHC+Jz1lln4c4778QjjzyCc8+t5Mq6++67MXToUJx00kmh+wkYEp8nnngCH3/8cTk748aNG/HII48AAL7yla9g2LBhOPfcc3H77bejV69emD59Ovbu3YvFixfj+eefL+frATwH5uuuuw4//OEP0dTUhNNOOw0vvPAC5s+fj29/+9sYNUow8kuQEZ8MACqFEmUJvEW+cDIlVeY/x/eNTTVuzF1ZLlfBR5YWwBw6juWi2msuMFm9ihTNAE5w0ZESXDvSujbvuG7P9fGmvb0KXBCf008/HaeddhpmzZqFXbt24aijjsKSJUvw1FNP4f7773fmxmJEfGbPno1//MPz2svlcnj44Yfx8MMPI5fLYdu2bTj88MNxzz334OSTT8bixYuxaNEi9OzZEyNHjsQDDzyAb3zjG772rr32WjQ2NmLBggX4+c9/jsGDB2POnDmYO3eu1cFkxCdNEPlTfEqyjYypCKKcunuuSg7YyuwAAgVK9eE643F84BUqrX3QqRRsJ7n1gu91SALrK7UVwHDeihlqDK6Cix577DHMnTsXP/rRj9Da2orRo0fjj3/8o08BCotcsYYcYp588klMnz4db731ljDk3RTBwTOKEhOu2xMdu0y2F/1msw0gJjExEp80KD6uTF1TSfRVFExvseB7/fDltNXqAoD5uRyTVVpExkVFSgFzU1eakHa1xaZ/qj6IShbE5WBt6xDtsg9m7bl6du+9915cdNFF2LNnD+rr6520GRXS/uQagaTJ5hU5jR/sgCkaQG2S7Nm2lyFVuEyWaVZ0HTth/9jKi5QG0YBqz1DLVmcHeP5KAL9WmU3x0gzJwrRcjsu5ohFeVF/3zM5N5t0DDjgg4Z6oUVPEhxQnJcVKo4FtyYUMVYEWyFUfHdy7lfkibpLK8zMJO8AHsx77C5WmS/HxlNp55X9fT/3NVpfPLdkmaEVUq8uUQGboHiDFoV0mbDR5bpOdfz7++GPU19dXRa3M9PfQAI2Nnm2DOF9nSCM+gdzcZYgwAogJ1m5AukhvHNmZzwTr4EyThvlLKorKvNz8VJq7CERV5gHgjeKkwPpH5p6XtCbyNcgyZmdIEslaOtra2spzcNpRU8Snb9++ALo78emAWyXBpqp1FWEbx5FzW1yPhc61chUBpJOqnmfu8RBUSYJV2tOLBlSctD3lqtL/0XjjgSDxeaM4yZex+cbclRr7kYWyZ6QoQ22jra2tPAenHTVFfOJRfFwTAddEJQM6V3qfH0rO64dpI3NPJLRfol4ES1cQJYROZMj3jyn9lhLVRx7N5SdAQFDd4SlAbJHS5ZjGjZ4Tg5AilnxGGdqeIUN8yBSfhEBO+u7duxPuSQb34GceBgB0pkmhkUFyDLFBHhLKm/SDZIcOcU+XA7QXyUWXqpgIfkg++e41sEpXhQg1cauyE7AZstVEiKe4kWSGbIROWCKchgCPDN0Ju3fvzohPEujVqxd69erVzU1d1QY674jsdrQhHEmpaavUq8QKeTbha4s/55au8Cshm6CT02d+LofrE8yQMY/5zwNRd9j+E9WPjt4iJMhzamZNXNcWf44bc1dyy4JMLL6ItTvG+77rGrTf8AgIplN/s6QoCQKdrrpQGdzApUKbKT4Jom/fvhnxSQVkvh9RhGonAZbgpKnv6grttJqxEJeUyQ/fsVdEegipaC4TjaRMXsTERUhP8bxSf8qKlYgA0SAkqAE8dYwQoebccVzyw6LHe5XQ3q6N/ah8TGEwrfTJKm1pI9wyZIpUrSEjPgmisbExIz5ChPVP4pUwEN1CafdbMjkX1eDAqyY6QNBXhQaf8NBEgJ1oPSJBHJ/pKK+4yQ/r1zMPyzBvidc30j8/AZIlLSSgw9aDJIg1bekQITzDhMcriZCJ0kKXq6CJRTURogzVira2NgwZMiTpbmihJolPcj4+rhWLOCKqCqVP12Ym14jabCVK059m6JEd2YQcjFYiE7GMFFQcg/2kwk8Q4yI/QWdmEsU1Frkly8p9DBKgM8E/Tt4zR0hQHUTn3SNCnmmLVnoAoGtjH37nn2kCXmW+kya2tAFbv4smU9VA6jNUAzIfnwTh2tRVLF5fpTV/3hZ8HxeJSWu02gbq7zT2Twb9qt4isrMc0zANyyXh2TLC0wDg6PK/gkRiIuIkP7nc/LInT8WfhyYtFfJD9/eN4rcodYv49JhkaSah6WLiSfv2sCRIiQUDK3+/DeAm0bNsA1ZBost40CpRNb4IZEgSWTh7gqhtUxf9Jug6mV61+dfQ4PW9ViqL65MdnQm2OXecoGApUTREUVqEIGzy5fSpkB7iN+Pffh7mYX6u5HPjkACx/jzef8sQ9N/xyE+lr17eHhK9FiRATdDPzEzn5lGRIMqkxZq7VJjDKeomJUMu/GfYCu3kmm5gV4wRmYN1mpH5+CSI/v37Y+fOnUl3wxCyFOcu05/XGtjSEEnDBXmcrl6FgYrwyMOsdSb5Srj3G8VvAed7fwdJD0EDgLHlkha03w8QLpKEp74SR+Z5SwCP5PDID/nNIzuE+AQJEOD58/DOi+zaykgQQ0R8fj0dfoWHhozbzGTIUKBMSlQYx/ybJiPV9rKROVi7QldXFz766CP0798/6a5ooeaIz+DBg7FpUxoeQPatWfaQ2VyGtJqSdGDa99dKnzV3u1IQER7BuSqpBj3G8OvSrd0xXhFKLSM8RLXx57chJGHbA4NxZO6u0re8CCmPYFx/HpBbMi/wKzFRAVCGvtNEhxi1gvDMNcXzvlwiYzzyA3j30WgAfvIDeMd2ZI52tiHOzDZ1udSmMB9onx4RCVJh5vDK3yTf5HJX46BunqbRzL/JmJe2F5So0H0VqZaWFnR2dmLw4MFJd0ULNTeTDB48GNu3b1evGAqyG1z0W3ciKrYoxLAPGkldk2nqVXjQNJHICQ+ZlFUlLPikB6DVER4p8Ed6VXxIgk608zAP8wQKzjwuwWFR8echTsxi8kPn7DnTdxykPAWJePP7P9HRXKZlJ+j1NU2WNAmaY0mCCKYxRERKhKJUP4Yz/yZj5Gvsihl8qB5Fisy5GfFJCEOGDEFrayv27NmD3r17J92dhJBmf51C6bNaSaAtJsM6Qk+D8KgT5elO2vRE3xrI5FwhPY8j6A/kJz1+Uxg59nBZniu+PKRd8jlWQn7Y814hP+SYri0uKv+bT4AAT8GxrblFl6bQJEGsHw/P14eF7J1vAqvIoFR8NykczfmOTPaFGPuRISzeffddAMjC2ZMCYZzvvfce8vl8sp3RRpqJigqivr8T477SBLqPbBixLD0Boz7RJo9XAYzib9W1sQ8wVeYHJpuoWdWHzVVTIVx+Hxg6AoqQn7HggzU5TSwTI88nRxcTMQ8eseFvJyI/54IfsUW+OxNAU5nk0HmOSJZmP4j5KkzRUZoEyUycDMpEqN1v2gqDCSWfnZbSv3lFe60QVq3Ic77rQEaI0gmi+AwaNCjhnugh7bOIMQjx2b59e4TEpxomYBFcm3dc5xypZhCiY3lvmPh3+BxkefszmZiDyfkIAViIS3Bj7pvUL/K8N0G1hwZLjtg6WhXfoHlLiIks6KtTITVyVMLWRfl6AK/8Q+X4b8xd6SM/JC1A0EHcBQEC/DW5DBzbaWdmVyQIAI5go7kAbBMlQIzTp4XsKy/4LYoXrQy62L59OwYMGIBevXol3RUtVOvsLQSR2oj0lsEVMoLDB6vqGELHfEEgzPJLE3HTidhPetjMzhXFoxVi8uCZjYrnHQFAFu1VgQ5x8aMSll4B68dDVJ9tvqgtMfkhhO0hsOQHQMwECPCTIAM/MF9EV2fQt0eFFslvOwEcxNzjH6YxE/RQwfcZIYoD7777btWYuYAaJD5NTU3o1atXDA7OtYhO+PKN+OD6VkmLs7dpP+iEbzY1x9qBOYI39LcB8HhQ2eFV1k/TiTcYcSQmPYCK9ABetBcAwMiEZYagKYuFF7nFhqwHyQ9rciSO2n4CxCaCnFh8UZAewCUBAiqFSFkCJPORKpmXaCdmQoJcZvggRIhOl9a50kHDUTjzilTUDgDvR7C/6oHLvFrbt2+vGsdmoAaJTy6Xw+DBgzPFRwoXhRK7EyaqV1HBxhyhVbpgqXqVANRh1n7S8xD4iQ3p/D68fDhBMxad/NANRA7MwZB1D55PjzhM3a/+EJJDEyBShZ3vUB4VAQKsogHLJKhkKprA5uFxhLop/n93iirIx2keUxGpQwXfd29CZIN3330Xo0YJHBFTiJojPoD7kPb0lq2QqRWus1enxa8pjn7Qfg4h92VKeKRkh77eNoQH4JEeVtUIkh4CfmRWkFwAlaiv4Pr24Jm7RiNYcNQfsv5GcRKl+uikBHgM7HniVWMniSODBKgTlcgt2+vEA00mLE2sdBSXigTJVCLZ8NIJ8EmaiAylDaJ7RPbC2H1z+ACe4nPqqacm3Q1tpGEmc454cvnYIArzjiyiJw2mpGpAB4CT3DVHzAsi5Zdn0vKFLsuuW3SEZzmmMWYckSrSAN7E5g91Z9cfjWBmX7EPkAzEh+eNB74lqChP+uDPRXRtcRFF6kTZmQnE5GftjvG+73u8d4AknUAUBAioVFwnBEimbggm5bUb/NvxHJudYhrTFxdEKE7CISvH4FpFr54cPsViMTN1pQFDhgzBc889l3Q3HMPVW3MGD7qSv0xhooisqUMpgXYByifUqwihJj0A67grIwXEFNQqMHGx8AjIG8VvYRsGU+uKiY/fh4c2Y5mfZzpPjz9EnUd+6GsdzMBMKrCzZULK6s/GPgIn9KgJUCdCO9rToeyRkyDAT56JhMRznE6TeUyEfoLva7/kUGtrK/bt25c5NyeNww8/HIVCAcViEblcLunuaIJMsN2J4MTt4ByBf0MkhIc+L24JDxAkPfJaXiwq/i900r8KTKqcs3453rmsOEnzzFgegiYsXtQWMbed6wtTF5Mf0XAYVH+6BgXJTxnPNAnITwcqIethrisLMlnTpMEFCaLGIjayS9UVLnQIDLufNEaQmaAf+MdNcnpV/3j/j3/8AwBw2GGHJdwTfdQk8RkxYgTa2tqwY8eOCBMqhfE1Ed3sNXk5IoDJubckJjLQvhE7ARwoWG87/OauQCFJWSbnsBNjkPSQiXrtjvGYMHAdAB7p0S1p4WEhLgnRRxsETVhi+M+vnPwAcodkPvkBGLWHgGTbFqYgiIIA0aAJg8g5X3Z9mTHKF8LeEXRmjgyTAXzCfBcsf2KGJPIPsWhgPmlUFxnavHkzAG/erRbU5ExLLsCWLVsSziTJG1iqx3YbRFocnFUQkR2ZwiTLqtwJTAgh/RtVzg7r98Avh8CqE/yoJHbiF2d2pkPfg9mNaRAzF88B2g2Cqo9eaRB+ZmYZ+MVHvfPYCTzTJ7BFRf0RTYBREyDATxQcRCgClfD12AgQDXIMxDzGK7sR5zjrel8NgjbTSYi2bNmCQw89FAceKHoDTB96JN2BKHDUUUchl8thy5YtMe2xU7BkiAfDmcURjjihsthg+SbgXllhSHogWw63pKdy//FMMmrSQ0CIrg7pIWHvYuIh9wUyQ7CtJgD/xlmz4sdD93c5pjFmP81q6qJzJVJ3nmlSZOXugOfvYlO81nScaUZ4xYTe/Ur/YgVZiBir9vAwjlmiQBqithrgPY/sYg6XOXy2bNlSVWoPUB2v78aor69HPp+PgfikndykJUmgDDZ9PNpiG02ISE4LgAGCbWhzl7QCNg8uIltERS870eO9+sC3+qSHIFjSIgjWSVilupj4AvG29ZSkiqlNFaZeyc/DS0xYgW4hUtr0RY0DU1uDRWVfLX0S8iNNWUDIT9Sh3zT5Ye95mbIgIwCdEOcdipM4jEPwGEyfyyRhOq/IpvHo56gtW7Zg3LioCGc0qEnFBwBGjhxZtj26gEuGnMEGR1MLoB6AReBsd9Dk0hLCnLV8kyHpWQU3jpsC0vNME/BMk9/vBDakx6+C6Ju4AI9oNAFowkJc4tAfqKHcLq9fOkSN9m0KkiAT5Ydz/qa2VtSfV4M/+9Uf0X2sowC5mtTWU4tLLGeWJEEc5fMQq8Jpf5G1gb0qpINisYjNmzdj5MiRke0jCtQs8RkxYkSMpi5XqIYHL64+5ksLTXYcgxCeMFi7oZQPRZeIuSI800sLZ7+s4lCCOenxY2LxRSwvTcbiJIdB+EkJISwNUJnFKuv/m3B9Ox+dIOzJDyAMURc6NsMjP8qitO3wnHtD3qNGWA/PZ4bnNxMWLkmQzDym4wvDmscdmsjLSIN5DIhyzN6xYwfa2tqqztRV08Tn9ddfR2dnNZCJDB7y1BIR6qbICY8sWy0p5kjIzlqTycEV4QGCVbypAZZDero29rEkPZXJnyYGy7V8UXTMYwSE0LDqkMx8FSRbctUnWMSVjWhzT346PLOWyLT1NrwitVqFalkCFNW4Rk/WLAFyqLKWCdCTED8XOv49LpFHLGNQDYGICxnxSQlGjBiBjo6Oco6BDGlFHl5l5bzgd4uMtOx2dVMqC2BfzWPbegXZ4fWHEB4XExVReQT75ik9ZbMLvX8z0sPCLO9PBWbqjCvIpX72WIKO4Lbkh7kXWPLDpnEKkB+RahG3AgRUCFCUfjKr4PblwAXypWUgvFTrPIIa54t1WhSkCrZs2YJcLoejjjoq6a4YoaaJDwCnfj7VibQ9LHnE9lZFk50w2LaeymqrO9C5HsQlKovIbBIwtXTCpnim2BEYpfbE5IKt+G4C3W1pQqXn6+M/B4T8kFIU4cmPwPQlU38AA/WHZGl2TYB0xopN1BIVyLPjOAJNagLTDRUfxiyuET2RcumvunnzZuTzedTXB4Mo0oyaJT6HHXYYevfu7dTPJ3NwJjB9OIdCruq4xDRoOYbKVJ+dqJCdbTYOnyLCYzOo0cfCmZhowkM70nL9S3TLJfBNXACrkPDy/tRBbeaS+wTZYjmmaZrhABH5IQhHfgBpXp7LFGVKZg4zKG6bhAJEoEOCXJjHmuEnQbIXirhz3QzlLN0HW7ZsqTrHZqCGiU+PHj0wfPhwvPoqL6wizagFnyTRQCAbBMMMkAKyY3MqP1zFZKk16U8z1CqPSacsVJ5X4Yz0APAV5LQ1cUUHfn4eP0wJSwVB8iNKGyACj/yU7vOb3uaXLaG/4pIf0f2TFPkh2EotUYKQnyicr11hKDz/tIGlJUlEO59s3rwZw4dH4RgeLWqW+ADAuHHj8OKLMpk+gzsMRLwPuqayo4sA4TF5c3QtycuOq0MeDXTZDgRJpB3poUtcBKFnMiMqjDv/HnXyNrlpjob/GHgV1r1zQCckdUF+KKiK1M4cbqD+TIR9ZmaXJvE4CBBQ8T2KKgoNsM9pRGMgs4gKmlYXdu/ejS1btuDYY49NuivGqMkEhgTjx4/Ho48+is7OTtTV1fShKhBFIsMk3mQsSI6sykYbgE5bP5xO2Oc9kXVqMqTlM+YM81QBnnuBz3+EXPNwpIdAX+0Rm8kqkD+LfHPVubDxT6r0SbStvwYXW3zUy4O0l9lmBsyqrBPyI7h/CfmZKfEZmTncq/2mlSuKkB+XZNwGdF/ZlBS2QQuy3+jnkS1b48K/xyVY8lN9VdxffPFFFItFjB/PezFKN2pa8Rk/fjz27NmD//3f/026KzWCQ6nFZWir7Dfiw3CSZZscaKfYFw2IUSR70/DVoB1faaFA6DSbNsd2GWxNUvzwdH3Vx98GUFF+yskfuXmRTJUfQJ6/pl1e02176XPaaG8JgHetdRWgOO6T16glDmxCPI7YrtAP3stOI7XQcHONXPqprlu3DvX19RgzZoyzNuNCTROf448/HrlcDuvWrXPWZjwOzmny86HJTpyIwGmTS3ZMBhSXhIe+xrzjpEiXKNrnbUiihJ6AfuHLGb7+sGpP2MSH6YAZseoaxKg8QvIjSi8ggiJ5n25BWy75EeEEBMtSRA3ZGPYagEJpiQtx+SAB7tQlmgTpFd6NE+vWrcO4cePQs2fayyIFUdPEp7GxESNHjnRKfGofTcwSJybCnvBICEzn8hAFFNsRjcIDVMKSJfuWhTcL/UNMKn37C5uGIz06Zq54sHbHeG5xVjHIcUkKDAfIj+1bOJu9mJkM793qJ0DbwUdZ/dHtB48AJa0KFuAnQDJi4kplpkmQbhFhk75EhQbOkhzWrVuHCRMmJNoHW9Q08QE8c9fatWuT7kYKIHtQGyEmOg4SCEq3mwi/JC/L1mqSedAkPb7oOKJ0mtQgeDOHBxPdEZRJD9t3W9LjoWtjn0B9r7CQh8C7RXPuOIEzNqBWfZi+8aLjyuSHPu+mqg+B4v68dytwr8JMsxPAhHHewgXvOUxCAVKhgHgVIN6+k9p/GPRkFhrRWQ52796NTZs2VaV/D9BNiM9LL72Ula7woRFiW3JcCBOBIoNoMjG5/iap+m2gKD2giuQJKD2kf+FID61omJe5CIbBkyV6REimeORnAe8FwZb8aEC3+K2Q/IgwGl4lc9cIO9YWkBwRofetiLhzCtkYY2IeE5Egt6hmx2agmxCf6nRwdk3U+iFZogN4gyxZbCFSfZbDbOJn0YFoFR5ArPJQ15pHeOjxV2jeCknOuD4sErNPIjB3gDZPRMiJ1qIrrZO0YNyUArbkR5Tlmzr3PPLDqyvnU390r13YZ9IVePdwAWoSFKUq/Ta16CBt0WMVZI7NFdQ88aleB+ew6McsQPRmK952IrIjM1vpFic0rfYsOg6dyI+wqo/Kb6lTrvK8DQnpMa14LajqTsOndJgRIDPfmiBotcimLV4+HjPwyA/HiZxLfsLklaLJD+d8L99UIUCyYroAcMQ44AiZOYs3CbsgQFES5XeoJQkQAlRIaP/pQTU7NgPdgPh0Lwdn4vCWdIKs0dQSBWQTvQlBiSPUlZj0FP1SRencuxX8ycqE8AB+VaLUJ1btMc78LFZRwpMQXYjNXX7y1AnzUHTde4qsF5b8PCtfRdf0BQjIj0p5iPLZFcH0xSJJAkTvP8l+yEhmtOpSNTs2A92A+ACeuas2iU96PPztyI6p6uOi0nkn7AmP6eDM+jBxtqfzsvCid9gIH9+AFob0lMD1V2FhkqwvPFyTJc9Z29RsR4qNUteMlzqgrPqw19ZRRnEuOoC1G7yFhxbm30ecoFB/RDB9npMwi74DYEdpSRKEANn0w5V/Tzz4+OOP8eqrr1atfw/QTYjPiSeeiL///e/Ys2dP0l0xBG8g0SE6UZmtWAwvLXnJdmEf3FUwr3QuOo648nhoOm7zVAIKSSAAACAASURBVB6a/AhzurTDCekhoH1XuGqPCuIcQH7El/+na9B+f4RawIfJJgEh+ORnjiiLuQ35IedSdM8z9zZLfljSQ+OIE4CDVASI9+wkoQDZYAf8JCiucZDdbgeCfakd/54XXngBXV1dOPHEE521GTe6BfE55ZRTsHfvXqxZs8ZZm/H5+aRJ1SEghCdKqAp+mrxduiQ8KtVHc2KRmba2Q5HIzrbMBgPWR+VVcBQgldpDCISnqLgOhTdHJ8zVnRAgblfCfEsm5IfXZ41rLVN/eDhI5G+murejeu6jyImzA8D7pSVp7IBXkqL6ylLwsHLlSjQ1NeGYY45JuivW6BbEZ9y4cTjooIOwcuXKpLuiibQRHcBTdWwHPtnbDmvuogt+hhkQOxCdwsPrF50bRTLpThgHTFC8PS/fJGnDhvRwJl+eY2659EW4iYgoLRUTk0vIorIUeXisVB9OpCBRfVhfcyH50UnIKTtP5JorrstaRZJN2iH6oMkSAsQD3T/eOJCm6D/A35+0ECCgQoB0iJBs3EzOv2flypU45ZRT0KNH9dKH6u25AXr06IFTTjkFK1asSLorErgmOy5k3jy1ANGFajbDrsK56DjirAkEaCeDo/OsiKJyfI6r7PGFIT0dMCM0HQjt21MmHvYKjN8/J6nwepb8dACXCSLsAuSHPBcykqFzTKsgfz5KbWxb7y0sRPcbIUB1NtnSXShAcWZAfh9qEhSFeUyGVngvfybJWZNDe3s7/va3v2HKlClJdyUUugXxAYApU6Zg9erVKfLzqWMWEZIY6POQ++24AsmbY1uNmQVLeKI8dx2QZ8Bl9s1LLsdORtxoHdrnwxQ8M0uHWO0xBqOYcHMB6cJPbqLIIO0HIYKmuXeo+1GUXqBMftiXAUUOJy1ovhzwyI8KdVMEP6j6mEew+nra0UotaUEbs7iBS7eM1atXY9++fRnxqRaceuqp2Lt3L1avXu2sTfMbSofoJIVh1GILXUXIZaLATsSv8ABevhOVylWaMHQy6kpDlB359ADepMzO11zSEyYZZJJgVCrW3LWgCXZv50/wt+ORn7cBzNQxe9kSc5b8CNoh6o8q5w89xwrJjwj0OTkaZgQoKrXH9LymkQQB3vhCFpNtosPKlSsxYMCAqvbvAboR8TnmmGPQ1NQUs7lLV9VJCiKyE5VJy5bwiAbIAtTJxFyrPgZJ3o44wUskJ8NOaJAem2NQONWaJKN1gnhD4gOgI9cC0FV9NCLp6HMqTEg5GcBJmvsUgZAf1b3RDnwoIc48YaFuCkWAbO49UwKUJogUlyjMXCZttsOOCLnFihUrqt6/B+hGxIf4+cTj4GxDdGxty6ZtDi0tohBc19hELa4GjwKSyZ7KIzGCQUgnb0o5Gkd07Gwm35Akjud8e9nbnP3rqD0k+3Np0UqCaAlrE1qpbyzh4WZcDgGi+mgTSVdkvBmApknrw1VyAsRD3RTICbTqmc1DbDJPi9ojQ5p9b6hnLyZ88sknWLNmTdWbuYBuRHwAz9y1evVqtLenLylUtBhKLS6gowi5zorcAXvC42IwlCk3zPlgSQ8vt0og/JgdwESTlOpYSD0wpk/CiCPR/g3xKrMkCo3EgwGEVH1uEkQPBlQf+vrZOBTzICI/nOeUJj+qOb3c1TD5iID4fAajQivSHY4uJkFR+PeceuqpztpMCt2K+EyZMgX79u3D3/72N2dtprdu18DSIiM7UWQM3Qo54bFRfQpIrj6ObvK2dnmGXJr8CHOukONXvZmLziE7kZauoYj0cKu866o9KnQgurdSy+SDNEKpPjT5oUwPotxLZfLDu25RkR/J8/vhKuDDlfLmAl2dBj8Bsrmm+dLi6gUsCbCh6HGYuVy0GR4rVqzAwQcfjM985jOR7icOdCviM2bMGBx88MExmLtcmqZM2ySEhyAOdWsr/DlzXE10BfgJT5h2ba4JTXgU+9bJidICjURzuseoezztfBOMsOCpBZyYjhiiFFCQoiRSpE1TZcPg2ZomC/t2SX7WA9BM1Nq50mIfLAHiNqzRTt5i3yrEHQG7C0n73MSJlStXYsqUKVXv3wN0M+JD/HyeeeaZpLviGP0QJDwuoDJpRZEgMIxJSwXdgVGk8ggmXEJ6VNEz29Yr+mCay0jDbELUBi2eY1oGQwCr0PiIISw1YUuiBOeKp/qQMiTSQrSToVXmxAmoY+aRH+VjYkMSOfuuChOY7v2RPgLk0hqxe/fumvHvAboZ8QGA6dOnY/Xq1fjggw+ctZmMuasftQD2kVg2A38BamJio17pVDqO2pnPsCYRq/SIyI8vrwrv+E2TN9JtKUgPAYnkCqX2mOa9SQoSk500os1kQheYJAn52Y5g4Vkh+XF1X5N2VOZUCp0rKwRIi/QQ8NQfG9Ulj/AmsLRkjtaJvKo+M9fTTz+Njo4OnH766ZHuJy50O+JzxhlnoFgsYtmyZRHvKSpzF0124kQBQcLj6mHTITyuIOozyUKrOqbS9ZGl/GfJDzeZHL0fW9IDeKTHNOKLNzDrqj025qY05AQqZVuOK4T/XomfW4D80OcyjOrDXhPD1BFWpi9Az/wF6N0z1ez/wyL58HNXWLp0KcaMGYOjjjoq6a44QbcjPgMHDsRJJ52EpUsTzitihEZqiQIqJ+dCBO2SidqG8IR9O2ZJgmHdIZ30/jshLh/g208Y0sNrrwRRDhmfKSbsgNzhmYtizQdkAg2SxnX8Dqv6lK6DLD/TtNEQ982G/IiOkc6dpXpuOuERYBEJVm0fJjcR3bZpBGpUaodLdTleEuTSCrF//378+c9/xowZDoIKUoJuR3wA4Mwzz8STTz7ptHyFe3OXiOzYRmKZPnBvQy/Lne2gswPhVB4Xg5JFrSGS2E0VCqyVM8WirIAPPALWqUl6CNqhr/ZICIHvdokjt4iOyS2Mz5KJwzG51gbKm8ofzIj86JxvlfrD9ps9d6p9kAeCpFNgYXNPEPITV86xONAJMQlKp5mLuIaceeaZke4nTnRL4jNjxgx8/PHHKYzuIkVKk8ryLCI7LieyHaUlaeQVvzPXzpfNtgQR+SmTHhnZjIL0lMDzLYkLZf+hqCOwLNvl+Tf5VB/bt/Jn+V/zVB9CepSlTHSK35qcB9O8WjL1RwURAeJBdgwq5aca1B4Z0m8KW7p0KQ499FBMnBiX83306JbEZ8yYMTjiiCNSYu5yXZVdBtED9g7CRVHpOOtFQXhsBic6nb6O9A+z+kUBpYd3ziMkPbT/CE1+RDlmwtQBUyZFNAETrk5z8Lfp36OC7QSkofLQ5IdVeoTkh7QpIz8254OXWFRFHlQ+WjL5czKA41WdkoD0LYqo1bghO8/tcJKdHe6tD0uXLsUZZ5yBAw44wGm7SaJbEp9cLocZM2bg8ccfR7FYdNau/g3Xs7TIyI5tJJbJ4G3iVGw76byPdCg8gF39IMVARI/5QvMWfU3Ckp6JEF4LXsSQM/XHNnzZEZySLApSa65KsWDvDQmJXL5JbN4KkB+2XR75CUsCCfnRdOYXqj8qm+8npU+RWmB6HDQBsi1lokJcao8InfCTINk1itbMtXXrVrz66qs15d8DdFPiA3jmrrfffhsvvvhixHuib0xCeHi/xYV3ICY8YfrDbvt+aQnbrgo6g5SsaKJse81olTZo+PS0IzzpoWGggCwnddIMJuoAbBURR7mBrMHZ/01bg4RHWFCUB9mbOe+calyrMvkRtatj9jKFTUmZMNdzIsx8l2TjRrWrP7rQIT/R4fHHH0fv3r0xbVrCLz6O0W2Jz+c//3n0798/JnMXS3hcwFT1IaamMG8zOtvShIdGUuQnr7FvdnuG8Kg271yp6APgOZeGfZPkTRqlNkX5YQI+JmTSNiE9tPIhiU5xmQ3aOWyjauhjtzFFUNdclbX7CF2fH5eKhCwBqWg/RP3RVXtYEALkYiw6tLS4QlRqj4vxT33/RWHmmjp1Kvr06eO03aTRbYlPz549MX36dOfEx/zGU9l9w8LUt8b2AW0Fn/AkhTzMssKSAU/wZsM7LXTiN18bLDZorKOC5E15gmHSRQChB+KZwzl+OHFD9y1U8zkSqj4nwex8rSotnGvNIz8tqNRyE9V6K8PmWquwFUECpOP/JiPPItJD0AZ54V/T+9M1AUozWFNYNGhpacFf//rXmjNzAd2Y+ACeuWv9+vV48803I95TEl77MsLjUvVpLS06iEv1yVvuWzGR0k0Ik73pJJEzPf8a5oGdCPqQCPPIsOHXrq5LO2dJAzgTtNDZm0ZEkwtNflo4v3PJD30+VcqQLcg5Mbk/CcmzxTiYH4+sf4ciOt8fW6Qlq7QZli1bhv379+OMM85IuivO0a2Jz5e//GXU19fjwQcfTLorEpg4ObOVg22h86CaEJ640AG5yiM6LmLa0hjwO6GR4VZVNoBeRweS68E6xirzw8j2wduPST4bWygIZ1lRioFYzSTZuyOerNZu4JMeAh/54R1n1ORHBt65ocmPjtrDgiZALnwNm2BOgJJ2ag4H12auJUuWYNKkSRgyZIjTdtOAbk18+vfvjxkzZuD+++932i7/BrSt3aIDG7Jj+5C3IRzhiWpCycOu4CE76arOy3LoHYPO+dVZh0yAnH2KQqF56k8ZqrdzjUnfyBFYBxxCk4QJTRoBZ0MAZVm5dZ8D2bgRBfnpBPBaaRH9LsIqACsV7av8glT3lum4ZUOAMrz//vt46qmncOGFFybdlUjQrYkPAFx44YV46aWX8PLLLyfdFQuEJSEysANcG9SDVhLIw094NHPzABArDTKnTl47LHi5UkSQRfywJg8D0riWOFOHyQZroHxwTUdhTCAx4N6tFbITWcJHlvxQ51RazgTAQTqRXOPgzu+HvdYi8qNC2DIso8E/JpNnm0UT5CV/0uzUHD8efPBB5HI5nHvuuUl3JRJ0e+Jz+umno6mpKSbVRwZdJ+c2uCMhOg97FITH1WCQF3yvM0CqHGPZNnTDeG1ChHn7E016pYlTlvU34EBLCJApESEqRzzOlNGCPnbF8QgrqIc1+3H2JyI/RLETFcItw9WELbq2r2msw+sLj/yoxhFW3RIRIBuQvkVZ8zA5uDZz3XfffZg+fToGDBjgtN20oNsTn169euHcc8/FAw88gK6uroj3FsYnIW7FZRfS58NDI2+5HUmlb3IuRaSHnQhsSQ+B5iR2xAli/xBVyHRodCZbEsMKnRCSHVkhUS5s0vY3B/dLgyU/rJlSSH7o+yUMQVARGpnpi9cXgmboqz+ysXE0jEvMKEEToEztobF161Y0NzfXrJkLyIgPAM/c9dZbb+Evf/lL0l3hgPg9yB7OMA+YyEE6atj2eWhpUQ1WvN/ZCURFfjqgVnrIcYQlPfQ+NZPV0WHQSjiqAs+WxCiToBgGea5vEY8USIiOLgKqT9iSGRrJK2W+WQHyw+tLFKHuNAqW2zXDzUubyP/HNE8XjUYA/ey6kxK4Vnvuv/9+NDY21lRRUhYZ8QEwadIk5PN53HfffU7btXdyTiIcWER4oox0MJ2U2IKFJuTH1ExBkrTpHL/LSX8cpNeeF+pMyI9TtcfgfJVVE5ZwxGEeM9mfjc9R1DXCStim0bcy+ZH1x5T8mF6fAuc7nedwA8RRjqpxjv59ONQO0Oy+ZSDH3w9uCVB1qj3FYhH33XcfzjnnHDQ0xFE/MhlkxAde7a4LLrgADz/8MPbs2ZNgTzpgT3ZsH7R2xKPwhAFReWzQAfkkznsTZVUe2eDJS/5mC9Zvh7kXZMntWgDxPeBI7bGGiBBJlns3+RUln7LkmFSx5i6ivAgTQ9pWqRapPqXrrCx7AqBOh5DqEgPb81eg/jZVXk3JuWg8JMc4zLA9FVwToOpCc3MzXn/99Zo2cwEZ8SnjggsuwEcffYRly5YlsHfdt0qXb54milKSqo+qJo+sb7q1gWjyIzJt8fbjivAA4tBkzWtU9hOJUGUROv3WCKRpAFyBJj+cZ1BEfmgXv7opkvbJtVeRn7D3SAH29z+t/oRRtcOGvsvOQT/Ii0jbtusWUTg1DxkyBFOmTHHabtqQEZ8SRo8ejRNOOCFGc1cH3MvoOg9cmjLqEvD6TVdhtrHh04RHlVAN8GYVlT8PvR+XpEcFxTXjRgYRAmSj9pA6SmHuz5SHspdROkZjwmOr+hAYPIPa7jHsc2LrE2OCguQ31b2j8nvSNYHlYR7soDumNMCeAFUXOjo68Mc//hHnnXceDjjggKS7Eyky4kPhwgsvxH//93+jtTXN0Uy2E1FYwhNnVlObysukfyKVR0V+VkE/MaFr0qNIRHfQZG+JXI0QwZYcpBGkhhZD6kT+UbK0AVYV0zsArJGvQqs+ItITUH1E9y5LflyRHno8KCh+54EOCrAJDOCNZXlUCFCY8Yq3rS4Bql6156mnnsIHH3xQ82YuICM+Ppx//vkoFotYvHix03b5N2hUDwjbbhoVHh464Vd5eL+roHoLF5EfNseLDK9p9kUXoyEdpNloHtYcI02CpxFJFIDiHO5kFuNw8CgQR1kNGrbmRPo6K3xdPlylVnrK5EfVF0J+oiA9BAXF7yqw91GYMSuv+D3MeahdBej222/Hsccei2OPPTbprkSOjPhQGDhwIL72ta/h97//fQw5fcJAZ2CJgvBEqfocqtG+bMAaBz2bAEt+eCYZF8ncdED7zBialeJWX4TKRwdnqTGUj51HeHRUH9PzUlpfWRMOAKZotpk32L8MsuMoQC/kXfTsEPVH5QCtYwIbBjvHZ93rVFsEqFAo4M9//jMuu+wy5HK5pLsTOTLiw2D27NnYunUrnnnmGaftxqf6tCPaRIeuJ7ZDS4sueOeMnpR1yY+qqjS7H14CtyiuH3V+VVl7t62CeBKwUXtcIWryE5ZsCfyeeOaucq4k22st6xtvgmfWl5Gfcpd0s5AfrVhPBd3z/I7ktzA5dwCz0HcgSIBcP7M0AapeM9cdd9yBxsZGnH/++U7bTSsy4sNg8uTJGDt2LH73u98l3RUF2AGiWkxaNHiEx2QSoys601CRH12HXzKQybLW2g52sgipDjXp8UX/uLr2YR12bcEjMh1Bs9pOsm4MYJNDylIJcFUfXUJGyI8hgQvcdrp152zJj+l5l5EfHbhMEwHoqT9h7q3qVYD27t2LhQsX4t///d/Rp0+fpLsTCzLiwyCXy2H27NlYunQp3nrrrRj2GPYtQTTpRfn24WLykak8OiYvVWVqEfkhpEf3GFSp+kl/TKARFm5VoYTcC7ZqjyA/jsjMJUyYmHTeIFt0AmvXG2TD5sFGgVKYdmjVR+pexJIfUT9MyY/t8/4O/ATIRu2hyY+p2sNiIOzzgalAji16AuRa7Xn00Ufxz3/+E7NmzXLabpqRER8OLrjgAvTp0wd33HGH03bd3rCdqD6FB/CqJDeFbGM07I6dnZBVA3rBoG1d8jNcvS4dscMjP8pEdy58beLKvGwIaZSVCci9YHCcUtVH4aQeFp0rNS8FIT+qvhwNPb8fF8f0DsKZuHTUH9V4QLfPS4ga5jh5x1Y9CtDvfvc7/Mu//AtGjRqVdFdiQ0Z8OGhsbMRFF12EO++8E/v27Uu6OwxMJqM0qT4s4bEZCNlqzarBjmYNIhVCdBwFRds86EbXSNblJaczUn9o9cCEACmyQhvVBUsCNmY6F8+HKzOjKjquE+o8U4B3vU/RWC9uQrsj5PadsK8VJkJU6g+NutLiBq7VnpdeegmrVq3C7NmznbabdmTER4BZs2Zhx44d+NOf/uS03XBOzil7+9aeVEUKjwn5sc0a3Aa16YU9joLlvoBIiWkb9MoaBBBCAWJVDh8JStv96AjSFAEEPMLjSo2iYXI/6V5jur28g/Z0sQN8AmQSzVlA8Pk0UXtYDEU4BVr32rgjPy7x+9//HkOGDMGMGTOS7kqsyIiPAGPGjMEpp5ySEifnMCaHpCensGatDshJj2rQ2wCzciAFjXV1YJJFt7SutBQBSr4eMhKjCgOOItycV2vLJrFfFWAngINOQHQmZlb14T27uiVVRM7OvDbzGu25BE1+bFNYFDT3pRslJsshZgt23+HUH9dqz65du3Dffffh0ksvRc+ePZ22nXZkxEeC2bNn47nnnsPGjRudtquv+rjysUjC5NVYWsLk5gH0Ci6KJiLW7KNCFKQA0DsG02sUhsSwYeAWRMV54kQJWsBRmniLKTT7aVTDK6zqQ1e6F4ElP6L7gCU/sjbzGu25xA6Ej/wqwH0WdRPyYzuuujV/2eLee+/Fnj178O1vfzvprsSOjPhI8NWvfhWDBg3Cb3/725j3nEKnUinYgbJR8TsL3rEOR4Uw6Lxh0+vQRRB1+1EoLVGRHxWmyS+7MJ8LITCmVa/p7V2mQrAhPQIik6Rf0bZVYsKjSjUQGjrZsAn5Ud2vhPzo3Id5xOP3ogtdU/jbit9FEN3zLtQfnXFEnwC5Vnu6urqwYMECnHXWWRgyZIjTtqsBGfGRoFevXrj88svxhz/8Ae+9957TtsWqT1SEJw4iRVQeG9D94ykkuuRHx9zDoqCxji2OhlqdYd7MY+W8tELRDh8JkkYx1SKY47dGHFXsOwA8obnuSQZtxgVyk79fWkS/i8D2VUR+RNAJje9XWnhw+ZDGrwAtXboUmzZtwhVXXBHrftOCjPgocNlll6FXr1741a9+lXRXUg5V6KbuoCozC6kGq03QG5DovhQ01rEFmy+F16bAF4PmwMrSBVHUy2o3NPGkATZEbQ2syA5X9XHpQyUz4Zjsg4QDqlQqus28Qfs24D2j7yt+18HbqBAglyRORoB4sN03n/y4VnuKxSJuuukmfOELX8CkSZOctl0tyIiPAgceeCBmzZqF3/3ud9i50+0s4PqGViMKKYEeFMLk6shDL7uqaIIynfw7oHaQDDN4ipLEGbbZWf5fcmCzJ2tFPaUVbIZoCbQj6HhtuVB9WPLD24/I2ZmXA8GErOVVnbOE7F7mKT88hM2/ZZsIUXecSzdWrFiB5uZmzJkzJ+muJIaM+GjgiiuuwL59+1IS4ZUmmLwFEagGWVN/HkAvGoaFrmOlDflRZcYlk42qxhJQmdhExxRG7ZE44kr9WNo5SxrBkhyHKkAbgLrJbtuUQrYflvzIkj7R11XV97zid1PoPJfvA2iV/K7rL8hmjdaF6l5ugNyc7zYRYhQvxzfddBOOP/54fOlLX3LedrUgIz4aGDRoEC6++GL86le/wiefsNW9w6E6VR+Xtu885zuTiVQ08cv6YTogRjG5TdNol53Q0uz0TghQRERDa9/s/iOAUTJJV6qPzrGQe0W3c7rnx5Wzs+l9KyM/JmCfdVckPYw/Y3JYu3Ytli9fjmuuuaZbVGEXISM+mrjqqqvQ0tKCRYsWJd0VBwgzeeqkYdc1eeUl6+jI0TrZblnYhtDqTqZ5wX7DtkuDEKAofHtcg6e6sAtPQWrnFyn15dJxpTZJHOI/XFUhO8a108KCXOfXNNd/UnO94zXXI/dlUpFeLPmxTY1B1J+CYnubRIiNit91EZ/aM3z4cJxzzjnO264mZMRHE0ceeSS+8Y1v4NZbb0VHh9u3yfhVHxvQtWfCHv8whAsX3Qr9N2F6QAmbNwSKfeaZ/YoGYp6JS9d3wxVIbSmXColtSH0aoXleVEkntXI4sTBV9kzWJ4q1qsQHe9xhyE+YFy1CfsLmAwsL2f4bkYa8PCq8+uqr+NOf/oSrr74aBxxwQNLdSRQZ8THANddcgzfffBNLlixJuisOYDJQ8FQe24GIdmC2cTLkOXyq0Ak3pEe2z7xk3zRkfj0mBMR14jZ6/5I+WJXMqBYIjl8ZVceDrVlStI1I9WHXl10f1kwvIj+i629DflwQklaEl9vIMYnKZrhQDm0Lk8aj9vzsZz/D4MGDMXPmTOdtVxsy4mOAsWPH4owzzsDNN9+Mrq4up20no/qoBiXVg2xKfnSitljQA5Joslf1I2yBRNU+84p1TQf/J+BGgbFE3ZQEzTthYZs5OYRa5VN9eIRHN3O36j6hyY9sfR75EfkmsuRHdd+ZkB/XKozoZrSJJqXHhDC1vnjbp68q+5tvvon77rsP//Ef/4H6+vqku5M4MuJjiDlz5mDTpk1YunSp87bTZfJy9fB2wiM8ItKjG8WlUjhEg1MUpEe1Tx50FQDaxCVqP6zaY+h06yNBCRKytKKccsBW4THdTmd9Qn4+gZj0EBDyo3ttdczUUZmeWPITJoWGSP0xgWj8Ii+NqhpY8ag9v/jFL9CvXz9ceumlztuuRmTExxCTJk3CKaecgvnz5ztXfZIB++CZyrUu3hBV5KcA88SEQLSkh2AoZ788TC4tphJKnNFR0PBZAcROyjaIoqJ5BOCZu4w4C0/1sSEHuo7OALDSYF0dMkwfsIz8RO1vQ1i4q8ivXZLfXDx7IvITT4TmO++8gzvuuAPf+9730Ldv31j2mXZkxMcCP/3pT/Hiiy/iwQcfdN52sqqPrcqj8gnQGTxE5KdA/W1CfuIgPXnOfnUgsiHJHJoJuYjCt0cTVr4uUcIlAdMEXU7MB52cTGwjtihorKN7Huh7UUZCef3lkZ80pVswuRd2IUiATE1cLOhzoVf9PIrxf/78+ejTp0+3LU/BQ0Z8LDB58mSceeaZuO6667Bv376ku+MAnQhv2lJFgdiQnwJnHRvlJwrkDfY7WfC9jQNN2ElzuIM2eLANrZeEudMmNnaJ5BrLckJ1OjhlecQTfUSfm2bJurr3n6zPNPmJk/SQYxSpNbb3h0z9MQHvXPREhQDFc642b96MRYsWYe7cuejXzybhbG0iIz6WuPHGG7Ft2zYsXLjQedvxqj7EtOXqQRwKsXlLl/wUIH+rVaW9f1+xTljkJb+xE4+I9BCQmVwnfF3XuVUXnZBIGN0cvPOiukYi1ScKJaog2A8PPPIjIj2s6qNzXwwE0KSxniuwx2mq1KiwC2ozWpgoML76E8W4f91112Ho0KGYNWuW87arGRnxscQxe8P8mwAAIABJREFUxxyDiy66CDfccAN2796ddHcswPPlCTv56Qx+qkHJtMoyDbbWT5KTeQfUpIdgFeJxjlVFF02pEQ4UJmOyqxMgMr2pypmYoMDsTxcqpYeQH93zEKfDu2hfPFOVLchxi86TiYlLb/soSM8LL7yARx55BDfccAN69+7tvP1qRkZ8QmD+/Pn48MMPcdtttzlvO1rVJ4pwSyJ567wJiQYvk8rK7OAiKnDoehY3CefViahhYZrHhd7OIWtJuxik5YRtg7A+VKcgXiKgm8izGR7B1jVv6SZepPd9qOY2ttA5zrDkh73ZXed0iK+u3Zw5czBmzBhceOGFse2zWpARnxD49Kc/jdmzZ+OWW25BS0uL8/bdk586qDOMms5yAxF0crR5uFmlx4T8qKo6uzTjAXqmCzo/ioz88HKuhGUbEbCVzuVIPxuKAyJzl4k5y4XqY2o+M1mXPL8q5YzXZlTkx6b/LkHIT1i1J4goXnKffvppPPPMM7jxxhu7fZZmHjLiExLXXnstisUibrzxxqS7ogBNeFylf5eFtJokBROZt2yUn7DricBTekT942XEtVV/OmEWwszbvhPqN3iTqCRe+2lBnP2JKZKMu18aBYP1VYka2edWRH5kxxy18qMDmzpuqvtG5fdjbuKKAl1dXbjmmmswadIknHnmmbHss9qQEZ+QOOSQQ3DllVdiwYIFePPNN5237+ZtIK11ZDqg9umRDbCtMMvlEcWEaDrp0eQnzvIPUU7QEYXY0+HzrNjUKfwyYnRAnVlbh0TmLfdtqvLw1ueRHxlRYMmPTh9cOjuHuW91yYbuvWNbGJe/TRRqzyOPPIL169fj5ptv7tYV2GXIiI8D/OAHP0C/fv0wb968pLvCQGbaslV9+lGLTa0tGu9I9kOD11fb5GU2k6PKr4eeYFTFHwEz9aeguZ4uklIpROCRF2pJRFBiiZyrcxaGoKn2XTBcn4bORE7Ij4kK64L8uLhPXSgtbD/YNtNh4uro6MDcuXPxla98BZ///Oedt18ryIiPA/Tt2xc//OEPsXjxYmzY4L5Ctd0DoqPymJIfmzwQokGHLhpqSn7CZmyNaiY1GaSbER8BEfmU6Cb+0wm1TxtsqqLTCEt2aNVHRHbyBv3QQYHaRgUyTumSgnbo9Zc9zjDkx+WzITtO1Xgg6kc79BSg+Bya77jjDrz++utV4HqRLHLFYrGYdCdqAfv27cOxxx6Lgw8+GM8991wkEmMuN19jLVOzlk5G0TrISY9OlBhZR1YlXafvrqtmyvZJ+zDp9I3OgdIoWU+UVI69FgWNfepA5UzLM83QfVERnzCmrjC+Ry77xU5+Og7IKpNWG9TmzILg+zCTvm7kYSf0Qv91CgWT9kQwfVmJ8oWAHq9sSQ8L2TgaT/j6P//5T4wYMQJf+9rXcOeddzpvv5aQKT6O0KtXL/z2t7/FX//6V9x///0J9cLGl0fnwVYpPbqFRmWkB1APQlGUChftk3XcVpko2MRvNn1Ni/kJiLUMRCJQ+QbZOpQDZiHQeebfcZxv+phVGbfZZ1tEVFXPronyE/Xxk2NyqfyKrls8pAcArrnmGuRyOdx0002RtF9LyIiPQ0ydOhVf//rXcdVVV2HXLlfJtCoQPzA6YeoyiAaaxtKim3FZhh2a7YgGoyhIj2qfuuuKahzxJkBZCQGgMoAWDPokg6uEeSwRqiZCFIcTtK+ehiVcnU/ZC4bJsYueaZb86LYpU0EJ4rqnwuQb090mPhPX6tWrsWjRIvz0pz/FwQcfHNt+qxWZqcsx3nrrLYwaNQrf+c538Mtf/jKSffBNXi4it2i5ljdI6ZjFeGYvumCoXrE+//FESXp4+5SF6bPrAnoVxhuhJj00CqVP3fMlQh7ye0MnAkllUhKpIzr3pCtTF2+SCqPaAGrS2Am97Nwycxfpd0GnQ4agTV4qcsKavHQm7a0a7RLQ10f0PMdJpOl+88YsF30Jnpso1J79+/dj4kQvqKK5uTnL26OBTPFxjMMOOww/+tGP8Otf/xovv/xyjHt2+RYrejOzGQzYKum6bajSxkeBTuiRHrJuJ/RID+Adh6mTKhBOWcmXPsMoHbakh92vaBEpSTrLKYhWeeIdmyvVKE7FTKevxORlEq49THM99jh540tSpAcIHm/1kB7Ac2hev349FixYkJEeTWTEJwJcccUVOProo3H55ZcjCkFN/AC5GIxVcrRqUKAHEZb06LZBECfpAbzEayYT0mjoTxIkikbVfkHwfRpCqrsrXJ4z0TXMO2ibhW66CAKTiFRy3+cV68nM6Kp1ooDofLg0S8X3bH3wwQeYO3cuvvnNb+Lkk0+Obb/Vjoz4RIBevXrhN7/5DZ577jksWbIkkn24f3sgRUtdDADtEJMeAtVgt0tjnShhmgfF9LyFITGu1IKT4MY3Ja1w4d+kmsR0klBORnI+UarnkID0TRUJx7vX84o2RWhENHUDRdDJrBx2/OPvIyq1Z86cOSgWi7j55psjab9WkRGfiHDaaafhnHPOwZVXXhmJo7MYNm8b7OBjUm6ChUmVZFE7u5h14pgweGn2ZfvlhQGLzpvsTZreR0GynmhbUR/zhm0BfhJUi0QoCZgkq8xH2A8RTJ4v2biQt2gzTtUxbACDPaIiPWvWrMFdd92Fn/zkJzjkkEMi2UetIiM+EeKXv/wlPvroI8yfr5N/xxxuTF6iNy4b8mND8Nh2RG1ESX5ktYV4g7gs90nc6g+9vetz9CTiLwmRNhQstvkEQcKjk9E7KpianHmqj+49berDB9glRo0aNvd6fM/H/v37cdlll+HYY4/Fd7/73dj2WyvIiE+EOPzwwzF37lzcdttteOmllyLZR7i3ibAyMz3IhVG1SDuqNqIgP7oFFU33TSYKXb8JVY4jXUStkHV3IsQDMXfZFKKl0QH95IOmMA0yoMmPLukxCQxgESX5sb1Pw6tEUak9t99+O9atW5c5NFsiC2ePGHv37sWECRPQs2dPrFmzBj17hg1P5sMsxL2nxjoEOuTIlWOgbjuuziFJqmaaCkA3kotEyui074r4AEGTCe98hQ3DBtyEYNtklSYIkxnZBHnJb4RA6Co6vJQGLAlxeS/QGMjZlwyd0De/0c+urPCwiky4dgtwRc5Vz3B8UVzbtm3D2LFjccEFF+D222+PZB+1jkzxiRj19fW4++67sWHDhgQc0HgPvSlpUJGRNsF+TPdhQp5cKxq66sXw0mJK9FTtRzXREaQ54aAq3J31O0rCB6lA/S1K4GiSo4lti0VUqo/JfRYm0aEozF2nzX5w5/DsUpGUtRWf8tnV1YVvfetbGDBgAG699dbY9ltryIhPDBg/fjyuueYa/PjHP46kiCmg4+/TE3zSo/PQiiZ6evKxffht1aKwkzgvhb7pYC/rO68UQBrMQ2SyDWuWqRbkHbThmjjGTULp++59zfUJCop1Rc8AS35sEh2GQRTPGa/N+E1cK1aswF133YV+/dLoG1UdyExdMSFZk5fOG5SOOYZuR/TGbWI2cmEiszmPqrpBvGOQZRjmnV9VDSSyH9dqT15jHZFZ5lPU3zph2gWNdVTIK35XmeSi6Kdtnh2VuYs8M3H5fYkmf5Ffm4ws5Jl/6z67byvapcGed9vxIY6XizrhfjITV/qRKT4xIVmTl84AYjJYuDAzuPILMn1D1CmWyCozqrIKrPqjQ3rIftIEOhopraaxKCAyX5mCZ+5KKj2A7N7iqT6qe7FA/W3y7Oo4PIvOu43JK65nKl7Sk5m43CIjPjEiWZOXC/LTDvUAHsZ0ZguTycrUudME7dAnPYB+crkkwSsVkSboOGnTsD2egsE+VGRHxznepHwKu53OfUuTH937vACzsHZybmU+S6rzb0J+0vYi4Q6ZicstMlNXzEjW5AWoBxKRqYoe8HT6rNNOFJD1jU6Tr3ve86VP3QGYhAHrtO+a+NATjGz/OhFIKmfdAuc7m3s5r/jdJvqMnUwLup1RIK/4vQP6EX+iFx928ja5R2wmfh0FlKADevW5RGSGNd2ZZkeXIXnSk5m4qgeZ4hMz4jJ52Ss/vAHEpoifTjtRQKc2EFlPdRx56m/datW67Uet9ogUDRekR2efOosJPpEsSSpS7H51lVyWIImUGh3VJ2qnefr4ZKHqZF3dNk0ge/GoXdKTmbiiQUZ8EkAcJq9woAcS0YRvSn7iID0EbN9khVdNSj4kkZVZFypzQhrNVDrkyEX0WT7k9gQFuCdZYUiLC8LTqvidd4wi8qM6H+QetT1vPPJTu6QHyExcUSEzdSUEYvLK5XJYs2YNGhqiKdZnb/IC9AYoHfNGUhMu6Zuq4jy9LkFesT57/rZy1+LvJwq1Rzf3C6s0/P/tnXlcVFX/xz8sCohrKIqCEqbIo5mK5J7gnihSpo8hhuaCC+X68EiISoqWiZkoiok7irvihpGCWolrlksqpZiEApI/FEVkmd8fPEPDMMu9d+42M9/36zUvdbxzz7l37pz7vt/zPedoOjdMIj6ZDMvTh6ue/+crQpXJYBsmMDnP+rq7lPLMNR+M75u9pu4uJr9Z1W4vpr9xJgsYM9mHEmnFR0jpuX37Njw9PamLSwAo4iMRNjY2SEhIQEZGBmbNmiVYOdy7vPgKWRdBusapBMykR7mt8lhcGWyvGv1hKj3KcqRC0w2ZFiblH21RXPWIoa4139RpDOG6tNSjPmwHC7CRHoBb0rYqyocO6aM9QlFUVISRI0fC2dkZ0dHRUlfH5CDxkZD27dtj5cqVWLduHfbs2SNYOezlh+3NWdv2Uj+Z1QX7mwXbY2fb9ZUL/te64nOm32cQL2fGVcB9y4UicOsiVUWMbkql/LApJ5PFturHb4j8lIK/ZWu4I2S0Z/bs2bh9+zZ2796N2rVrC1aOuULiIzGTJk3CyJEjMWHCBNy9e1ewcpjJjyENLJPJx8SUH/X+cKZlO4PdjSoT0udp8ClRmqIVch/SLiZMJxW8AmbXkK6oj/q5ZrqgLle4TPWg73zo+i1xHa6vRJj0ACYIKT179+7F2rVrsXLlSrRvz3SUIMEGyvGRAQUFBejUqRNee+01/Pjjj6hZs6ZgZWnP+WG7UKc2akB/g89XWdrQlQSoq2xNQ3X1Na6ZLMpgslwAl3Oj6QaiaT9MGlEmyfZ8zTbtymAbNjk+um7cfM6QrS3Cpi6fTLuy1HN9dB0Hk2vIEPR1DWsTbE3nhOnDA9OcH21lizlwQljpuXfvHjp06IABAwZg9+7dsLCwEKwsc4YiPjKgXr16SExMxC+//IKwsDCpq2MgXIe6i4W26Ii2+Ul0PbFmsiyDCWw/q+2pWb07TY5PjkxGdelapFT5kjIKxUe3pVKQmByH0FEfXXlehixcqms7JqOVdJUtXuRHSOl59eoVRo0aBQcHB3z77bckPQJC4iMTvLy88OWXX2LFihU4cuSIYOXoX8zUENjcAIRK1GQ65JPLrMyqZDIsg83ikNo+ayil4CfXxJjhMw/qLzD7fpiM2lJ+J3LvPtR3rKoRNTbSw0fZQIX8CBtJFlJ6AOCzzz7DlStXkJiYiPr16wtalrlD4iMjZsyYgSFDhiAoKAhZWfomCuOOMPJjSISDL+zAPk+hFMxmowW4i0Mp2M2Qq/5ZviWxCJpFSI5zSjFBjtEsfRgioWJGfdhce5ngLj3aHljEXtVdGo4ePYro6Gh8+eWXePttJt27hCFQjo/MyM/PR4cOHeDq6orU1FRYWwv3FMNfvg8fN2VDj1NTuJvJyA9lJIDNKBGllDL9jL5J4tiiPFdMk0PZ5pro+y74ypfhY04cJWLmJgHczr0+KWAyLYLyt8b3NaUJNl1IquKh79zoOg9PVf7OVXr4jyQLGe3JyspChw4d0LVrVxw+fJi6uESAIj4yw8HBATt37sS5c+cwf/58QcviJ/LDZ3cM3wg1HwmX/fOJ8ilcqPJLNbwIbvwGfroY1b8HrhFEJih/C0/1baiyvSq6kpX1nQdl5MeQSA+/D4tCSk9JSQkCAgJga2uLzZs3k/SIBImPDOnZsyeioqKwdOlSQef3AYTO+WEL1zJ1PZnqkgNNEQd9MqHeBalPmIR6Mld2d4g1xFx54+UzasIEPrrfhJA3fSORuHwvrTS8J7Z4Cjl5IVP5Yxpp0lU2P/IjdF7PjBkzcO7cOSQmJqJhw4aClkX8A4mPTAkNDcWoUaMQFBSEn3/+WdCyuMmPUA0y230yXXpDvZFksrYVG+Qwv42mmy3XIdXGgvKGqilKJYU0iLmGF59RH2311hb10XeMqmLIJtql3K++YfVMzrFh8iO09MTFxSE2NhZr1qxBz549BS2LqArl+MiYFy9e4J133kFubi4uXryIxo0NnepdN8xzfsS4mTBptLgMY1Xm5TAd5aOax8Mm4bwGhM3DYJrg2grMzhNT8eFznTE+c5Skqj9TydEUzVFHKQiZLOpg6DXGtP6qycdsxI7pKEtt+9Q0vJ6tWLJvr4SWnjNnzqBv374IDg7G6tWrBS2LqA5FfGRMrVq1cPDgQbx69QrDhw/Hq1evBC2PWeRHrCdofU+8XOfuYPtErtye7Sg7NuuEsYXtqB5to7gI9qhGkvjq9jPkezEk6iNGVJNJnpCufar/hrhE09hFfoSWnszMTAwfPhy9evXC119/LWhZhGZIfGSOs7MzDhw4gIsXL2LatGkQOkCnW36MKe9HG43/92J7o2F7/HJe9JNEiDlCdZtJuYaXUALDBSb1qAPDuxCZyY/Q0lNYWIhhw4ahTp062L17N2rUkH7NMXOExMcI6NatG9atW4cNGzZgzZo1gpcn9I+fPeo3HL5mamV641F9sjemEU5MuleAfybQM/e1uITMD1IOVWdyzbnq+X/174hpZNHQ75arpHHNEzK0XHV0y4/Q7V55eTnGjh2Lu3fvIikpiZKZJUToRZMInhg3bhyuXbuGGTNmwMPDA3379pW6SiJTiorL1RDp0ZZTUsRyv8oborafj5CRntdQ9YbMx09Y0/wxmm5KYswdIxZ8SU0u9Hc9Kssy9AbO1wLCUvAU3PKElOfMDkJGJ8V42Fu0aBH27duHgwcPol27doKXR2iHkpuNiNLSUvj6+uLixYu4ePEiWrZsKWh52pOdpcKQaemZJtKqCxCTPA71OgktPkzrwTTiw2TiPED3sht8La6qCT4X+xRikU9N4qNNrFxZ7Dfzf38ylQQ+EoGZYMjDB5vPahIdQ+RH83cihvTs27cPH3zwARYvXozw8HDByyN0Q+JjZDx58gRdunRBjRo18NNPP6FevXqClicv+VFtNNneaNmMiFOWwzZ51RrC5/UwTWZlKgsAP+LDFmuwFx8+Zj4G+JcfpfgwiSK5MtxnEdgn1Ktee0JHebjIj/L8MOma0/Vdc102pjpiSM/Vq1fRo0cPDBkyBImJiTRJoQygHB8jo0GDBjh06BCys7Ph7++Ply9fClqefPJ91BtaNl0VdcGusTRkTS4hbzhsRvDIPYFZea6YvOR6DEr4zAdSPVama8gp4SMJmClcBgcw3a++fbOVLumk5969e3j33XfRpk0bbNy4kaRHJpD4GCEeHh5ISkpCeno6xowZg7KyMkHLk15+tDV0TBpT1bwCNo11DtjfRJRJnHJMENYmQkyjJIR22OQ+ZWp4jw9BlevSIprqpS0qyub4mY6Gkk568vLyMHDgQNjb2+PYsWOwt7cXvEyCGSQ+RkqvXr2QmJiI/fv349NPPxVlmLs0AqTv6Y5tg8/lBiPmsGGxUJ/xmBAfJteivqiP+vfHdMJAPtBXd13Xlbr8MP1Nqh6vLvnRfl2L0Y4VFhbC19cXBQUFOHHihOCTzxLsIPExYoYNG4Z169YhNjYWUVFRopQpffRHG5oaOV03AV0NrbbZfXVJjK55TvgQIDbdGK4c9q9vKLcQScHmSAkqomx8L1oqJ5jWTSk/bKRHHU3yo71sMdov5YSzt27dQnJysuCDUAj20HB2I2fixInIyclBREQEGjdujIkTJwpepkKxQKSkZ2tU3CjYhrWtwezJV3WorBImi0+CRZ34+qym/fCxL13I9cZqCEId09/QnoPFZ9SPafeuUBMOqqM+FQSb81uCivPGJGeH6X6llZ7y8nJ8/PHHSEtLw/Hjx9GxY0fByyTYQ+JjAoSHh+PRo0eYPHkyGjVqBH9/f8HLFF5+VC9NtsLA9ubGdh4foGqd2N5k2B6PrlEwYkmQsSAXWeNTdhqDv+UxhED5+2ErPUzRt98a/9uftNKjUCgwZ84c7NixA4mJiejTp4/gZRLcoOHsJkJZWRlGjRqFw4cPIyUlBb169RKlXGHkR5ePM7m5q4oCWxng8qTMx4gjXfVks+aXMieEichlstivEJMXMh2l5spwu0wWZQtxPFy+J12oygEX8REr6lMKdg8OmqTHkAEMgK7foFjd81999RVCQ0MRExODkJAQUcokuEE5PiaClZUVtm/fju7du8PPzw/Xrl0TpVz+GxV9QUh9T4qGLGoolfQA/C8XIffh7OaOtjl6tF0HzVjsW8ykejazUuuql/rn2eQwSS89W7duRWhoKMLDw0l6jACK+JgYBQUF8Pb2Rm5uLs6ePQs3NzdRyuUv8sOm91U9SqLviVtXVIXr07GQUlED7Fd4Zzr3ix2kj5CYWsQH4B71YSIq+qI+mvYh1PWpTUq0RW6YihjbLjPppefIkSPw9/fH2LFj8e2339JcPUYARXxMjHr16uH48eOwt7eHj48PMjMzRSmXn0aGbcoZ26dabdvLUXqAf/IWhMhbMcWFSV2lrgBLMsHu3GuL+ogd4WF7PbKpG9PflO5IpljSk5ycjOHDh8PPzw/r1q0j6TESKOJjomRlZaF3795QKBRIS0tD8+bNRSmXe+TH0Dx7NrMaA1WjP3IVH0Dz6DRt54rNTL9Ml0OoAeGiI3xHfADmUR+pIj7q0sCmCwv4J+ojhEzog6nwKKM+XGVMX16e9MIDACkpKRg6dCgGDBiAvXv3ombNmqKVTRgGRXxMFGdnZ6SmpkKhUKBPnz7IymK77g83uDU8hkoPl5WblY2ysUkPoH2uHSEwhUgQUD26xbYLkSmaZiXm8/t6DeJ/J1wmCTWkjro+Kw/pOXXqFPz8/NC3b1/s2bOHpMfIIPExYZo3b47U1FSUlJSgT58+yM7OFqVcdrM88zmjAttEXt1DYI0DMSVIjqgncWtb70tsmH4vTEdrGZKkbshq6lJdV+zylcSUntOnT2PIkCHo3bs39u3bBxsbG9HKJviBxMfEcXV1RWpqKl6+fAlvb2/RIj+AWI2Rpkad6Q1C+XTO5YlWbjSD9Dd7trBdpFTTSy6o1pfvBUulOE5DZZpvYZKH9Jw6dQrvvvsuevTogQMHDsDW1la0sgn+IPExA9zc3HD69Gm8evUKvXv3xp9//ila2cI2SrqeZLncMOQUOeFjzSU5i5BQXU1iIeS51Xftcln3iWnUh8/fgKH7UZVezYid0+Pr64t33nkHSUlJsLMzJJJGSAklN5sRmZmZ6NOnDxQKBVJTU+Hq6ipa2fqTnrl0eTFteDRtp22FaFU01UnMp2828sMmSVbZvcL0nDM5V2zhe7I/JWwimmyOi63csBVXtsn5+pZW0YSua1co4TekK1vamZhVSU5Ohr+/P/r27Yt9+/ZRpMfIoYiPGeHq6oq0tDRYWVmhd+/e+OOPP0QrW3/eD9uGl83TlvpTI9MbnlTRH7v/vYSO1uhbmBQQRnrkjlg5Qobsn6+oj9BRTq77lo/0HD16FMOGDcOAAQOwf/9+kh4TgMTHzGjevDlOnz4NGxsb9OzZE7/88ouo5euXHyYNZQ0W26rCJVqjWo6UOSW6bpJsh0TrQl2ETD2cL3YitJy6HsUSezbl6P5diy09CQkJ8Pf3x+DBg7F3715KZDYRSHzMkGbNmuHs2bNo2rQp3nnnHZw+fVrU8vU3Xvw1lNX5G9xuOnLK/xHz5vkUmuVATjdwXWiKaokhdarTJDA5V2y6rgwd4SXFtcykPN3biC09K1asQGBgIAIDA7F7924asm5CkPiYKY0bN0ZaWhq8vLwwcOBA7N+/X9TyucuPtsnNuMwmy2USOKFvGlxWiWdzIxRilW+mo7MMlSUm+20M+Qzx51MMNY3wYptHJIdzogn9UR4xpae8vByhoaGYPXs2wsLCsHHjRtSowXaxY0LOUHKzmVNcXIygoCDs3r0bsbGxmDx5suh1YJ74zLTx0ZZQqUsOdO1b1+f4nIcIYC8+unI9NO2LrfgItcK3UInbALtjFHIFcz6+S31Cq6/+2oRCqm5b9d+LvKI8JSUlGD9+PLZt24aVK1di+vTpopZPiAPfrTZhZNjY2GDHjh1wdHTElClT8OjRIyxYsEDUNWcUigV65KcU7C5VTdvra+iVT+Zsn+yUDbccf0rqx2zq+TrGTg4q5IeNlNSFZvmRY2QHqPrblJf0PH/+HCNGjMD333+PnTt3YtSoUaKWT4iHHFtrQmQsLS3xzTffwMnJCZ999hkePXqENWvWwMrKSrQ6KBs57QJUCnZSwlVIStTKYXoTUm3Euf6shBYTTUsJCLFiPVFBEdh/p4ZGYtispyVV1EdewgMAjx8/hq+vL27evImjR4+if//+oteBEA/K8RGQkydPIigoCK1bt4a9vT2cnZ3h7++PK1euVG5TXl6O6Oho9OvXD02bNoW9vT3+9a9/ISwsDAUFBdX2aWlpqfG1bNmyatuuX78eLi4ucHJywpIlS3TW1cLCAmFhYYiPj8eGDRswcuRIvHz50vCTwBLdjR7XSQm5rONlSG6GWLkUXIY0q2MMycpC5CXJFS6yWRfyzd9hhxTSc//+ffTs2RP37t1DWlqaTum5evUqfH190aJFC9SqVQsODg7o3r07EhISqmz3ww8/YMKECfD09ISNjQ0sLS21ThwrVJtOaIciPgISFxeHvLw8zJw5E23btkVeXh6io6PRtWtXnDhxAj5q5RpTAAAbiElEQVQ+Pnjx4gUWLlyIUaNGYfLkyXB0dMTly5exePFiHD58GJcuXao2b8SIESMwe/bsKu+5uLhU+feFCxcQHh6OVatWwc7ODrNmzYKHhwfee+89nXX++OOP0ahRI4wcORIDBw7EoUOHUL9+fX5OCEP0d30pRUaMCIkhyK0bjOnIITnKD5+oHp+UkQ8+MOS7ks+xSyE8AHDt2jUMGjQINjY2+PHHH9GqVSud2xcUFKB58+YYPXo0mjVrhsLCQiQkJGDMmDHIzMxEeHg4gIqlLU6ePIlOnTqhXr16ekfOCtmmE9Wh5GYByc3NhaOjY5X3nj9/jjfeeAPt2rVDSkoKysvLUVBQgAYNGlTZbt++fRgxYgS2bduG0aNHV75vaWmJkJAQrFq1SmfZ0dHRyM7ORnR0NAAgJiYGGRkZej+n5KeffsKQIUPQqFEjJCUlwd3dndHn+EZ/4rM++TE0csM32iSIq8Sxifqwne1XyJsi2xFJbGem5jJiTyjYfrf6zo2mY+M6C7X0USKppOfw4cMICAjAG2+8gePHj6NJkyac99WtWzdkZ2fj/v37AACFQlGZJ7l8+XKEhoYiMzMTzZs3r/ZZsdp04h+oq0tA1KUHAOzt7eHh4VG5WKilpWU16QEALy8vANC4qCgTV23ZsiWOHTuGjIwMZGVlYffu3WjTpg3junfv3h3nz5+HpaUlunTpguTkZMaf5RP9jaKu7i+5SY9yv3ztuy7ku2inPtgOf9e1UKmmFdqNGW3dXbq6IpksAaLp89JGI6WQHoVCgSVLlmDYsGHo378/zp49a5D0AICDgwOsrf85l2wHh4jRphP/QOIjMgUFBbhy5Qratm2rc7tTp04BgMbtEhISUKtWLdja2qJz587YvHlztW38/f3h5eUFd3d3NG/eHA4ODpg0aRKrurZq1Qrp6eno2bMnfH19ER0dzegHyjfM5vEwpps+YPhcM9qiAsYqQoRm+Mi9kl/ulthz8yh58eIFPvzwQ4SHhyMiIgJ79+5F7dq1We9HoVCgtLQUeXl5iI2NxYkTJzBnzhzO9RKrTScqoK4ukQkMDMSePXuQnp6Ojh07atzmr7/+gqenJ1q0aIHz589X+/yQIUPg4uKCnJwcxMfH4/jx45g3bx4+//zzavt69OgRSktL4ezMZqHHqpSVlSE8PBxffvklPvroI8TFxUm2Xo3+ri/gn3WuuCJ2+F+1K4TpEzjbriIuSbNCihOf8xWpI6cuPSVsjld5/XFZwV61y4vNb0Cca16qbi0AePDgAfz9/XHr1i1s2bIFH3zwAed9TZ48GevXrwcAWFlZYfny5Vrn/NHX1SVFm27ukPiISEREBKKiorB69WpMnTpV4zZ///03fHx8kJubi3PnzjFaQd3Pzw8nTpzAX3/9hYYNG/Jc639ISEjA+PHj0aFDBxw4cABOTk6ClaUPYVZ7VyKl+Kii6xiEFh+55b2Yuvhouua4ig/XJVmERUrp+emnn/D+++/DxsYGhw4dQocOHQza34MHD5CXl4fc3FwkJSUhLi4OS5YswX//+99q2+oTH02I1aabK9TVJRKRkZGIiorCkiVLtErPkydP0L9/fzx8+BApKSmMpAcARo8ejZKSkirD5IVg9OjROHv2LB48eIDOnTvj4sWLgpanC2ZLXnBpzOUiPYD27jC20mNu8DHMXwz4XlqjCNyFX7hcH6m6tZRs3LgRPj4+aNWqFS5evGiw9AAVI646deqEQYMGITY2FsHBwYiIiEBeXh4PNRavTTdXSHxEIDIysvI1d+5cjds8efIE/fr1w/3795GSkoJ27dqJXEtmeHl54dKlS3BxcUGvXr2qzV8hJswaVNOY30Q+60+ZImLOaC2E7Mg3n0tK4SktLcXMmTMxfvx4fPTRRzh58qTGASd84OXlhdLSUty7d0+Q/RP8IpcJRkyWRYsWITIyEhEREYiIiNC4jVJ6MjMzkZKSgrfeeotVGVu3bkXNmjXh6enJR5X14uTkhLS0NAQHByMwMBCXL1/GF198IdnqxfpnfQaYzaljLFKh2pVBiyfK9aZfgfKaegb2gvUMmru79B0v1/l5rMHXb0BK4QGAnJwcjB49GmlpaYiJicG0adMEXYYnNTUVVlZWaNmyJS/7E7tNNzdIfAQkOjoaCxYswKBBgzB48GCkp6dX+f+uXbuiqKgIAwcOxNWrV7Fy5Uq8evWqynaOjo5wc3MDUNFXfPv2bfTt2xdNmjRBbm4u4uPjkZKSgsjISDg4OIh2bLa2tti8eTM6duyI0NBQ/PDDD9i1axdef/110eqgjv6JDwF5TSrIR6SByRIUprD0hJzlRh2hBFr+50Bq4QGA77//HoGBgQCA7777Dn369OFt35MmTUK9evXg5eWFxo0b4/Hjx9izZw92796N0NDQyjb48ePHSEtLA1AxSSIAHDt2DA0bNoSjoyPeeecdAPJq080JSm4WEB8fH5w5c0bjEHALCwuUlZUhMzMTbm5usLCw0Ljd2LFjsXHjRgDAkSNHsGzZMty+fRtPnjyBnZ0dOnbsiKlTp2LkyJGCH482Lly4gFGjRuHvv//Ghg0bDBotwRfMRn8BTBdMFAZrsI/YsE14rQH5jegC+J/UTx0xk7nZrI/FlhKOnxPjeKoitfSUlpZi4cKFWLJkCfr27Ytt27YZPD+POps3b8amTZvw22+/4f/+7/9Qu3ZtdOjQARMmTEBAQEDldmlpaVWES7V99/b2rpyuRK5tuqlD4kPwQkFBASZOnIg9e/ZgypQpWLFihWRD3lVhJkD8hfjZoSnqpE+EDB3irK1cVeQwykkduYiPodcJk+NWj+JxjQyKIz9SCw9QMdHrhx9+iHPnzuHzzz/H3LlzYWlJKayEZkh8CN5QKBRYv349pk+fDnd3d+zatUs2M4syjwCJBZOuNk0SxIf4aEK1PiQ+FQghw9qOW9cQdLbnSllvrnNZMTtuOQgPULH0xNixY1GrVi3s3LkTPXv2lLpKhMwhJSZ4w8LCAsHBwTh//jyKi4vRuXNnbNu2TepqAZB+SC031GftFUp6gKqjjcxtBFkpKiRT7GNnOiszm2iUar3ZdqMyO265/JZevXqFWbNmwc/PD927d8fVq1dJeghGUMSHEITCwkKEhIRgy5YtCAoKwurVqzlNDS8E8oj+cEmutlP7kwlsFq9UwibiwzVJXMyID1OJESt52JDvXh19x8Yk6mNcER4AuHv3LkaNGoWrV69i2bJlmD59uqCjtgjTgsSHEJStW7di6tSpcHZ2xpYtW9ClSxepq1SJdALEtyzokgihxYcL1uAmPmy6buR43Er4EB82ESlt58248niAiu70bdu24ZNPPkHDhg2xa9cudO7cWepqEUYGdXURgvLRRx/h0qVLqFOnDrp3747//ve/ePnypdTVAiCfkL3hGNvCpNq61XS95LXIpvgUgb9uOOb7kNNvJDs7G0OHDkVQUBD8/Pxw5coVkh6CExTxIUShtLQUX331FRYuXIiWLVti06ZNsor+AGJFgMTqGlJSBHa5HmKJE9vjETK/SYmcIz4A98kqldJofF1aQEWUZ+vWrZgxYwZsbW0RFxcHPz8/qatFGDEU8SFEwdraGmFhYbh8+TLs7e1lF/0B5PV0yw/KG3mJhhchHWImjatGinQjx+tfGeUZO3YshgwZghs3bpD0EAZDER9CdIwh+qOE/yiQmBEfthGMGhw+wxU5RnwA04n6MD8OuckOQFEeQlgo4kOIjjFEf5Tw+xQsh2UydGFukSBTi4Sxy/OSY4QHoCgPITwU8SEkRTX64+bmhk2bNqFr165SV0srhkWA+HrKZ7ofLtELsbphhI74cE30FrMbisv1oCniw+445Sg7AEV5CPEg8SFkwfXr1zFu3DhcvnwZEyZMwJIlS9CwYUOpq6UVbgIkVPeGpv1y7bIR48bP5zw2ujBF8QG4dknKVXiAit9/SEgITp8+jcDAQHzzzTd47bXXpK4WYaKQ+BCyobS0FOvWrcO8efNgZWWFqKgoTJw4EVZWVlJXTSfsJYjtDY/taB5rkPgApic+7OsmZ9kBKtb4W7hwIWJiYtCyZUvExMRgwIABUleLMHFIfAjZkZubi7lz52LTpk3w9PTEmjVrZJv8rIowUSBDhzCzQaybPolPBUzOg3GulK4PhUKBhIQE/Oc//8HTp08RERGBmTNnwsbGRuqqEWYAiQ8hW86dO4dp06bh559/xvjx47F06VI0atRI6moxgj8JIvGpwBTFB9B8LkxTdpT8+uuvCAkJwdmzZzFixAhER0fDxcVF6moRZgSJDyFrysrKEBcXh/DwcFhYWGDx4sUIDg6WffeXEsOToUl8KhBLfABpoj6mLTtARbfWggULsHr1arRq1QoxMTHo16+f1NUizBASH8IoyMvLQ1hYGOLj49GpUyesXr0a3bp1k7parDBMgtjc+LkOyzY18QGMI+rDHmMSnvLycmzfvh2hoaEoLCzEggULMH36dNSsWVPqqhFmCokPYVScP38e06ZNw+XLl/Hee+8hKioKHh4eUleLNYZPjKhLBMxNfHTVW86RL3YYk+wAFXk8J06cQFhYGK5evYpRo0Zh+fLlaNasmdRVI8wcEh/C6CgrK8OOHTsQERGBBw8eYNy4cViwYIHR5gnwMzu0qhiYoviYYpeffoxNdpScP38ec+fORVpaGnr27IkvvvgCPXr0kLpaBAGAxIcwYoqLixEXF4fFixfj6dOn+OSTTxAWFmbU83/wt0SG8EOluSPmIp3GJT7GKjpKfvvtN4SHh+PAgQNo164dli5dCl9fX1hYWEhdNYKohMSHMHqePXuG6OhoREdHw8rKCqGhoZg+fTrs7e2lrprB8LtWmC7hIPH5B3HFx9hlBwAePHiAyMhIbNq0CS4uLli0aBECAgKMZhACYV6Q+BAmQ25uLqKiorB27Vo4ODhg/vz5mDBhAmrU4NpNIi/4XzAVMHRUkWFlskUs8QGEPB+mIDpK8vPz8cUXXyAmJgZ169bFvHnzEBwcTPPxELKGxIcwOTIzMzF//nxs374dbm5u+OyzzxAYGGhyo0iEESExkKP4kOiwIT8/HzExMfj6669RXl6OOXPmYNasWahTh+2aagQhPiQ+hMly7do1zJ8/HwcPHoSzszNmz56NiRMnmkQXmCaMR4SkEB/xIlqmKDpKsrKysGLFCsTFxUGhUCA4OBhhYWFwdHSUumoEwRgSH8LkuXHjBpYtW4aEhATUr18f06dPR0hICBo0aCB11QRHnjIk7uKcQmLKkqPKnTt3sGzZMmzduhX29vYICQnBp59+ajQzqROEKiQ+hNmQmZmJ5cuXIz4+HtbW1pg8eTJmzpyJpk2bSl010ZGnEMkbc5EcVX7++WcsXboUe/fuRePGjTFr1iwEBwejbt26UleNIDhD4kOYHTk5Ofjmm2+wZs0avHz5EmPHjkVoaChatmwpddUkhWSoAnMUHFUUCgXOnDmDpUuX4sSJE3Bzc0NoaCiCgoJga2srdfUIwmBIfAizpaCgAGvXrsXXX3+Nx48fY/jw4QgJCUGvXr1o3hENmIoYmbvYaKO4uBj79+9HTEwMzp07h/bt22Pu3LkYMWIErK25dk8ShPwg8SHMnqKiImzatAmrVq3C7du30bZtW0ydOhWBgYEU0jcAMUWJZIY7f/75J+Li4rBhwwbk5ubC29sbc+bMweDBg+kBgDBJSHwI4n8oFAqcOnUKsbGxOHToEOzs7DBmzBhMmTIFb775ptTVIwjeKC8vR0pKCmJjY3HkyBHUrl0bQUFBmDJlilGufUcQbCDxIQgNZGVl4dtvv8X69evx6NEj9OrVC1OnTsX7779vcvMBEeZDfn4+Nm/ejLVr1+KPP/5A+/btMW3aNAQEBKB27dpSV48gRIHEhyB0UFJSgoMHDyI2NhZpaWlwdHTExIkTMWHCBLi6ukpdPYLQi0KhwIULF7Bu3TokJiaivLwcI0aMwNSpU9GtWzfqziLMDhIfgmDIzZs3sXbtWmzZsgXPnj1Dz549ERgYiBEjRhj1wqiEafL7778jISEBCQkJyMjIQIsWLTB58mR8/PHHNOEgYdaQ+BAESwoLC3Ho0CFs374d3333HaysrDB48GAEBgZiyJAhNOSXkIy8vDzs2rULCQkJSE9PR+3atTF8+HAEBgbCx8eHFg0lCJD4EIRBPHr0qPJGc/HiRdStWxcffPABAgMD0bt3b1haWkpdRcLEefHiBZKSkrB9+3YkJyfDwsICgwYNQmBgIIYOHYpatWpJXUWCkBUkPgTBE7dv367sWrh79y6aNWuGgIAA/Pvf/0anTp0ol4LgjeLiYqSmpmLnzp3Yv38/CgsL0a1bNwQGBmLkyJFo2LCh1FUkCNlC4kMQPKNQKJCeno7t27dj165dyM/PR7NmzeDn54ehQ4fCx8eHusMI1uTn5+P48eNISkpCcnIynj17htatWyMwMBABAQFmP/M4QTCFxIcgBKSkpARnz55FUlISkpKScO/ePdjb22PgwIHw8/ODr68vPZ0TWsnIyKi8dn744QeUl5fDy8urUqLbt29PkUSCYAmJD0GIhEKhwM2bNytvZOfPn4eFhQW6d+8OPz8/+Pn5wd3dXepqEhJSVlaG9PT0ymvk1q1bsLW1Rb9+/SpF2RwX1SUIPiHxIQiJyMnJwdGjR5GUlITvvvsORUVFaNmyJXx8fODj4wNvb2+6yZk4CoUCv/32G1JTU5GWlobU1FTk5+fD0dERQ4cOxdChQ9GvXz/Y29tLXVWCMBlIfAhCBhQVFeHkyZNITk5GWloabty4AQBo1apVpQR5e3vDyclJ4poShqBQKHDr1q1KyUlLS0NeXh5q1KiBLl26wNvbG76+vnj77bdpRCBBCASJD0HIkNzcXJw5c6by5njz5k0AgLu7e6UE9e7dm0RI5igUCty+fRtpaWmVr5ycHFhbW+Ptt9+ulNpu3bpRVIcgRILEhyCMgJycHJw+fboyUnDr1i0AgIuLCzw9PdG5c2d4enrC09MTjRo1kri25olCocC9e/dw+fJlXL58GZcuXcKVK1fw5MkTWFtbw8vLq1Jae/ToQaJDEBJB4kMQRsijR49w5swZXLp0qfJGW1BQAOAfGVJ90RIF/KKUHNXzr5QcAHB2dq4892+//TZ69OhBi4AShEwg8SEIE0ChUODu3btVbsTqMvTWW2/B3d0drVu3RuvWreHu7o4mTZrQcGgdlJaW4v79+7hz5w5u376NO3fu4NatW7h69apGyVG+GjduLHHNCYLQBokPQZgoShlSdrtcv34dd+7cwb1791BeXg4AqF27dhURUv69devWqFu3rsRHIA4KhQK5ublV5Eb5+v3331FSUgIAsLGxQatWreDu7o633nqLJIcgjBQSH4IwM4qLi3H37t0qN3jlDT8nJ6dyu3r16qFp06ZwcnKq/FPT3+Waq6JQKPDkyRM8fPgQDx8+RHZ2drW/K/8sKioCAFhYWMDV1VWjDLq4uNBIK4IwAUh8CIKopKCgABkZGbhz5w6ysrKqCUJ2dnalJCipU6cOnJycUK9ePdSuXRt16tSpfGn7t729PaytrWFpaVnlZWFhgfLy8mqvkpISFBYW4tmzZ3j27FmVv2v6t1J4iouLq9S1fv36VeRN+WeLFi3g7u4ONzc3Wk6EIEwcEh+CIBijUCjw9OnTahGUhw8f4unTp1rlpLCwEIWFhbzUwdLSUqdU1alTp1Jw1CNUdnZ2vNSBIAjjhcSHIAhRKC8vx/PnzyslSDWiU1ZWVvl3KyurKlEgKysrWFlZVcqNnZ0dJWQTBMEZEh+CIAiCIMwGytQjCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCIIgCMJsIPEhCMJgTp48iaCgILRu3Rr29vZwdnaGv78/rly5UmU7S0tLrS8PD49q+42JiUGbNm1ga2sLNzc3fP755ygtLa223fr16+Hi4gInJycsWbJEsOMkCML4sZa6AgRBGD9xcXHIy8vDzJkz0bZtW+Tl5SE6Ohpdu3bFiRMn4OPjAwBIT0+v9tn09HTMmDED77//fpX3o6KiMH/+fISFhWHAgAG4cOEC5s2bh7/++gtxcXGV2124cAHh4eFYtWoV7OzsMGvWLHh4eOC9994T9qAJgjBKLBQKhULqShAEYdzk5ubC0dGxynvPnz/HG2+8gXbt2iElJUXrZ8eNG4etW7ciIyMDbm5uAID8/Hw4Oztj7NixWLt2beW2S5cuxbx583D9+vXKCFF0dDSys7MRHR0NoCJKlJGRgVWrVvF9mARBmADU1UUQhMGoSw8A2Nvbw8PDA1lZWVo/9+zZM+zZswfe3t6V0gMAycnJKC4uxrhx46psP27cOCgUChw8eLDyvZYtW+LYsWPIyMhAVlYWdu/ejTZt2vBwVARBmCIkPgRBCEJBQQGuXLmCtm3bat0mMTERL168wIQJE6q8f/36dQDAm2++WeX9Jk2aoGHDhrhx40ble/7+/vDy8oK7uzuaN28OBwcHTJo0iccjIQjClCDxIQhCEKZNm4aioiKEh4dr3SY+Ph4NGjTA8OHDq7yfn58PGxsb2NnZVftMgwYNkJ+fX+W9rVu3Ijs7G3/++ScOHjwIa2tKXyQIQjPUOhAEwTsRERHYsWMHVq9ejY4dO2rc5saNG7hw4QJCQkJQs2ZNg8ts0qSJwfsgCML0oYgPQRC8EhkZiaioKCxZsgRTp07Vul18fDwAVOvmAgAHBwcUFxfj5cuX1f7v77//hoODA38VJgjCrCDxIQiCNyIjIytfc+fO1brdq1evsG3bNnTu3Bnt27ev9v/K93799dcq7z969Aj5+flo164dvxUnCMJsIPEhCIIXFi1ahMjISERERCAiIkLntklJScjPz8f48eM1/v+gQYNga2uLzZs3V3l/8+bNsLCwgL+/P1/VJgjCzKAcH4IgDCY6OhoLFizAoEGDMHjw4GoTFXbt2rXKv+Pj42Fvb4+AgACN+2vQoAHmzZuHiIgIvPbaa+jfvz8uXryIyMhITJw4kYarEwTBGZrAkCAIg/Hx8cGZM2egqTmxsLBAWVlZ5b8fPHiA119/HWPGjMGmTZt07jcmJgZr1qxBZmYmnJycMG7cOISHh8PKyor3YyAIwjwg8SEIgiAIwmygHB+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMwGEh+CIAiCIMyG/wd3Z9CL/9osiAAAAABJRU5ErkJggg=="
      ],
      "text/plain": [
       "Figure(PyObject <matplotlib.figure.Figure object at 0x7fa9ba2f0b90>)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "PyObject <matplotlib.collections.PolyCollection object at 0x7fa9b9744350>"
      ]
     },
     "execution_count": 39,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "L= 1.0 # length of the domain\n",
    "Nx= 30 # number of cells\n",
    "Ntheta= 5*20\n",
    "m= createMeshRadial2D(Nx, Ntheta, L, 2π) # create the domain and mesh\n",
    "m.cellcenters.x+=0.1\n",
    "m.facecenters.x+=0.1\n",
    "BC= createBC(m) # create the boundary condition\n",
    "BC.left.a[:]= 0.0 # left side Neumann term equal to zero\n",
    "BC.left.b[:]= 1.0 # left side Dirichlet term coefficient equal to one\n",
    "BC.left.c[:]= 0.0 # left side Dirichlet term value equal to zero\n",
    "BC.left.c[1:5:Ntheta]= 1.0\n",
    "BC.right.a[:]= 0.0\n",
    "BC.right.b[:]= 1.0\n",
    "BC.right.c[:]= 0.0\n",
    "BC.bottom.periodic=true # periodic boundary condition\n",
    "D_val= 1.0 # value of the transfer coefficient\n",
    "D= createCellVariable(m, D_val) # assign D_val to all the cells\n",
    "# harmonic average of the transfer coefficient on the cell faces:\n",
    "D_face= harmonicMean(D)\n",
    "D_face.yvalue[:]=0.003*D_val\n",
    "(Mbc, RHSbc)= boundaryConditionTerm(BC)\n",
    "Mdiff= diffusionTerm(D_face)\n",
    "M= -Mdiff+Mbc # matrix of coefficients\n",
    "RHS= RHSbc # off course!\n",
    "phi= solveLinearPDE(m, M, RHS) # solve the linear PDE\n",
    "visualizeCells(phi) # visualize the results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Not really bad!\n",
    "That's it for now. I will come back with more real life examples."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 0.3.6",
   "language": "julia",
   "name": "julia 0.3"
  },
  "language_info": {
   "name": "julia",
   "version": "0.3.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}
